The Influence of Biopolymer Coating Based on Pumpkin Oil Cake Activated with Mentha piperita Essential Oil on the Quality and Shelf-Life of Grape

Abstract

This work aimed to determine the influence of biopolymer coatings based on pumpkin oil cake, with and without the addition of Mentha piperita essential oil, on the quality and shelf-life of the Afus Ali variety of grapes, stored at room temperature and in the refrigerator. Furthermore, a 10% (w/w) aqueous solution of composite pumpkin oil cake (PuOC) with the addition of 30% glycerol was prepared at 60 °C and pH 10. The active biopolymer coating was prepared similarly by adding 1% (v/v) Mentha piperita essential oil. The quality of packed grapes was tested by determining the dry matter content, total sugar content, total acidity, alcohol content, total phenolic compounds content, and total flavonoid content, as well as by determining the antioxidant activity, through the application of the DPPH, FRAP and ABTS tests. Additionally, microbiological parameters were investigated: total aerobic microbial count, yeasts, and molds. The obtained results proved that in all tested samples, over a certain period of time, the content of dry matter, content of phenolic and flavonoids substances and sugar content decreased as a consequence of the spoilage of grapes, that is, the consumption of sugar for the production of alcohol, which consequently leads to the total acidity increasing. The application of lower storage temperatures and active coating (with Mentha piperita essential oil) had a positive effect on all inevitable reactions. Grapes’ antioxidant potential may be enhanced or maintained by applying PuOC coating with or without Mentha piperita essential oil, which is best observed in the case of the DPPH test. The uncoated sample stored at room temperature had the largest decrease in DPPH values during storage, with changes ranging from 2.119 mg/g to 1.471 μmol mg/g. The samples, coated with PuOC and PuOC with the addition of essential oil, had uniform DPPH values throughout the entire storage period. Additionally, regarding phenolic content, at the end of storage period the highest phenolic content was observed in samples with active coating stored at room temperature (734.746 ± 2.462) and at refrigerator temperature (680.827 ± 0.448) compared with untreated samples and with samples with plain PuOC coating. The presence of active essential oil in the applied coating significantly affected the microbiological profile of grapes during the storage period. Besides the positive impact of the applied lower storage temperature, the effectiveness of the applied active packaging is even greater (microbiological results were in the order of PuOC+essential oil < PuOC < Control). The developed artificial neural networks were found to be adequate for modeling the microbiological profile, antioxidant activity, phenolic and flavonoid content.

1. Introduction

The focus of researchers on creating natural bioplastics from biodegradable biopolymers for food packaging applications has expanded recently. Fossil fuels are a finite resource, and it is essential to discover not only new energy sources but also creative materials to replace plastics made from petroleum [1]. Additionally, the disposal of packaging made of petroleum-based materials has resulted in enormous waste issues and subsequent environmental degradation [2]. The biobased and biodegradable natural polymers are the most intriguing sustainable possible replacements for the fossil sources of plastic products [1,3,4]. Because of their accessibility, environmental friendliness, and capacity to disintegrate via direct consumption, biodegradable polymers, made from natural resources, are regarded as promising alternatives to non-biodegradable synthetic polymers [5].

Even if there is plenty of novel research on this topic, the number of bioplastics generated today is still a small portion compared to oil-based products [2]. On the other hand, the agricultural industry produces a large number of different by-products with biomacromolecules, such as proteins and polysaccharides, so biopolymers from agricultural sources are an intriguing option for producing biodegradable/edible polymers [3].

Oilseed crops, regarded as foods high in energy, are farmed worldwide primarily for edible oil production. Oilseed crops are rich in fibers, antioxidants, vitamins (vitamin E, niacin, and folate), minerals (phosphorus, iron, and magnesium), and monounsaturated and polyunsaturated fatty acids [3]. One of the main byproducts of the oil extraction from oilseeds is seed oil cake, accounting for about 50% of the original seed total weight. Traditionally, either a screw press or solvents are used to extract the oil from oilseeds. The co-product that comes directly from the expeller is referred to as “cake,” but the co-product that has gone through an extra, mostly organic solvent-based de-oiling process is referred to as “meal.” However, there is some ambiguity in how both names are used [6].

According to USDA [7], in 2021/22 world production of major oilseeds (copra, cottonseed, palm kernel, peanut, rapeseed, soybean and sunflower seed) was 604.61 million metric tons. However, despite the enormous amount of agricultural biomass produced worldwide, only a small portion of it is currently used for applications other than those related to human nutrition or animal feed.

The aim of this work was to contribute to the application of biopolymer coatings based on oil seed cakes for food packaging. So far, various biopolymers have been applied on grapes with the aim of testing their efficacy, keeping up their quality and/or prolonging the shelf-life of grapes. These uses include: coatings based on starch/gelatin [8], chitosan and polyvinyl alcohol blending with salicylic acid [9], multilayer films composed of chitosan, sodium alginate and carboxymethyl chitosan-ZnO nanoparticles [10], edible coatings composed of alginate, galactomannans, cashew gum, and gelatin [11], starch-based films reinforced with cellulosic nanocrystals and essential oil [12], natamycin-incorporated nano-TiO₂/poly(butylene adipate-co-terephthalate) (PBAT)/poly(lactic acid) (PLA) biodegradable active film [13], poly(lactic acid) nanofiber packaging containing essential oils [14], epigallocatechin gallate grafted with pectin [15], etc. The novelty of this work is the application of composite biopolymers obtained from waste. These products are also applied as active packaging with the aim to keep up the quality of grapes and prolong their shelf-life.

After cold-pressing the oil from pumpkin seeds (Cucurbita pepo L.), a byproduct known as pumpkin oil cake (PuOC) is produced. Unfortunately, the majority of this cake is either consumed as animal feed or exposed. The purest form of pumpkin seed protein, PuOC, contains 63% proteins, 12% carbohydrates, 4.5% crude fiber, 8.4% oils, and 13% of other components [16]. Pumpkin oil cake was examined by Popovic et al. [17] and Hromiš et al. [18] where the influence of process parameters on the properties of the obtained biopolymer films was examined. The possibility of forming pouches was the next area of investigation into potential uses for PuOC films. Because PuOC films do not have the ability to seal in heat, they were laminated using zein, a material that possesses the aforementioned properties [19]. Earlier work [20] pointed at the antioxidant activity of PuOC films without incorporating any active components. It was further proven that such films can be used as active packaging material and that the addition of essential oils could improve their properties [21,22]. Since it has been proven that pouches made of pumpkin oil cake-based films can preserve modified atmosphere conditions [23], they have been used for flaxseed oil packaging. It was proven that pumpkin oil cake-based packaging could ensure good oxidative stability without inducing significant changes in oil composition [24].

In this work, the biopolymer coating based on the pumpkin oil cake, native and activated with Mentha piperita essential oil, was applied by the immersion method to grapes of the Afus Ali variety and the quality of the grapes was monitored during the storage under certain conditions. Additionally, the principal purpose of this study was to investigate the possibility of anticipating the microbiological features (number of aerobic bacteria—TNAB, number of yeasts and molds) and antioxidant parameters (DPPH, FRAP, ABTS, total phenols content—TPC—and total flavonoids content—TFC) of a coating according to time, data of treatment and spread type (introduced as categorical variables), thus developing an artificial neural network model (ANN).

2. Materials and Methods

2.1. Experimental Material

