In Vitro and In Silico Evaluation of the Antimicrobial and Antioxidant Potential of Thymus pulegioides Essential Oil

The study was designed to analyze and evaluate the antioxidant and antibacterial properties of the essential oils of Thymus pulegioides L. grown in Western Romania. Thymus pulegioides L. essential oil (TPEO) was extracted by steam distillation (0.71% v/w) using a Craveiro-type apparatus. GC-MS investigation of the TPEO identified 39 different compounds, representing 98.46% of total oil. Findings revealed that thymol (22.89%) is the main compound of TPEO, followed by para-cymene (14.57%), thymol methyl ether (11.19%), isothymol methyl ether (10.45%), and beta-bisabolene (9.53%). The oil exhibits good antibacterial effects; C. parapsilosis, C. albicans, S. pyogenes, and S. aureus were the most sensitive strains. The antioxidant activity of TPEO was evaluated by peroxide and thiobarbituric acid value, 1,1-diphenyl-2-picrylhydrazyl radical (DPPH), [2,2′-azinobis(3-ethylbenzothiazoline-6-sulfonic acid) diammonium] (ABTS) radical scavenging assay, and beta-carotene/linoleic acid bleaching testing. The antioxidative data recorded reveal, for the first time, that TPEO inhibits primary and secondary oxidation products, in some particular conditions, better than butylated hydroxyanisole (BHA) with significant statistical difference (p < 0.05). Moreover, TPEO antioxidant capabilities in DPPH and ABTS assays outperformed alpha-tocopherol (p < 0.001) and delta-tocopherol (p < 0.001). Molecular docking analysis revealed that one potential target correlated with the TPEO antimicrobial activity was d-alanine-d-alanine ligase (DDl). The best scoring ligand, linalyl anthranilate, shared highly similar binding patterns with the DDl native inhibitor. Furthermore, molecular docking analysis also showed that the main constituents of TPEO are good candidates for xanthine oxidase and lipoxygenase inhibition, making the essential oil a valuable source for protein-targeted antioxidant compounds. Consequently, TPEO may represent a new potential source of antioxidant and antibacterial agents with applicability in the food and pharmaceutic industries.


Introduction
Food additives are substances deliberately added to food products for technological purposes (for example, preservation) that can be incorporated directly or indirectly into them, affecting their characteristics [1]. Additives such as butylated hydroxyanisole (BHA), butylated hydroxytoluene (BHT), gallic acid esters (e.g., propyl gallate, octyl gallate, and dodecyl gallate), sorbates, sulfites, and benzoates are widely employed to prolong the shelf-life of foodstuffs and to improve their overall quality and safety [2,3]. However, despite their benefits, food consumers need to be aware that some synthetic additives can be responsible for several life-threatening conditions (e.g., allergic reactions, asthma, nausea, or in some circumstances, even carcinogenesis) [4][5][6]. Therefore, the food industry has now focused on discovering natural alternatives to synthetic additives to manage food safety preservation without threatening consumers' health.
The plant kingdom represents a rich source of valuable natural bioactive compounds with multiple applications in the food and pharmaceutical industries. For instance, natural extracts such as essential oils (EOs) could be an appropriate alternative to synthetic antimicrobial medicines, especially antibiotics [7][8][9]. Furthermore, the ability of EOs to inhibit the growth of undesirable microorganisms and lipid oxidation in foodstuffs could solve the substitution problem for synthetic additives and meet the increasing demand from food customers for developing "clean label" products [10].
Consequently, the main aim of this research was to analyze: (i) the chemical profile of TPEO by GC-MS technique; (ii) the antimicrobial and antioxidant properties of the TPEO; and (iii) a potential protein-targeted mechanism of action correlated with the biological antimicrobial and antioxidant activity of TPEO, to find new sources of green preservative and/or antioxidant with applicability in the food and pharmaceutic industries.