The following materials were procured for the experiment: umpkin oil cake (PuOC) obtained after cold pressing of the oil (Linum, Serbia), glycerol (>99%) (Fisher Chemicals, USA), demineralized (distilled) water (Alfapanon, Serbia), NaOH (Lach-Ner, Czech Republic), Mentha piperita oil extract (Kirka Corporation, Serbia), “Afus ali” grapes purchased from a local store one day before packing. The grapes were harvested in October 2022 from the local vineyard (Južnobački district, Serbia) one day prior to purchase. Damaged grains were discarded. Only grains of similar size, color, and shape were considered for this study.

2.2. Synthesis of a Biopolymer Coating Based on Pumpkin Seed Oil Cake (PuOC)

A 10% (w/w) PuOC aqueous solution was prepared by using distilled water and grounded PuOC, followed by the addition of 0.3 g glycerol/g PuOC. The pH value of the film-forming solution was adjusted to pH 10 by the gradual addition of 50% aqueous NaOH solution on a magnetic stirrer (IKA, Germany). After achieving a pH of 10, the film-forming solution was heated for 20 min in a water bath at 60 °C and finally filtered through nylon mesh. The active biopolymer coating was prepared in the same way with the addition of 1% (v/v) Mentha piperita essential oil at the very end, so that the obtained biopolymer coating was cooled to room temperature. The active biopolymer coating was homogenized using a homogenizer at 166.67 Hz for 1 min (The SilentCrusher M Homogenizer, Heidolph Instruments, Germany).

2.3. Preparation, Packaging and Storage of Grape Samples with and without a Biopolymer Coating Based on Oil Pumpkin Cake

Grapes of the “Afus Ali” variety were immersed in the prepared solutions: PuOC and PuOC with the addition of Mentha piperita essential oil. After soaking for two minutes, the grapes were left to drain off the excess biopolymer coating by using a sieve and then spread on plastic trays and left in room conditions for 2 h to dry. Grapes without an applied coating were used as a control sample. Grape samples were packed into polystyrene trays measuring 17.5 cm × 13 cm and covered with polyethylene stretch film. Half of the samples were stored at room temperature (23 °C), and the other half at refrigerator temperature (4 °C). The sampling of the dynamics of grape samples stored at room temperature was performed on the 2nd, 4th, 6th, 8th and 10th days, and the samples stored at refrigerator temperature were sampled on the 3rd, 6th, 9th, 12th and 15th days. Different sampling days were chosen based on the assumption that the shelf life of samples stored at refrigerator temperature will be longer, and that therefore it is not necessary to sample these samples as often as samples stored at room temperature. In this way, the experiment itself would be simplified.

Table 1. Sample labeling and experimental design.

Label	Control Grape Sample (Uncoated)	Grape Sample with Biopolymer Coating	Grape Sample with Active Biopolymer Coating (Mentha Piperita Essential Oil)	Storage Period (Days)									Storage Temperature (°C)
				2	3	4	6	8	9	10	12	15	4	23
S2	+			+										+
SB2		+		+										+
SU2			+	+										+
S4	+					+								+
SB4		+				+								+
SU4			+			+								+
S6	+						+							+
SB6		+					+							+
SU6			+				+							+
S8	+							+						+
SB8		+						+						+
SU8			+					+						+
S10	+									+				+
SB10		+								+				+
SU10			+							+				+
F3	+				+								+
FB3		+			+								+
FU3			+		+								+
F6	+						+						+
FB6		+					+						+
FU6			+				+						+
F9	+								+				+
FB9		+							+				+
FU9			+						+				+
F12	+										+		+
FB12		+									+		+
FU12			+								+		+
F15	+											+	+
FB15		+										+	+
FU15			+									+	+

S—control grape sample, stored at room temperature; SB—grape sample with biopolymer coating, stored at room temperature; SU—grape sample with active biopolymer coating (with the addition of essential oil), stored at room temperature; F—control grape sample stored at refrigerator temperature; FB—grape sample with biopolymer coating, stored at refrigerator temperature; FU—grape sample with active biopolymer coating (with the addition of essential oil), stored at refrigerator temperature.

explains the sample designations used, while Figure 1 shows the appearance of the samples at the beginning of storage.

2.4. Examination of the Quality of Packaged Grapes

2.4.1. Preparation of Liquid Extract

The grape samples were turned into a mash using a stick mixer. About 5 g of the prepared sample was weighed into an Erlenmeyer flask and 50 mL of 95% methanol was poured as an extraction agent. The flasks were covered and placed on a laboratory mixer (Unimax 1010, Heidolph Instruments GmbH & CoKG, Germany) for 24 h in a dark place. After extraction, the samples were filtered into vials and stored in a refrigerator until analysis. Methanolic extracts were used to determine the content of total phenols, total flavonoids and DPPH.

2.4.2. Dry Matter Content

The dry matter content was determined gravimetrically by drying the samples at a temperature of 105 ± 0.5 °C to a constant mass. The percentage of dry matter content is equal to:

$$DM=\frac{M_2-M_0}{M_1-M_2}\cdot 100(\%)$$

(1)

wherein:

M₀—mass of vessel and auxiliary material (filter paper, sand, glass rod, lid), in g;

M₁—mass of the same container with the tested sample before drying, in g;

M₂—mass of the same vessel with the tested sample after drying, in g.

2.4.3. Content of Total Sugars

Total sugars were determined by the Luff–Schoorl method [25]. The method is based on the principle that, under certain conditions, reducing sugars (natural inversion) will convert copper sulfate (CuSO₄) from Luff’s solution into copper oxide (Cu₂O). The unused amount of cupric ions is retitrated with a thiosulfate solution. The amount of sugar from the table is read from the difference between the consumption for the blank and the test. The non-reducing disaccharide (sucrose) must first be inverted, that is, the reducing monosaccharides must be hydrolyzed with an acid, and then determined using Luff’s solution. The difference between the obtained total invert and the natural invert gives the amount of reducing sugars formed by sucrose inversion.

2.4.4. Total Acidity

The total (titratable) acidity of the sample was determined according to ISO 750:1998. It is a volumetric method, using a sodium hydroxide (NaOH) standard solution, with phenolphthalein as an indicator. Total acidity was expressed in g of tartaric acid per 100 g sample [25].

2.4.5. Alcohol Content

Ethanol, separated by distillation, is oxidized with potassium bichromate in the presence of sulfuric acid, and excess potassium bichromate is retitrated with ammonium ferrous sulfate in the presence of the iron-ortho-phenanthroline indicator. The method is applied for the determination of ethanol in fruit and vegetable products where the amount of ethanol does not exceed 5% (m/m).

2.4.6. Content of Total Phenolic Compounds

The content of total phenols in liquid methanolic extracts was determined by the spectrophotometric method according to Folin–Ciocalteu [26]. The reaction mixture for determining the content of total phenols in the sample was prepared by mixing 2 mL of the sample, 2.5 mL of the Folin–Ciocalteu reagent, and 7.5 mL of Na₂CO₃ in a measuring vessel of 50 mL. The measuring vessel was filled up to the mark with distilled water. Gallic acid was used as a standard. Absorbances were measured at 750 nm. Based on the measured absorbance, the concentration (mg/mL) of phenolic compounds was read from the calibration curve of the standard gallic acid solution, and then the content of total phenolic compounds in the sample was expressed as the gallic acid equivalent (mg GAE/100 g).

2.4.7. Content of Total Flavonoids

The content of total flavonoids in extracts of fresh and dried grapes was determined by the colorimetric method, in accordance with Markham [27]. The reaction mixture was prepared by mixing 1 mL of the extract with 4 mL of the distilled water and 0.3 mL of 5% the NaNO₂ solution. The mixture was then incubated for five minutes at room temperature, after which 0.3 mL of a 10% AlCl₃ × 6H₂O solution was added. After six minutes, when the solution became intensely yellow, 2 mL of NaOH solution was added. The reaction mixture was supplemented with distilled water up to 10 mL and the absorbance was measured at 510 nm. The content of total flavonoids is expressed in catechin equivalents per unit mass of the sample (mg CAE/100 g). A standard catechin solution was used to create the calibration curve.

2.5. Determination of Antioxidant Activity

2.5.1. 2,2-Diphenyl-1-Picryl-Hydrazyl-Hydrate Assay (DPPH)

The DPPH assay was performed using a modified method originally presented in the study by Brand-Williams et al. [28]. A methanolic solution of the DPPH reagent (65 μM) was adjusted by adding methanol to obtain an absorbance of 0.70 (±0.01). The previously prepared extract and DPPH reagent were mixed (0.1 mL + 2.9 mL) in 10 mm glass cuvettes and incubated at room temperature for 60 min. Measurement of the neutralization of free radicals was carried out at 517 nm. A UV/VIS spectrophotometer (LLG-uniSPEC 2 Spectrophotometer) was used for the spectrophotometric measurement of absorbance. Results are expressed as mg Trolox equivalents per g (mg Trolox/g).

2.5.2. Ferric Reducing Antioxidant Power Assay (FRAP)

The FRAP test was performed by a modified method originally presented in the study by Benzie and Strain [29]. The FRAP reagent was prepared from 300 mM acetic buffer (pH 3.6), 10 mM 2,4,6-tripyridyl-s-triazine (TPCTZ) in 40 mmol/L HCl solution and 20 mM/L FeCl₃ aqueous solution in the ratio 10:1:1 (v/v/v). Previously prepared extracts and FRAP reagent were mixed (0.1 + 2.9 mL) and incubated in the dark at a temperature of 37 °C for 10 min. After incubation, the absorbance of the sample was measured at 593 nm. The results are expressed as mg equivalent of Fe²⁺ ions per g (mg Fe²⁺/g).

2.5.3. ABTS Radical Scavenging Assay (ABTS)

The ABTS test was performed by a modified method described in the study by Re et al. [30]. The ABTS reagent solution was prepared by mixing a 7 mM aqueous solution of 2,2’-azino-bis(3-ethylbenzothiazoline-6-sulfonic acid) diammonium salt (ABTS) and 2.45 mM potassium persulfate in a 1:1 ratio (v/v) and incubating them in the dark at room temperature for 16 h. The ABTS reagent was diluted with 300 mM acetate buffer (pH 3.6) to adjust the absorbance at 0.734 to 0.70 (±0.01). The previously prepared extract and ABTS reagent were mixed (0.1 + 2.9 mL) and incubated in the dark at room temperature for 5 h. After incubation, the absorbance of the sample was measured at 734 nm. A UV/VIS spectrophotometer (LLG-uniSPEC 2 Spectrophotometer) was used for spectrophotometric measurement of absorbance. Results are expressed as mg Trolox equivalents per g (mg Trolox/g).

All measurements were performed in three replicates.

2.6. Microbiological Examination

The preparation of samples for microbiological testing was carried out in accordance with standard ISO 6887-1:2017 [31]. Standard methods were used to enumerate selected groups of microorganisms: total aerobic microbial count [32], total molds and yeasts count [33]. The diluent (buffered peptone water) and all culture media were acquired from HiMedia (Mumbai, India). The number of microorganisms in the tested samples is expressed as colony forming units (cfu) per gram.

The number of microorganisms in the tested samples is expressed as the mean of three measurements with standard deviation.

2.7. Statistical Analyses

Statistical processing of the data was presented by the STATISTICA 10.0 software package (StatSoft Inc., Tulsa, OK, USA). The results were exhibited as mean ± standard deviation of triplicate analyses for all measurements.

Principal component analysis (PCA) was applied to reveal the feasible correlations between measured parameters and employed to categorize objects.

2.7.1. Artificial Neural Network (ANN) Modeling

A multi-layer perceptron model (MLP), which consisted of three layers (input, hidden and output layers), was engaged in the establishment of a model building. Prior to the computation, experimental data were normalized to enhance the behaviour of the ANN. The Broyden–Fletcher–Goldfarb–Shanno (BFGS) algorithm was employed as an iterative method for solving unconstrained nonlinear optimization problems in ANN modeling [34].

The experimental database for ANN was randomly divided into training, cross-validation and testing data (with 60%, 20% and 20% of experimental data, respectively). The training data set was used for the learning cycle of ANN and also for the evaluation of the optimal number of neurons in the hidden layer and the weight coefficient of each neuron in the network [35].

Coefficients associated with the hidden layer were grouped in matrices W₁ and B₁, while coefficients associated with the output layer were grouped in matrices W₂ and B₂. The neural network is usually presented using matrix notation (Y is the matrix of the output variables, f₁ and f₂ are transfer functions in the hidden and output layers, respectively, and X is the matrix of input variables) [34]:

$$Y=f_1(W_2\cdot f_2(W_1\cdot X+B_1)+B_2)$$

(2)

The elements of matrices W₁ and W₂ were determined during the ANN learning cycle and during the iterative procedure, with an optimization algorithm being used to minimize the error between network outputs and experimental results. The coefficients of determination were used as parameters to check the performance of the obtained ANN model [36].

2.7.2. Global Sensitivity Analysis

Yoon’s global sensitivity equation for the obtained ANN model was exploited to estimate the relative impact of the input parameters (time, temperature and applied coating) on output variables (microbiological profile, antioxidant activity, phenolic and flavonoid content), depending on the designed ANN model weight coefficients [37]:

$$R{I}_{ij}(\%)=\frac{{\displaystyle \sum _{k=0}^n(w_{ik}\cdot w_{kj})}}{{\displaystyle \sum _{i=0}^m\left|{\displaystyle \sum _{k=0}^n(w_{ik}\cdot w_{kj})}\right|}}\cdot 100\%$$

(3)

where: w—weight coefficient in ANN model, i—input variable, j—output variable, k—hidden neuron, n—number of hidden neurons, m—number of inputs.

2.7.3. The Accuracy of the Model

The numerical confirmation of the obtained ANN and RFR models was performed using statistical tests, such as coefficient of determination (r²), reduced chi-square (χ²), mean bias error (MBE), root mean square error (RMSE) and mean percentage error (MPE) methods. These commonly used parameters were calculated according to Puntarić et al., [38]:

$${\chi }^2=\frac{{\displaystyle \sum _{i=1}^N{(x_{\exp ,i}-x_{pre,i})}^2}}{N-n}$$

(4)

$$RMSE={\left[\frac{1}{N}\cdot {\displaystyle \sum _{i=1}^N{(x_{pre,i}-x_{\exp ,i})}^2}\right]}^{1/2}$$

(5)

$$MBE=\frac{1}{N}\cdot {\displaystyle \sum _{i=1}^N(x_{pre,i}-x_{\exp ,i})}$$