Extraction of TPEO
The collected raw material was air-dried in the dark at room temperature (middle of 2020/beginning of 2021), then cut into approximately 2 cm long pieces, and stored in a paper bag until use. The dried plant material (100 g) was submitted to steam distillation for 4 h at 100 • C temperature using a Craveiro-type apparatus [32]. A glass boiler (3000 mL) filled with water (refilled as needed) and heated with electrical resistance produced the steam, introduced at the bottom of the glass extraction vessel (1000 mL). After steam passed through the plant material (positioned on a perforated plate set some centimeters away from the extraction vessel base), the steam and vaporized oil were condensed in a water-cooling system. Finally, TPEO and hydrosol (aqueous phase) were collected in a glass receiver (250 mL) equipped with a hydrosol overflow outlet and a water-cooling jacket to decrease the formation of artifacts due to overheating [33,34]. After separation by decantation, the TPEO was dried over anhydrous Na 2 SO 4 and stored at −18 • C in sealed vials for future analysis (yielding 0.71% v/w). TPEO (%) yield was calculated and expressed in mL/100 g of dried plant material.

Gas Chromatography-Mass Spectrometry Analysis
Chromatographic analysis was carried out on a gas chromatograph (HP6890 Gas-Chromatograph) coupled with mass spectrometry (HP5973 Mass Spectrometer). One µL from the diluted TPEO (1/1000 v/v, in n-hexane) was injected in splitless mode. The sample was carried out at a helium flow rate of 1 mL/min through a Bruker BR-5MS fused silica capillary (Bruker, Billerica, MA, USA). The column was heated in a GC oven at a temperature range of 50 • C to 300 • C with 6 • C/min, the final hold was 5 min, and the solvent delay was 3 min. MS source was set at 230 • C and 150 • C for MS Quad, and ionization energy was 70 eV. The scanned compounds were weighted between 50 to 550 amu. All peaks from the obtained chromatograms represented chemical compounds found in the analyzed sample. The chemical constituents of TPEO were identified by comparison of their linear retention indices with those reported in this literature [35] and by computer matching NIST2.0 (USA National Institute of Science and Technology) library software.  [36], these strains generated the majority of foodborne outbreaks across the EU member states.

Antibacterial Activity Assay
In the first stage, the antimicrobial effect of TPEO was determined by the Kirby-Bauer method based on the indications of the Clinical and Laboratory Standards Institute [37] and our previous studies [38]. From each standardized microbial suspension (0.5 MacFarland or approximately 10 8 bacteria/mL and 2 × 10 6 yeasts/mL), 100 µL were inoculated on Mueller-Hinton (MH) agar, or MHF (MH supplemented with defibrinated horse blood and β-nicotinamide adenine dinucleotide-bioMérieux, Marcy-l'Étoile, France). On the Antioxidants 2022, 11, 2472 4 of 21 inoculated culture media, sterile blank disks (∅ 7 mm; BioMaxima, Lublin, Poland) were placed, impregnated with 10 µLTPEO; then, these plates were incubated for 24 h at 37 • C (for bacteria) and 28 • C (for Candida strains). The antimicrobial activity was reflected by the diameters of the inhibitory zones formed around the disks impregnated with TPEO. Gentamycin and fluconazole disks served as positive controls, while DMSO (dimethylsulphoxide) was used as the negative control. Experiments were performed in triplicate for each tested strain.

Minimum Inhibitory Concentration (MIC)
Minimum inhibitory concentration (MIC) was determined by the broth microdilution method used according to the Clinical and Laboratory Standards Institute (CLSI) guidelines [39] for bacterial strains and the European Committee on Antimicrobial Susceptibility Testing (EUCAST) [40] for the Candida sp., respectively. The standardized inoculum was diluted, resulting in approximately 1-5 × 10 5 colony-forming units/mL (CFU). Serial dilutions of TPEO in DMSO were obtained with concentrations of 40, 20, 10, 5, 2.5, 1.25, 0.625 mg/mL. Then, 100 µL of each TPEO dilution was added to 400 µL MH or MHF broth and 500 µL microbial suspension, resulting in a final microbial inoculum of 0.5 × 10 5 CFU/mL. Following a 24 h incubation at 37 • C or 28 • C, the lowest concentration without visible growth was considered as MIC value. Tubes containing 100 µL DMSO, 500 µL microbial suspension, and 400 µL broth were used to control microbial growth. Experiments were performed in triplicate for each tested strain.

Minimum Bactericidal Concentration (MBC) and Minimum Fungicidal Concentration (MFC)
The method described by Danciu et al. and Jianu et al. [38,41], with minor modifications, was employed to determine the MBC and MFC. Up to 1 µL from each test tube with no visible growth was inoculated on Columbia agar supplemented with 5% sheep blood and incubated for 24 h at 37 • C for bacterial strains or on Sabouraud Dextrose agar supplemented with chloramphenicol and incubated at 28 • C for Candida strains. The lowest TPEO concentration at which 99.5% of the inoculated microorganisms were killed was considered MBC and MFC, respectively [42,43]. Experiments were performed in triplicate for each tested strain.

Sample Preparation
A total of 200 ppm (v/v) TPEO was added to 10 mL cold-pressed sunflower oil. Separately, 200 ppm (w/v) of two synthetic antioxidants (BHT and BHA), used as the positive controls, were added to 10 mL cold-pressed sunflower oil, according to Directive 2006/52/EC [44]. As a negative control, 10 mL of cold-pressed sunflower oil without any additive was utilized.

Peroxide Value (PV)
The PV (meq of oxygen kg −1 ) was carried out by the potentiometric method every four days up to 24 days using the ISO 27107:2010 [45]. Experiments were performed in triplicate.

Thiobarbituric Acid Value (TBA)
The TBA value (g malondialdehyde g −1 ) was carried out every four days up to 24 days, based on the previously described methods by Jianu et al. [46]. Briefly, each sample (2 g) was mixed with benzene (5 mL) and 0.67% TBA aqueous solution (4 mL). The mixtures were constantly shaken for 1 h using a mechanical shaker. Subsequently, the supernatant was incubated for 45 min in a boiling water bath. After cooling to ambient temperature, its absorbance was measured at 540 nm (Specord 210 Analytik Jena spectrophotometer). Experiments were performed in triplicate. With minor modifications, the scavenging effect on DPPH radical was carried out by the method reported earlier by Brand-Williams et al. [47]. Briefly, serial dilutions of the TPEO (ranging from 1.5 mg/mL to 2.93 µg/mL) were prepared in pure methanol, then 100 µL of each of them were mixed with 10 µL from a DPPH stock solution of 1 mg/mL. The samples were added to the wells of a 96-well plate and stored at 25 • C for 30 min in the dark. The absorbance was registered at the wavelength of 515 nm (Tecan i-control, 1.10.4.0 infinite 200Pro spectrophotometer). The following formula computed the percentage of DPPH inhibition: % inhibition = (A blank − A sample ) 100/A blank (A blank and A sample are the absorbances of the control and the absorbance of the test sample, respectively). Experiments were performed in triplicate with BHA, alpha-tocopherol and delta-tocopherol serving as positive controls and methanol as the negative control. BioDataFit 1.02 program (Chang Broscience Inc, Castro Valley, CA, USA) was employed to determine the IC 50 (µg/L).

Beta-Carotene/Linoleic Acid Bleaching Assay
The assay was performed as reported by Rădulescu et al. [48]. Briefly, a mixture of 0.5 mg β-carotene in 1 mL of chloroform, 200 mg Tween 40 and 25 µL linoleic acid was prepared. The chloroform was then evaporated under vacuum (at 40 • C for 5 min). Subsequently, the residue was diluted with 100 mL of 3% hydrogen peroxide aqueous solution to form a transparent yellowish emulsion. Finally, TPEO and BHA (serving as positive control) were individually dissolved in ethanol (2 mg ml −1 ), then 350 µL transferred in test tubes with 2.5 mL of the beta-carotene and linoleic acid mixture, stirred exhaustively, and incubated for 48 h at room temperature. The absorbance values were registered at 490 nm. Experiments were performed in triplicate. The Relative Antioxidant Activity (RAA %) was calculated as follows: RAA = A TPEO /A BHA , where A TPEO and A BHA are the absorption of TPEO and the absorbance of the BHA, respectively.

In Silico Prediction of Bioactivity and Molecular Docking Studies
Molecular docking simulations were carried out following in detail, the previously described method [42]. In brief, protein structures from the RCSB Protein Data Bank [50] ( Table 1) were optimized as docking targets. Discovery Studio 4.1 was used to create 3D structures of the 39 TPEO components (Dassault Systems BIOVIA, San Diego, CA, USA). Molecular docking was carried out with the PyRx v0.8 virtual screening software (The Scripps Research Institute, La Jolla, CA, USA) and Vina's encoded scoring function [51]. The target protein was docked with ten conformers of each input molecular structure. The calculated root means square deviation (RMSD) between the predicted and experimental native ligand docking poses for each case could not exceed a 2 Å threshold in order to validate our protocol. The docking grid box coordinates and size were chosen to best fit the active binding site (Table 1). Docking scores were generated as ∆G binding energy values (kcal/mol). Accelrys Discovery Studio 4.1 was also used to analyze protein-ligand binding interactions (Dassault Systems BIOVIA, San Diego, CA, USA).

Statistical Analysis
The antioxidant properties of TPEO were tested using ANOVA procedure followed by post-hoc analysis. Two-way ANOVA with main and interaction effects was performed in the case of MDA and peroxide values with sample and incubation period as main effects. Because the number of observations of sample per incubation period was low (only nine observations), the Shapiro-Wilk test for normality was applied on the ANOVA residuals. For the homogeneity of variances across groups, the Levene test was used. All the assumption checks regarding normality and homogeneity of variances, were successfully completed after the removal of two extreme outliers. Because the interaction effect was highly significant, the post-hoc analysis was performed at each incubation period in turn using Tukey's parametric test. One-way ANOVA was applied on DPPH and ABTS with antioxidant type as grouping variable in both situations. The assumptions regarding the normality of data and homogeneity of variances were checked with the same procedure as described above. In post-hoc analysis, Tukey's test was applied in the case of DPPH. In the case of ABTS, due to non-homogeneity of variances (Levene's test, p < 0.001) the Games-Howell test was applied. A similar approach was used to test the antimicrobial properties of TPEO: one-way ANOVA applied on the diffusion disk diameter with bacteria as main effect followed by post-hoc analysis using the Tukey test. The significance level considered for all the above-mentioned testing procedures was 0.05. The statistical software used was Jamovi (Version 2.2.5).

TPEO Chemical Composition
The steam distillation of T. pulegioides' dried aerial parts provided a brown yellowcolored oil with a 0.71% (v/w) yield and an agreeable odor reminiscent of T. vulgaris oil. Consistent results were obtained by Pavel et al. [14], that reported a yield of 0.7% (v/w) for EO obtained from T. pulegioides growing wild in Romania. Beicu et al. [52] reported that yields range from 0.44 to 0.49% (v/w) for Romanian T. pulegioides. Similar oil content has been reported for T. pulegioides growing wild in in Italy, 0.50 to 0.87% (v/w) [53] and Lithuania, 0.23 to 0.96% (v/w) [54]. Furthermore, the results compromise the European Pharmacopoeia requirements for Serpylli herba that mentioned a yield of at least 0.3% [55].

TPEO Antimicrobial Activity
The in vitro antimicrobial activity of the TPEO was qualitatively and quantitatively evaluated against eight strains of pathogenic bacteria and Candida spp. by the presence (or absence) of inhibition zones, diameters of inhibition zones, and MIC, MBC, and MFC values with the standard antimicrobial drugs gentamycin and fluconazole. DMSO was used as negative control, and the recorded antimicrobial effect of this compound was absent, therefore was not reported as separate set of results. The halo diameters of the inhibition zone caused by TPEO against the tested microorganism strains ranged between 15.3 ± 0.58 mm and 45.7 ± 2.52 mm ( Table 3), suggesting that the oil exerts low to moderate antimicrobial effects. The results showed that C. parapsilosis, C. albicans, S. pyogenes, and S. aureus were more sensitive than the other microorganism tested. Our results agree with previous investigations [7,13,14,52], which recorded that TPEO displayed antimicrobial activity against C. parapsilosis, C. albicans, S. pyogenes, and S. aureus. The tested strains' recorded MICs, MBCs, and MFCs were 1.25, 5, 10, and 20 mg/mL, respectively ( Table 3). The TPEO exhibited a moderate MIC for C. parapsilosis, C. albicans, S. pyogenes, and S. aureus and showed low activity against the rest of the analyzed bacteria. Furthermore, our data is consistent with previous reports in which EOs exhibited comparable or higher antimicrobial activities, even when tested on drug-resistant bacterial strains such as methicillin-resistant S. aureus or carbapenem-resistant P. aeruginosa [61][62][63]. Overall, the Gram-positive strains appeared to be more susceptible than the Gram-negative. These differences might be generated by the presence of lipopolysaccharide molecules in the outer membrane of the Gram-negative bacteria cell wall that forms a hydrophilic permeability barrier providing protection against the effects of macromolecules and highly hydrophobic compounds [64][65][66]. The ability of these compounds to disrupt the permeability barrier of the bacterial cell walls, accompanied by the loss of chemiosmotic control, may represent the most likely explanation for the volatile oil's lethal action [67].

TPEO Antioxidant Activity
Hydroperoxides are the primary products of lipid oxidation. Therefore, peroxide value can be employed as a direct marker for the primary phases of lipid oxidation. A higher peroxide value indicates lower oxidative stability of fats and oils during storage [68]. The antioxidant activity of TPEO has been carried out on cold-pressed sunflower oil. Figure 1A displays the progressive increase in PVs throughout the storage period of the studied sunflower oil samples. The PV for the sample treated with TPEO, at zero days incubation period, is not significantly different from samples treated with BHT (p = 0.979) and BHA (p = 0.402) according to the Tukey test. At four days, PV value for TPEO is not significantly different from the control sample (p = 0.407) and from BHA (p = 0.957). However, more promising results are found at 20 days of incubation, where PVs for the sample treated with TPEO are lower and significantly different than BHA values (p < 0.001) and lower but not significantly different than BHT values (p = 0.64).  The propagation phase generates considerable secondary reaction products, particularly malondialdehyde (MDA). Therefore, MDA concentration is a direct marker for the second phase of lipid peroxidation [69]. In addition, it has been suggested that MDA in food and edible oils has harmful effects on health [70]. Figure 1B   The propagation phase generates considerable secondary reaction products, particularly malondialdehyde (MDA). Therefore, MDA concentration is a direct marker for the second phase of lipid peroxidation [69]. In addition, it has been suggested that MDA in food and edible oils has harmful effects on health [70]. Figure 1B displays the progressive increase in MDA throughout the storage period of the studied sunflower oil samples. The statistical analysis shows significant differences in MDA values between the sample treated with TPEO on one part and samples treated with BHA and BHT on the other part for all storage periods from 4 days through 24 days. Unfortunately, the MDA values for the sample treated with TPEO are higher than those treated with BHT and BHA. However, better results are recorded for the TPEO sample at zero days of incubation. More precisely, MDA values for the sample treated with TPEO are lower but not significantly different than MDA values for BHA (p = 0.252).
DPPH radical scavenging assay is a broadly used method for determining the antioxidant activity of natural products or synthetic compounds [71]. The procedure is based on reducing DPPH in solution with a hydrogen-donating antioxidant and is dependent on the formation of the nonradical form DPPH-H in the reaction [72]. TPEO reduced the stable free radical DPPH with an IC 50 value of 19.13 ± 0.61 µg/mL (  The ABTS radical scavenging method is broadly used to analyze the antioxidant activity of single compounds and complex mixtures of various plants [82][83][84]. In the ABTS assay, TPEO indicated a strong antioxidant activity with an IC 50 of 1.66 ± 0.1 µg/mL (Table 4), which was significantly (p < 0.001) more pronounced than that of alpha-tocopherol (IC 50 2.08 ± 0.1 µg/mL) and delta-tocopherol (p < 0.001) (IC 50 1.99 ± 0.1 µg/mL) according to the Games-Howell test. Our results are comparable with the findings of Taghouti et al. for T. pulegioides hydroethanolic extracts [85] and better than the findings for EO of T. algeriensis, T. polium, T. vulgaris, T. broussonnettii, T. willdenowii, T. citriodorus [80,[86][87][88].
In the beta-carotene-linoleic acid bleaching assay, beta-carotene undergoes rapid discoloration without an antioxidant. Antioxidants can delay the extent of beta-carotene destruction by "neutralizing" the linoleate free radical or additional free radicals developed within the system [89]. Table 4 shows the inhibition of beta-carotene bleaching by the TPEO and the positive control (BHA). TPEO exhibited strong antioxidant activity (89.62 ± 0.14%) in the beta-carotene-linoleic acid test but lower than that of BHA (100%) ( Table 4), with no significant differences (p > 0.05) observed. No previous studies were recorded in the literature regarding TPEO activity in the beta-carotene/linoleic acid system to allow us to make direct comparisons.
Given the preceding results, it is hard to ascertain the precise impact of each constituent on the antioxidant and antimicrobial activity of the essential oil. However, several of the TPEO's major components, including thymol, thymol methyl ether, beta-bisabolene gamma-terpinene, p-cymene, or carvacrol, have previously been shown to exhibit potent antioxidant effects and antimicrobial properties on a variety of microbial strains [90][91][92][93][94], which can explain the overall synergistic effect of the TPEO observed in the present study.

In Silico Prediction of Bioactivity and Molecular Docking Studies
Computational methods are now an important tool for reducing the time required to unravel the action mechanisms of pharmacologically active substances. Molecular docking is a computational technique that allows the user to dock candidate molecules into the active site of a biological target and then classify the compound set based on their binding affinity, which is calculated using a scoring function [95]. Given that the essential oil demonstrated an antibacterial effect against two tested bacterial strains (S. aureus, S. pyogenes) and two tested fungal strains (C. albicans, C. parapsilosis), molecular docking was employed to establish potential target hits for the 39 components of TPEO. In addition, given its broad antioxidant versatility, using this strategy, we also intended to deduce potential targets that can be related to a protein-targeted antioxidant effect. For our current purpose, we used two sets of target proteins. The first set of proteins pertains to druggable targets that are frequently used for antimicrobial drug design and includes Isoleucyl-tRNA synthetase (IARS), DNA gyrase, D-alanine: D-alanine ligase (DDl), Streptococcus pneumoniae (DHPS), type IV topoisomerase, dihydrofolate reductase (DHFR), DNA gyrase subunit B, penicillinbinding protein 1a (PBP1a). The second set of targets consists of proteins that play critical roles in the metabolic generation of reactive oxygen species as byproducts and whose inhibition can reduce metabolic oxidative stress. This set includes lipoxygenase, CYP2C9, NADPH-oxidase, and xanthine oxidase. Table 5 displays the docking scores obtained for the 39 docked compounds.  The goal of our in silico-based method was to identify protein targets that can be inhibited by multiple components of our TPEO, or at least by the essential oil's major constituents. Finding correlations within docking scores is difficult because each protein has different binding site characteristics and native ligands give different docking scores, resulting in each score set having different control values. To mitigate these inconveniences, we calculated each docking score as a percentage of its respective native ligand score (considered 100%). These percentages were represented as a radar graph, with the scores of each compound representing plot lines and the protein targets forming the corners of the radar chart. The final result should reveal the most likely targeted protein as a graph line stretch of the vast majority of all plotted series towards the protein targets corner. The graphical representation of docking scores for the antimicrobial protein targets are depicted in Figure 2A. When the scores of all 39 TPEO compounds are graphically represented, it is difficult to notice a tendency of a majority of the TPEO components to target a specific protein (Figure 2A). However, if we represent only the 7 major components (which account for more than 75% of the essential oil content) in the same graph, we can clearly see a tendency of the compounds to inhibit DDl (2I80) ( Figure 2B). In this case, when compared to the native ligand (3-chloro-2,2-dimethyl-n-[4-(trifluoromethyl)phenyl]propenamides, −7.2 kcal/mol), the most active compounds were 24 (Linalyl anthranilate, −7 kcal/mol), 36 (beta-Sesquiphellandrene), and the major component 33 (beta-Bisabolene, -6.8 kcal/mol). components to target a specific protein (Figure 2A). However, if we represent only the 7 major components (which account for more than 75% of the essential oil content) in the same graph, we can clearly see a tendency of the compounds to inhibit DDl (2I80) ( Figure  2B). In this case, when compared to the native ligand (3-chloro-2,2-dimethyl-n- [4-(trifluoromethyl)phenyl]propenamides, −7.2 kcal/mol), the most active compounds were 24 (Linalyl anthranilate, −7 kcal/mol), 36 (beta-Sesquiphellandrene), and the major component 33 (beta-Bisabolene, -6.8 kcal/mol).  Due to their crucial and universal involvement in bacterial cell-wall peptidoglycan production, which has been a proven target for antibiotics, DDl enzymes are a promising target for further chemotherapeutic research [96]. Based on known DDl structures, the conserved His, Val and Glu residues (His96, Val19, Glu16) have an essential role in the catalytic formation of D-alanyl-D-alanine. The native ligand of the DDl structure used for docking is an allosteric inhibitor that interacts with the aforementioned region within the binding site, rendering the proper catalytic activity of DDl ineffective. The only hydrogen bond (HB) formed by the native ligand is with Pro311 [97]. Similarly, compound 24 (linalyl anthranilate) interacts with the same binding site region and features the same HB with Pro311 making it a good candidate for DDl inhibition (Figure 3). Recent literature also revealed that cell membrane disruption and cell content leakage were two mechanisms related to linalyl anthranilate's antibacterial effect against K. pneumoniae [98]. Terpenerich essential oils exhibit antimicrobial activity and, due to their lipophilic structure, their bactericidal mechanism of action includes cell membrane disruption or impairment of cell membrane proteins [99]. Beta-bisabolene (compound 33) was also previously reported as an active compound against S. aureus, which can also restore the antibacterial effect of ampicillin against β-lactam resistant S. aureus strains [100]. While it is unclear that the mechanisms of action are related to bacterial cell membrane impairment, one of the TPEO mechanisms of action may correlate to the inhibition of DDl.
Using the same docking-based method, we also wanted to highlight a potential protein-targeted biological antioxidant effect of TPEO. By plotting all the 39 compounds' docking scores, the graph shows that the majority of compounds tend to target two proteins: xanthine oxidase (3NRZ) and lipoxygenase (1N8Q) ( Figure 4A). This targeted tendency towards the two proteins becomes much more significant when the chart only records docking scores of the 7 major compounds, previously mentioned ( Figure 4B). essential oils exhibit antimicrobial activity and, due to their lipophilic structure, their bactericidal mechanism of action includes cell membrane disruption or impairment of cell membrane proteins [99]. Beta-bisabolene (compound 33) was also previously reported as an active compound against S. aureus, which can also restore the antibacterial effect of ampicillin against β-lactam resistant S. aureus strains [100]. While it is unclear that the mechanisms of action are related to bacterial cell membrane impairment, one of the TPEO mechanisms of action may correlate to the inhibition of DDl. Using the same docking-based method, we also wanted to highlight a potential protein-targeted biological antioxidant effect of TPEO. By plotting all the 39 compounds' docking scores, the graph shows that the majority of compounds tend to target two proteins: xanthine oxidase (3NRZ) and lipoxygenase (1N8Q) ( Figure 4A). This targeted tendency towards the two proteins becomes much more significant when the chart only records docking scores of the 7 major compounds, previously mentioned ( Figure 4B).  The oxidation reactions of polyunsaturated fatty acids catalyzed by lipoxygenase enzymes are one of several processes that generate reactive oxygen species (ROS) as byproducts. Lipoxygenases play an important role in pro-inflammatory signaling by regulating ROS levels [101]. Therefore, lipoxygenase inhibition can be therapeutically beneficial in both reducing ROS levels and the correlated inflammatory response. From the TPEO major components, compounds 12 (para-cymene), 16 (gamma-terpinene), and 26 (carvacrol) exhibited better docking scores when compared to the lipoxygenase native ligand. The most active in silico recorded compound was carvacrol. Ligand-protein interaction analysis shows carvacrol interacting with Glu802 and Ala1079 through HB and also forms other hydrophobic interactions with Phe914, Phe1009, and Ala1078 ( Figure 5). The oxidation reactions of polyunsaturated fatty acids catalyzed by lipoxygenase enzymes are one of several processes that generate reactive oxygen species (ROS) as byproducts. Lipoxygenases play an important role in pro-inflammatory signaling by regulating ROS levels [101]. Therefore, lipoxygenase inhibition can be therapeutically beneficial in both reducing ROS levels and the correlated inflammatory response. From the TPEO major components, compounds 12 (para-cymene), 16 (gamma-terpinene), and 26 (carvacrol) exhibited better docking scores when compared to the lipoxygenase native ligand. The most active in silico recorded compound was carvacrol. Ligand-protein interaction analysis shows carvacrol interacting with Glu802 and Ala1079 through HB and also forms other hydrophobic interactions with Phe914, Phe1009, and Ala1078 ( Figure 5). These results are in line with previously reported studies that highlighted the anti-lipoxygenase activity of essential oils due to their carvacrol and para-cymene content or the anti-lipoxygenase activity of carvacrol alone [102][103][104][105]. Among the same major TPEO components, carvacrol also exhibited the highest in silico inhibitory activity against xanthine oxidase. Xanthine oxidase is a widely distributed enzyme that catalyzes the oxidation of hypoxanthine to xanthine and xanthine to uric acid. Mammalian xanthine oxidase is another physiological source of ROS, such as superoxide ion, hydrogen peroxide, and nitric oxide [106]. Ligand-protein interaction analysis revealed that carvacrol is tightly anchored within the protein binding pocket, mostly by hydrophobic interactions. Carvacrol also forms 3 HB with Gln514, His518, and Trp519 ( Figure 6). Our results are in line with previously reported studies that highlighted both the in silico and in vitro inhibitory potential of carvacrol against xanthine oxidase [107,108]. Regarding those exposed, TPEO can be a significant source of chemicals with an antioxidant potential characterized by the suppression of ROS-producing proteins such as lipoxygenase and xanthine oxidase. Among the same major TPEO components, carvacrol also exhibited the highest in silico inhibitory activity against xanthine oxidase. Xanthine oxidase is a widely distributed enzyme that catalyzes the oxidation of hypoxanthine to xanthine and xanthine to uric acid. Mammalian xanthine oxidase is another physiological source of ROS, such as superoxide ion, hydrogen peroxide, and nitric oxide [106]. Ligand-protein interaction analysis revealed that carvacrol is tightly anchored within the protein binding pocket, mostly by hydrophobic interactions. Carvacrol also forms 3 HB with Gln514, His518, and Trp519 ( Figure 6). Our results are in line with previously reported studies that highlighted both the in silico and in vitro inhibitory potential of carvacrol against xanthine oxidase [107,108]. Regarding those exposed, TPEO can be a significant source of chemicals with an antioxidant potential characterized by the suppression of ROS-producing proteins such as lipoxygenase and xanthine oxidase.

Conclusions
The research assessed the chemical composition, antibacterial and antioxidant activity of TPEO. The analyzed data shows that the volatile oil belongs to the thymol chemotype and exhibited good antimicrobial effects against C. parapsilosis, C. albicans, S. pyogenes, and S. aureus. The antioxidative data recorded reveal for the first time that TPEO, if compared to BHA in some specific conditions, inhibits primary and secondary oxidation products with the same efficacy, i.e., no statistically significant difference (p > 0.05), or even better with significant difference (p < 0.05). Moreover, TPEO antioxidant capabilities in DPPH and ABTS assays were better than those of alpha-tocopherol (p < 0.001) and deltatocopherol (p < 0.001). Molecular docking-based in silico analysis revealed that one of the antibacterial mechanisms of TPEO could be attributed to DDl inhibition by both minor and major constituents. Docking determinations also revealed that major components of TPEO may potentially be biologically active antioxidant compounds by inhibiting ROSproducing enzymes such as lipoxygenase and xanthin oxidase. In conclusion, the results recommend TPEO as a new source of natural antioxidant and antibacterial agents with applicability in the food and pharmaceutic industries.

Conclusions
The research assessed the chemical composition, antibacterial and antioxidant activity of TPEO. The analyzed data shows that the volatile oil belongs to the thymol chemotype and exhibited good antimicrobial effects against C. parapsilosis, C. albicans, S. pyogenes, and S. aureus. The antioxidative data recorded reveal for the first time that TPEO, if compared to BHA in some specific conditions, inhibits primary and secondary oxidation products with the same efficacy, i.e., no statistically significant difference (p > 0.05), or even better with significant difference (p < 0.05). Moreover, TPEO antioxidant capabilities in DPPH and ABTS assays were better than those of alpha-tocopherol (p < 0.001) and delta-tocopherol (p < 0.001). Molecular docking-based in silico analysis revealed that one of the antibacterial mechanisms of TPEO could be attributed to DDl inhibition by both minor and major constituents. Docking determinations also revealed that major components of TPEO may potentially be biologically active antioxidant compounds by inhibiting ROS-producing enzymes such as lipoxygenase and xanthin oxidase. In conclusion, the results recommend TPEO as a new source of natural antioxidant and antibacterial agents with applicability in the food and pharmaceutic industries.