(6)

$$MPE=\frac{100}{N}\cdot {\displaystyle \sum _{i=1}^N(\frac{\left|x_{pre,i}-x_{\exp ,i}\right|}{x_{\exp ,i}})}$$

(7)

where x_exp,i were experimental values and x_pre,i were the model predicted values, N and n are the number of observations and constants, accordingly.

3. Results and Discussion

3.1. Dry Matter Content

In all tested samples, over time, the content of dry matter decreased (

Table 2. The change in the dry matter content, total sugars/alcohol content and total acidity of packed grapes.

Sample	DM	Sugar	TA	ALC
0	20.440 ± 0.025 h	19.678	0.568	0.124
S2	19.690 ± 0.248 fgh	12.716	0.453	0.648
SB2	18.666 ± 0.089 bcdef	14.249	0.553	0.316
SU2	19.434 ± 0.346 efgh	14.150	0.571	0.303
S4	18.712 ± 0.353 bcdef	15.944	0.700
SB4	18.232 ± 0.760 abc	15.938	0.801
SU4	18.471 ± 0.533 bcde	17.648	0.786
S6	19.312 ± 0.175 defg	15.640	0.723	0.114
SB6	19.112 ± 0.142 cdefg	16.287	1.041	0.194
SU6	19.079 ± 0.374 cdefg	17.497	0.970	0.093
S8	18.722 ± 0.143 bcdef	14.515	0.943
SB8	18.423 ± 0.329 bcde	13.503	1.021
SU8	18.675 ± 0.066 bcdef	17.969	1.228
S10	18.286 ± 0.169 bcd	17.142	1.197	0.040
SB10	17.169 ± 0.190 a	10.698	1.382	0.031
SU10	18.003 ± 0.048 ab	15.667	1.401	0.688
F3	19.189 ± 0.201 cdefg	12.707	0.452	1.684
FB3	19.072 ± 0.400 cdefg	17.332	0.467	0.054
FU3	19.255 ± 0.118 cdefg	12.141	0.480	0.685
F6	19.315 ± 0.205 defg	12.040	0.663
FB6	20.096 ± 0.142 gh	15.244	0.842
FU6	19.598 ± 0.223 fgh	14.057	0.993
F9	19.640 ± 0.258 fgh	16.980	0.736	0.644
FB9	19.382 ± 0.141 efgh	13.828	0.661	0.654
FU9	19.887 ± 0.081 gh	19.197	0.700	1.354
F12	20.120 ± 0.034 gh	17.572	0.671	0.720
FB12	19.593 ± 0.273 fgh	16.948	0.692	0.725
FU12	19.932 ± 0.023 gh	13.566	1.937	0.732
F15	19.730 ± 0.048 fgh	8.945	1.495
FB15	19.405 ± 0.007 efgh	16.608	2.014
FU15	19.790 ± 0.015 gh	11.977	1.969

Means in the same column with different superscript are statistically different (p ≤ 0.05).

) from an initial 20.440% to 17.169% for SB10 sample; 18.003% for SU10; 19.405% for FB15 and 19.790% for FU15 sample A greater reduction or loss in dry matter was observed in grapes stored at room temperature, while there were smaller changes in grapes stored at refrigerator (4 °C), which is in accordance with the findings of de Souza et al. [11]. According to Melo et al. [39], water evaporation in the fruit, which can be very large at higher temperatures during storage, is linked to weight loss. Additionally, the coating’s porosity can let water evaporate and lead to weight loss [40], as can the fact that biopolymer structures tend to swell in the presence of humidity, with all these factors causing poor resistance to moisture and water loss [41]. On the other hand, there were higher values of dry matter observed in samples with an active PuOC coating compared to the applied biopolymer itself, which is a consequence of the hydrophobic structure of the added essential oil.

3.2. Content of Total Sugars/Alcohol Content and Total Acidity

It is noticeable that the sugar content changes during storage in both temperature regimes. The decrease in sugar content compared to day 0 is a consequence of the spoilage of grapes, that is, the consumption of sugar for the production of alcohol (Table 2).

Table 2 shows data related to the change in total acidity content, which is expressed in relation to the dry matter of the sample. A trend of increasing total acidity can be observed in samples stored at room temperature (23 °C), as well as those stored in a refrigerator (4 °C). The increase in total acidity was expected due to the grape spoilage during storage. The application of lower storage temperatures contributed to better values of total acidity of the grape samples. Previous research has demonstrated that organic acids and soluble solids are used by respiratory mechanisms to preserve the fruit’s normal activity during storage [42]. Cold storage is considered an additional process that may slowdown the physiological processes in grapes, which in turn may slow the consumption of organic acids [39].

Higher values of total acidity are observed in the sample with biopolymer coating (both with and without the addition of essential oil) compared to the untreated samples in both temperature regimes. On the 10th day of storage for samples stored at room temperature, total acidity was 1.197%, for samples with a coating of 1.382% and for samples with an active coating of 1.401%. The same trend was observed for samples stored at refrigerator temperature: uncoated samples had a total acidity of 1.495%, coated samples 2.014%, and active coated samples 1.969% (whose increase in total acidity value occurred on the 12th day of storage (1.937%)). An increase in total acidity is a consequence of the appearance of spoilage and, given that the grape samples were also visually assessed by the presence of color changes, of spots with molds.

3.3. Determination of Antioxidant Activity

Results of the antioxidant potential of packed grapes are presented in

Table 3. The change of antioxidant activity, total phenolic compounds and total flavonoids of packed grapes.

Sample	DPPH	FRAP	ABTS	TPC	TFC
0	2.268 ± 0.018 p	0.666 ± 0.005 hi	4.227 ± 0.034 r	640.100 ± 1.301 j	1.583 ± 0.005 o
S2	2.119 ± 0.010 m	0.775 ± 0.009 lm	3.473 ± 0.059 m	659.981 ± 1.350 l	1.470 ± 0.015 m
SB2	1.895 ± 0.013 k	0.776 ± 0.004 lm	3.388 ± 0.025 l	683.697 ± 0.475 op	1.450 ± 0.015 m
SU2	1.708 ± 0.008 h	0.655 ± 0.005 gh	3.066 ± 0.008 i	639.622 ± 1.368 j	1.252 ± 0.002 j
S4	1.812 ± 0.010 i	0.807 ± 0.009 no	3.642 ± 0.008 o	689.813 ± 1.421 pq	1.316 ± 0.012 k
SB4	1.862 ± 0.015 jk	0.916 ± 0.005 q	4.015 ± 0.042 p	749.782 ± 1.458 u	1.392 ± 0.015 l
SU4	1.605 ± 0.015 g	0.666 ± 0.013 hi	2.981 ± 0.042 h	698.297 ± 2.879 r	1.005 ± 0.005 de
S6	1.408 ± 0.010 cd	0.513 ± 0.009 a	2.524 ± 0.059 b	560.532 ± 0.918 e	0.897 ± 0.005 ab
SB6	1.441 ± 0.018 de	0.628 ± 0.011 efg	3.219 ± 0.042 k	577.991 ± 2.319 fg	1.054 ± 0.015 fg
SU6	1.320 ± 0.008 a	0.495 ± 0.002 a	2.897 ± 0.042 efg	615.128 ± 1.393 h	0.888 ± 0.015 a
S8	1.532 ± 0.008 f	0.577 ± 0.000 bc	2.820 ± 0.051 cd	650.627 ± 1.420 k	1.032 ± 0.007 ef
SB8	1.292 ± 0.015 a	0.549 ± 0.002 b	2.507 ± 0.042 b	552.976 ± 1.924 d	1.069 ± 0.020 fg
SU8	1.620 ± 0.010 g	0.749 ± 0.013 kl	2.905 ± 0.017 fg	693.767 ± 2.372 qr	0.995 ± 0.005 de
S10	1.471 ± 0.018 e	0.611 ± 0.002 de	2.778 ± 0.059 c	622.801 ± 2.907 i	0.932 ± 0.015 bc
SB10	1.365 ± 0.013 b	0.591 ± 0.007 cd	2.481 ± 0.034 b	575.685 ± 2.064 f	1.005 ± 0.005 de
SU10	1.630 ± 0.015 g	0.689 ± 0.011 ij	3.134 ± 0.042 j	734.746 ± 2.462 t	1.091 ± 0.012 g
F3	2.316 ± 0.005 q	1.083 ± 0.013 t	4.142 ± 0.051 q	721.047 ± 0.924 s	1.832 ± 0.015 q
FB3	2.225 ± 0.015 o	0.956 ± 0.013 r	4.312 ± 0.051 s	726.979 ± 0.465 s	1.744 ± 0.015 p
FU3	2.182 ± 0.013 n	0.834 ± 0.011 o	3.651 ± 0.034 o	657.112 ± 3.683 l	1.538 ± 0.015 n
F6	1.902 ± 0.010 k	0.640 ± 0.005 fgh	2.956 ± 0.051 gh	526.757 ± 3.211 b	1.208 ± 0.012 i
FB6	1.393 ± 0.015 bc	0.686 ± 0.011 ij	2.854 ± 0.017 def	522.188 ± 0.441 b	0.966 ± 0.015 cd
FU6	2.278 ± 0.008 pq	0.793 ± 0.009 mn	3.490 ± 0.025 mn	675.671 ± 0.452 mn	1.450 ± 0.015 m
F9	1.817 ± 0.010 i	0.736 ± 0.007 k	2.905 ± 0.034 fg	534.419 ± 2.256 c	1.250 ± 0.010 j
FB9	1.844 ± 0.013 ij	0.735 ± 0.002 k	2.930 ± 0.025 gh	584.114 ± 1.829 g	1.167 ± 0.015 h
FU9	1.532 ± 0.018 f	0.617 ± 0.007 def	2.320 ± 0.025 a	524.583 ± 0.891 b	0.897 ± 0.005 ab
F12	2.961 ± 0.005 t	1.239 ± 0.009 u	5.820 ± 0.051 t	694.766 ± 0.881 qr	2.219 ± 0.005 r
FB12	1.549 ± 0.010 f	0.698 ± 0.013 j	2.490 ± 0.025 b	493.260 ± 3.167 a	0.966 ± 0.015 cd
FU12	1.948 ± 0.015 l	0.775 ± 0.013 lm	2.837 ± 0.017 cde	620.019 ± 2.223 hi	1.074 ± 0.015 g
F15	2.618 ± 0.015 s	1.001 ± 0.007 s	4.100 ± 0.042 q	582.735 ± 2.695 g	1.707 ± 0.022 p
FB15	2.396 ± 0.015 r	0.936 ± 0.007 qr	3.549 ± 0.051 n	673.958 ± 1.827 m	1.450 ± 0.015 m
FU15	2.094 ± 0.015 m	0.873 ± 0.002 p	3.693 ± 0.042 o	680.827 ± 0.448 no	1.245 ± 0.010 ij

Means in the same column with different superscript are statistically different (p ≤ 0.05).

and they differ in all evaluated methods (ABTS, DPPH, and FRAP). These differences result from the various mechanisms used in radical antioxidant responses [43]. Although these differences are undeniable, for all antioxidant responses it can be noted that the values are higher in grape samples stored at a lower temperature, regardless of whether an (active) coating based on PuOC was applied.

The uncoated samples stored at room temperature had the highest decrease in DPPH values during storage, ranging from 2.119 mg/g to 1.471 μmol mg/g. The samples coated with PuOC and PuOC with the addition of essential oil had uniform DPPH values through the entire storage period. The uncoated grape sample stored at refrigerator temperature and samples on the 12th day (sample labeled as F12) had the highest DPPH value (2.961 mg Trolox/g). The same sample also had the highest ABTS (5.82 mg Trolox/g) and FRAP (1.239 mg Fe²⁺/g) values.

In both groups of samples, there was a slight decrease in capacity for ABTS radical scavenging, although in some samples, the decrease is not linear, which may be a consequence of sample inhomogeneity. The ABTS values on the second and third days of storage are uniform for all samples and range from 3.066 mg Trolox/g to 4.142 mg Trolox/g. Later, during storage, significantly lower values were obtained for samples stored at room temperature, compared to samples stored at refrigerator temperature.

For the FRAP method, the Fe³⁺ reducing power in all treatments undergoes a slight decline until halfway throughout the storage period, followed by an increase in values, so that the final values are the same as the initial values or slightly higher.

The findings indicate that grapes’ antioxidant potential may be enhanced or maintained by applying PuOC coatings, with or without Mentha piperita essential oil. In all of the samples, the coating application of PuOC (with or without the addition of Mentha piperita essential oil) preserved the antioxidant potential of the grapes, which is in the agreement with Tahir et al. [44]. The reason is that biopolymer coatings regulate the ripening process as well as the hydrolysis reactions and reduce changes in phenolic compounds, effects which can have an impact on grapes’ antioxidant potential [45].

3.4. Content of Total Phenolic Compounds and Total Flavonoids

Table 3 shows data related to the change in the composition of total phenolic compounds. The presence of phenolic compounds was found in all tested samples. Phenol substances are positively correlated to the quality of grapes and with the antioxidant activity of grapes. As can be noted from Table 3, the values for the content of phenolic substances in grapes varied during storage with a slight decreasing trend. The underlying mechanism is related to the activity of polyphenol oxidase and peroxidase present in grapes [46]. Polyphenol oxidase can oxidize polyphenols into quinones in an aerobic environment, which has the effect of promoting grape browning and polyphenols content reductions [47]. At the end of storage period, the highest phenolic content was observed for samples with active coating stored at room temperature (734.746 ± 2.462) and at refrigerator temperature (680.827 ± 0.448) compared with untreated samples and with samples with plain PuOC coating.

The presence of flavonoids was also found in all tested samples. Similar to the phenolic content, a discrete decrease in the flavonoid content can be observed during the storage period, which is in agreement with the findings of Lo’ay et al. [9]. In general, higher values were obtained for samples stored at refrigerator temperature compared to samples stored at room temperature. According to the obtained results presented in Table 3, the contribution of PuOC coating to the preservation of flavonoid content is significant compared to that of untreated samples and samples with an active PuOC coating with the addition of Mentha piperita essential oil.

3.5. Microbiological Examination

Table 4. The change of microbiological profile of packed grapes.

Sample	TNAB [×107]	Yeasts [×107]	Molds [×105]
0	0.013 ± 0.004 a	0.002 ± 0.001 a	0.002 ± 0.002 a
S2	1.203 ± 0.454 cdef	0.507 ± 0.210 abcdef	0.048 ± 0.016 a
SB2	0.477 ± 0.087 abcd	0.480 ± 0.288 abcde	0.045 ± 0.010 a
SU2	2.133 ± 0.757 g	0.103 ± 0.032 abc	0.008 ± 0.002 a
S4	0.413 ± 0.101 abcd	0.217 ± 0.076 abc	0.037 ± 0.015 a
SB4	0.433 ± 0.076 abcd	0.533 ± 0.076 abcdef	0.038 ± 0.020 a
SU4	0.163 ± 0.032 ab	0.203 ± 0.127 abc	0.029 ± 0.015 a
S6	0.247 ± 0.050 ab	0.190 ± 0.050 abc	0.056 ± 0.017 a
SB6	0.330 ± 0.125 abc	0.520 ± 0.075 abcdef	0.034 ± 0.006 a
SU6	0.230 ± 0.062 ab	0.213 ± 0.078 abc	0.011 ± 0.001 a
S8	0.217 ± 0.057 ab	0.132 ± 0.010 abc	0.627 ± 0.110 ab
SB8	0.163 ± 0.057 ab	0.303 ± 0.135 abcd	0.320 ± 0.080 a
SU8	0.133 ± 0.031 ab	0.217 ± 0.031 abc	0.055 ± 0.015 a
S10	1.520 ± 0.495 efg	0.150 ± 0.040 abc	0.473 ± 0.110 ab
SB10	1.157 ± 0.319 cdef	0.393 ± 0.179 abcde	0.467 ± 0.125 ab
SU10	0.777 ± 0.493 abcde	0.627 ± 0.583 abcdef	0.120 ± 0.010 a
F3	0.713 ± 0.090 abcde	0.026 ± 0.019 a	0.006 ± 0.003 a
FB3	1.253 ± 0.647 defg	0.220 ± 0.060 abc	0.010 ± 0.001 a
FU3	0.250 ± 0.100 ab	0.002 ± 0.000 a	0.003 ± 0.001 a
F6	0.213 ± 0.099 ab	0.173 ± 0.050 abc	0.017 ± 0.006 a
FB6	1.467 ± 0.569 efg	0.793 ± 0.555 bcdefg	0.016 ± 0.006 a
FU6	0.500 ± 0.100 abcd	0.160 ± 0.061 abc	0.009 ± 0.001 a
F9	0.391 ± 0.494 abcd	0.068 ± 0.014 ab	0.043 ± 0.012 a
FB9	0.157 ± 0.038 ab	1.267 ± 0.473 fg	0.019 ± 0.006 a
FU9	1.800 ± 0.361 fg	1.467 ± 0.723 g	0.010 ± 0.001 a
F12	0.163 ± 0.047 ab	0.127 ± 0.040 abc	0.867 ± 0.084 ab
FB12	1.020 ± 0.164 bcdef	1.090 ± 0.271 efg	0.930 ± 0.200 ab
FU12	1.187 ± 0.359 cdef	1.007 ± 0.261 defg	0.423 ± 0.090 ab
F15	0.053 ± 0.013 a	0.723 ± 0.357 abcdefg	1.800 ± 0.458 bc
FB15	0.540 ± 0.053 abcd	0.833 ± 0.176 cdefg	1.000 ± 0.092 ab
FU15	0.433 ± 0.123 abcd	0.253 ± 0.084 abcd	0.600 ± 0.075 ab
F18	0.063 ± 0.007 a	0.133 ± 0.042 abc	1.300 ± 0.265 ab
FB18	0.390 ± 0.066 abcd	0.530 ± 0.108 abcdef	1.027 ± 0.155 ab
FU18	0.145 ± 0.048 ab	0.283 ± 0.126 abcd	0.683 ± 0.070 ab
F21	0.065 ± 0.023 a	0.113 ± 0.035 abc	3.167 ± 1.795 cd
FB21	0.207 ± 0.116 ab	0.450 ± 0.087 abcde	3.600 ± 1.709 d
FU21	0.237 ± 0.146 ab	0.225 ± 0.158 abc	0.837 ± 0.111 ab

Means in the same column with different superscript are statistically different (p ≤ 0.05).

provides results related to the microbiological profile of grape samples. High initial values for total aerobic microbial count were observed, which affected all other values during storage. Although the obtained results are uneven, it can be noted that the application of an active PuOC coating with the addition of Mentha piperita essential oil is the most important factor for the microbiological stability of grapes packaged at room temperature. SU10 sample had 0.777 × 10⁷ cfu/g, compared with SB10 sample (1.157 × 10⁷ cfu/g) and compared with the uncoated S10 sample (1.52 × 10⁷ cfu/g). The grape quality preservation implies the application of lower temperatures, which was confirmed because the values on 21st day for samples stored at refrigerator temperature were significantly lower (in the range 0.065 × 10⁷–0.237 × 10⁷ cfu/g) than those for samples stored on 10th day at room temperature (0.777 × 10⁷–1.52 × 10⁷ cfu/g).

The results of yeast determination showed that the influence of storage temperature is negligible. On the other hand, the application of an active coating based on PuOC is more significant because, for all tested samples, on each sampling day, lower values were obtained compared to untreated samples, as well as for samples with only a coating based on PuOC. This would mean that the biggest contribution to the low presence of mold is the application of Mentha piperita essential oil. The total number of yeasts in each sample group was in the order of PuOC+essential oil < PuOC < Control.

The same results were obtained when determining molds. In each tested sample group, the lowest values were obtained for samples coated with added essential oil, i.e., for which an active coating was applied. The first significant increase in the mold value of the samples stored at room temperature was observed on the sixth day and for the untreated sample (S6) was 0.627 × 10⁵ cfu/g. In samples stored at refrigerator temperature, the increase in the number of molds was observed at the 12th day, and it was 0.93 × 10⁵ cfu/g for sample FB12. This fact favors the use of lower storage temperatures. By the 21st day of storage, there was a significant increase in the presence of mold, 3.167 × 10⁵ cfu/g for sample F21, 3.6 × 10⁵ cfu/g for sample FB and 0.837 × 10⁵ cfu/g for sample FU21.

Figure 2 displays the appearance of the samples at the end of the experiment (18th and 21st day), when microbiological damage is already visible.

Mentha piperita has various biological activities: antioxidant activities, cytotoxicity activities, anti-inflammatory properties, as well as antimicrobial activities [48]. According to [49], Mentha piperita common major components are menthol (oxygenated monoterpene), menthone (oxygenated monoterpene), carvone (oxygenated monoterpene), anethole (phenylproprenoid), 1,8-cineole (oxygenated monoterpene) and common minor components are menthyl acetate, limonene (monoterpene hydrocarbon), α-pinene (monoterpene hydrocarbon), β-pinene (monoterpene hydrocarbon) and myrcene (monoterpene hydrocarbon).

The structural functional groups of major components play an important role in the biological activity of essential oils. Menthol and menthone are cyclic and oxygenated monoterpenes that play essential roles in the disorganization of cell membrane structures, causing depolarization and physical or chemical alterations, thereby disrupting metabolic activities [50]. These major active components penetrate the cell membrane and target the ergosterol biosynthesis pathway, thus impairing its biosynthesis. Simultaneously, they react with the membrane itself with their reactive hydroxyl moiety, and the extensive lesion on the membrane is a combined effect of the two events [51]. Minor components also significantly influence the antimicrobial properties of the Mentha piperita essential oil through synergistic interactions [52].

According to the available literature, Mentha piperita has very strong antimicrobial potential against various bacteria, yeasts and molds [51,53,54,55]. As such, numerous applications in the food industry have been conducted [50,56,57]. The results of this research support the fact about the antimicrobial effect of the Mentha piperita essential oil when it is incorporated into a biopolymer coating.

3.6. PCA Analysis

The points displayed in the PCA graphic, which are numerically in close vicinity to each other, demonstrate the similarity of patterns that portray these data. The direction of the vector explaining the variable in factor space discloses a rising trend of these variables, and the longitude of the vector is relative to the square of the correlation values among the fitting value for the variable and the variable itself. The angles, amidst corresponding variables, denote the degree of their correlations (minor angles corresponding to elevated correlations).

The PCA of the microbiological data explained that the first two components accounted for 81.63% of the total variance (47.85 and 33.78%, respectively) in the three-variable factor space (microbiological parameters). Considering the mapping of the PCA performed on the data, molds (which contributed 19.4% of the total variance, based on correlations) exhibited positive scores according to the first principal component, whereas TNAB (50.8%) and yeasts (29.8%) showed negative score values according to the first principal component (Figure 3a). A positive contribution to the second principal component calculation was observed for: yeasts (40.0% of the total variance, based on correlations) and molds (60.0%).

The PCA of the antioxidant data explained that the first two components accounted for 94.44% of the total variance (85.62% and 8.82%, respectively) in the six-variable factor space (antioxidant data). Considering the mapping of the PCA performed on the data, DPPH (which contributed 12.6% of the total variance, based on correlations), FRAP (12.3%), ABTS (13.2%), and TFC (14.2%) exhibited negative scores according to the first principal component (Figure 3b). A positive contribution to the second principal component calculation was observed for TPC (90.4% of the total variance, based on correlations).

According to Figure 4a, there is a positive correlation between total aerobic microbial count and yeast count (r = 0.342). On the other hand, correlation between TNAB and mold content is negative. There are positive correlations between DPPH, FRAP, ABTS, TPC and TFC (Figure 4b). The highest positive correlations were found between ABTS and TFC (r = 0.928), DPPH and TFC (r = 0.920), DPPH and FRAP (r = 0.857), (Figure 4b).

3.7. ANN Model

The calculated optimal neural network model for microbiological parameters, such as the number of aerobic bacteria (TNAB), yeasts and molds count showed adequate generalization capabilities for the modeling of experimental results: The optimum number of neurons in the hidden layer of ANN model was 10 (network MLP 7-10-3) (

Table 5. Artificial neural network model summary (performance and errors), for training, testing and validation cycles.

Network	Performance			Error			Training Algorithm	Error Function	Activation
	Train.	Test.	Valid.	Train.	Test	Valid.			Hidden	Output
MLP 7-10-3	0.742	0.659	0.792	5.4·107	2.1·107	8.8·107	BFGS 106	SOS	Tanh	Logistic
MLP 7-10-5	0.982	0.956	0.960	38.339	96.566	79.667	BFGS 226	SOS	Tanh	Tanh

), while the r² values were: 0.742; 0.659; and 0.792, accordingly, during the training, testing and validation cycles for output variables, for the training, testing and validation cycles for output variables.

The optimal neural network model for antioxidant parameters, such as DPPH, FRAP, ABTS, TPC and TFC, showed quite good generalization capabilities for the modeling of experimental results: the optimum number of neurons in the hidden layer of ANN model was 10 (network MLP 7-10-5) (Table 5), while the r² values were: 0.982; 0.956 and 0.960, accordingly, during the training, testing and validation cycles for output variables, for the training, testing and validation cycles for output variables.

The obtained r² values during the testing cycle were: 0.652; 0.799 and 0.780 for TNAB, yeasts and molds count modeling, while the obtained r² values for DPPH, FRAP, ABTS, TPC and TFC were: 0.978; 0.967; 0.954; 0.984 and 0.995.

The goodness of fit between experimental results and model-calculated outputs, represented as ANN performance (sum of r² between measured and calculated TNAB, Yeasts, Molds, DPPH, FRAP, ABTS, TPC and TFC), observed during training, testing and validation steps, are shown in

Table 6. The “goodness of fit” tests for the developed ANN model.

	χ2	RMSE	MBE	MPE	r2	Skew	Kurt	Mean	StDev	Var
TNAB	1.8 × 107	4.2 × 103	−4.0 × 102	0.266	0.685	−1.1 × 10	1.1 × 102	−4.0 × 102	4.2 × 103	1.8 × 107
Yeasts	1.8 × 106	1.3 × 103	1.3 × 102	0.505	0.702	1.1 × 10	1.1 × 102	1.3 × 102	1.3 × 103	1.8 × 106
Molds	7.273	2.685	−0.084	0.084	0.805	−0.492	0.172	−0.084	2.695	7.266
DPPH	7.091	2.648	−0.035	117.584	0.976	−0.378	0.155	−0.035	2.663	7.090
FRAP	7.091	2.648	−0.035	290.328	0.962	−0.379	0.155	−0.035	2.663	7.090
ABTS	7.091	2.649	−0.033	65.939	0.942	−0.380	0.155	−0.033	2.663	7.090
TPC	7.130	2.656	−0.055	0.337	0.979	−0.356	0.112	−0.055	2.670	7.127
TFC	7.091	2.648	−0.035	173.988	0.995	−0.378	0.155	−0.035	2.663	7.090

r²—coefficient of determination, χ²- reduced chi-square, MBE—mean bias error, RMSE—root mean square error and MPE—mean percentage error.

.

The ANN model predicted experimental variables (TNAB, yeasts and molds, DPPH, FRAP, ABTS, TPC and TFC) reasonably well for a broad range of the process variables (as seen in Figure 5, where the experimentally measured and ANN model predicted values of TNAB, Yeasts and Molds are presented).

The efficiency of the ANN model in modeling TNAB, yeasts and molds is graphically illustrated by scatter plots (Figure 5). In most scatter plots, data are distributed with large dispersion, indicating low prediction accuracy.

The results obtained from the database were fitted to the developed ANN model. Reduced chi-square (χ²), root mean square error (RMSE), mean bias error (MBE), mean percentage error (MPE), and coefficient of determination (r²) were calculated statistical parameters applied for the determination of fitting quality between database and the developed model. The particularly high values of r² and low values of χ², RMSE, MBE and MPE suggested adequate fit (Table 6). The ANN model showed better fit to DPPH, FRAP, ABTS, TPC and TFC data, according to relatively low χ², RMSE, MBE, and MPE, as well as the high r² values (Table 6).

The ANN models satisfactorily modelled experimental variables for various process variables.

For the ANN model, the model calculated TNAB, Yeasts and Molds, were not too close to the experimental values in most cases in terms of r² values, while the sum of squares (SOS) values acquired using the ANN model were of the same order of magnitude as experimental errors for the outputs mentioned in the literature [58,59,60].

The ANN model predicted experimental variables (DPPH, FRAP, ABTS, TPC and TFC) reasonably well for a broad range of the process variables (as seen in Figure 1, where the experimentally measured and ANN model predicted values of DPPH, FRAP, ABTS, TPC and TFC are presented).

The efficiency of the ANN model in modeling DPPH, FRAP, ABTS, TPC and TFC is graphically illustrated by scatter plots (Figure 5). In most scatter plots, data are distributed with large dispersion, indicating low prediction accuracy.

The developed ANN model for TNAB, yeast and mold modeling consisted of 113 weights-bias coefficients, while the developed ANN model for DPPH, FRAP, ABTS, TPC and TFC modeling consisted of 168 weights-bias coefficients showing the high nonlinearity of the system [61,62,63].

Table 7. Elements of matrix W₁ and vector B₁ (presented in the bias column).

	1	2	3	4	5	6	7	8	9	10
Time	−18.323	0.258	−14.374	−22.471	−14.094	−10.103	8.658	3.175	−3.196	−1.702
Treatment(C)	3.538	8.000	−6.033	−12.017	−4.536	−2.369	−1.394	4.059	6.748	−4.653
Treatment(F)	1.833	5.212	1.067	−0.556	−4.657	2.060	−5.744	−2.246	−2.978	−0.275
Treatment(S)	3.561	−7.635	11.311	13.295	10.044	−4.592	−2.174	−0.839	0.714	7.602
Spread(B)	−0.760	1.646	0.453	1.468	−14.473	4.619	5.329	0.063	−3.376	−0.170
Spread(C)	12.301	1.096	5.841	11.777	5.212	−5.080	−7.528	5.368	2.277	−1.214
Spread(U)	−2.669	2.745	0.004	−12.523	10.037	−4.370	−7.103	−4.475	5.645	4.143
Bias	8.892	5.623	6.253	0.725	0.857	−4.826	−9.264	0.861	4.466	2.752

presents the elements of matrix W₁ and vector B₁, while

Table 8. Elements of matrix W₂ and vector B₂ (presented in the bias column).

	1	2	3	4	5	6	7	8	9	10	Bias
TNAB	−9.068	3.236	5.956	10.256	3.307	−3.085	5.517	−7.085	−6.790	2.167	5.501
Yeasts	−7.647	−14.211	8.626	−10.375	−3.223	−1.412	−3.777	0.832	−0.424	1.886	−6.377
Molds	0.799	−0.104	−1.480	−3.937	−2.262	−7.018	8.270	2.259	1.484	1.906	−8.626

presents the elements of matrix W₂ and vector B₂. These were derived during the ANN model development using Equation (1). The goodness of fits between experimental and model-calculated results were shown in Table 3.

Table 9. Elements of matrix W₁ and vector B₁ (presented in the bias column).

	1	2	3	4	5	6	7	8	9	10
Time	7.100	0.257	−17.721	−0.075	8.369	4.650	14.590	−6.199	−0.817	19.669
Treatment(C)	−1.072	1.133	−1.250	0.815	−1.449	0.740	0.173	0.885	1.067	0.115
Treatment(F)	−2.672	−0.859	4.899	0.275	−0.666	−3.011	−7.300	−0.007	−0.370	−2.419
Treatment(S)	2.284	0.237	0.461	−0.230	0.524	2.230	5.004	0.983	−1.201	−1.500
Spread(B)	−0.307	0.350	0.425	−3.578	−0.255	0.318	0.697	4.467	−2.723	−1.019
Spread(C)	−1.359	−0.638	2.212	5.292	−1.463	−0.651	−4.015	−2.191	0.386	−1.378
Spread(U)	0.288	0.525	1.566	−0.860	0.144	0.249	1.175	−0.303	1.641	−1.371
Bias	−1.343	0.394	4.114	0.897	−1.678	−0.084	−2.205	1.913	−0.549	−3.760

presents the elements of matrix W₁ and vector B₁ (presented in the bias column), and

Table 10. Elements of matrix W₂ and vector B₂ (presented in the bias column).

	1	2	3	4	5	6	7	8	9	10	Bias
DPPH	2.078	−0.827	−0.268	−0.452	−0.234	0.866	−2.572	−0.200	0.377	−0.334	0.782
FRAP	1.303	−0.621	−0.297	−0.479	−0.041	0.526	−1.655	−0.289	0.139	−0.492	0.621
ABTS	1.432	−0.245	−0.241	−0.216	−0.144	0.432	−1.760	−0.131	0.054	−0.310	0.439
TPC	−0.470	0.307	−0.525	−0.656	0.502	0.692	−0.553	−0.966	−0.345	−1.278	0.580
TFC	2.325	−0.439	−0.303	−0.280	−0.321	0.393	−2.533	−0.114	0.113	−0.341	0.588

presents the elements of matrix W₂ and vector B₂ (bias) for the hidden layer, used for calculation in Equation (1).

The quality of the model fit was tested, and the residual analysis of the developed model was presented in Table 8. The ANN model had an insignificant lack of fit tests, which means the model satisfactorily predicted the pig carcass compositions. A high r² is indicative that the variation was accounted for and that the data fitted the proposed model.

Global Sensitivity Analysis—Yoon’s Interpretation Method

The effects of analytical method parameters (time, temperature, and applied coating) on the determination of output variables (microbiological profile, antioxidant activity, phenolic, and flavonoid content) was analyzed by employing Yoon’s global sensitivity equation corresponding to the weight coefficients of the obtained ANN model [64,65]. Following the global sensitivity analysis of a displayed ANN model, the graphical illustration of Yoon’s interpretation method results was shown in Figure 6. Time was the most positively influential parameter influencing yeasts and molds count, with an approximately relative importance of +39.83% and +46.52%, respectively. On the other hand, the time influence on the TNAB count was quite the opposite −16.71%. The most negative effect on yeasts and olds count was observed for spread (c) (−23.52% and −12.21%, accordingly), as shown in Figure 6a–c.

Furthermore, time was the most negatively influential parameter for antioxidant parameters (DPPH, FRAP, and ABTS), total phenols content and total flavonoids content, with approximate relative importance of −42.24%; −45.89%; −43.56%; −69.45% and −42.29%, respectively. On the other hand, sample stored at refrigerator temperature generated the enhanced antioxidant parameters (DPPH, FRAP and ABTS), and total phenols content and total compared to other treatments, expressing the positive influence of storage at refrigerated temperatures with the following relative influences: +20.79%; +20.90%; +20.70%; +15.96% and +20.85%, respectively, Figure 6d–h.

According to the global sensitivity analysis, it can be concluded that the most influential analytical method parameter was time.

4. Conclusions

The results of this paper prove that pumpkin oil cake can be used for the synthesis of biopolymer coating, but also as a carrier for the active component (+essential oil) in order to obtain active packaging. The positive impact of the pumpkin oil cake-based applied coating based on the sustainability of the treated grapes was confirmed, but the effectiveness of the applied active packaging was even greater, especially in terms of microbiological stability.

The exposed results are inspiring because as the approaches using mathematical modeling microbiological and antioxidant parameters during storage were found to be effective tool. The influence of time, treatment and coatings on microbiological and antioxidant parameters was assessed in the sensitivity analysis. The outcomes of this study reveal that number of aerobic bacteria, number of yeasts and molds, as well as antioxidant parameters (DPPH, FRAP and ABTS), total phenols content and total flavonoids content of biopolymers coating based on pumpkin oil cake, activated with Mentha piperita, can be modeled based on time, treatment and coatings, and during storage of packed grapes. This modeling should take place according to reasonably depleted χ², RMSE, MBE, and MPE, and the increased r² values. It was confirmed that the artificial neural network is appropriate for the modeling of output variables.