Next Article in Journal
Synthesis, Characterization, and Anticancer Activities Evaluation of Compounds Derived from 3,4-Dihydropyrimidin-2(1H)-one
Next Article in Special Issue
Essential Oil Composition and Biological Activity of “Pompia”, a Sardinian Citrus Ecotype
Previous Article in Journal
Isolation and Purification of Bioactive Compounds from the Stem Bark of Jatropha podagrica
Previous Article in Special Issue
Volatiles Profiling, Allelopathic Activity, and Antioxidant Potentiality of Xanthium Strumarium Leaves Essential Oil from Egypt: Evidence from Chemometrics Analysis
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Machine Learning Analyses on Data including Essential Oil Chemical Composition and In Vitro Experimental Antibiofilm Activities against Staphylococcus Species

1
Rome Center for Molecular Design, Department of Drug Chemistry and Technology, Sapienza University, P.le Aldo Moro 5, 00185 Rome, Italy
2
Alchemical Dynamics s.r.l., 00125 Rome, Italy
3
Department of Public Health and Infectious Diseases, Sapienza University, P.le Aldo Moro 5, 00185 Rome, Italy
4
Faculty of Natural Sciences and Mathematics, University of Montenegro, 81000 Podgorica, Montenegro
5
Department of Drug Chemistry and Technology, Sapienza University, P.le Aldo Moro 5, 00185 Rome, Italy
6
Master Course in Cosmetic Sciences, Department of Life Sciences and Biotechnology, University of Ferrara, 44100 Ferrara, Italy
*
Authors to whom correspondence should be addressed.
These authors contributed equally to this work.
Molecules 2019, 24(5), 890; https://doi.org/10.3390/molecules24050890
Submission received: 10 February 2019 / Revised: 25 February 2019 / Accepted: 26 February 2019 / Published: 3 March 2019
(This article belongs to the Special Issue Biological Activities of Essential Oils)

Abstract

:
Biofilm resistance to antimicrobials is a complex phenomenon, driven not only by genetic mutation induced resistance, but also by means of increased microbial cell density that supports horizontal gene transfer across cells. The prevention of biofilm formation and the treatment of existing biofilms is currently a difficult challenge; therefore, the discovery of new multi-targeted or combinatorial therapies is growing. The development of anti-biofilm agents is considered of major interest and represents a key strategy as non-biocidal molecules are highly valuable to avoid the rapid appearance of escape mutants. Among bacteria, staphylococci are predominant causes of biofilm-associated infections. Staphylococci, especially Staphylococcus aureus (S. aureus) is an extraordinarily versatile pathogen that can survive in hostile environmental conditions, colonize mucous membranes and skin, and can cause severe, non-purulent, toxin-mediated diseases or invasive pyogenic infections in humans. Staphylococcus epidermidis (S. epidermidis) has also emerged as an important opportunistic pathogen in infections associated with medical devices (such as urinary and intravascular catheters, orthopaedic implants, etc.), causing approximately from 30% to 43% of joint prosthesis infections. The scientific community is continuously looking for new agents endowed of anti-biofilm capabilities to fight S. aureus and S epidermidis infections. Interestingly, several reports indicated in vitro efficacy of non-biocidal essential oils (EOs) as promising treatment to reduce bacterial biofilm production and prevent the inducing of drug resistance. In this report were analyzed 89 EOs with the objective of investigating their ability to modulate bacterial biofilm production of different S. aureus and S. epidermidis strains. Results showed the assayed EOs to modulated the biofilm production with unpredictable results for each strain. In particular, many EOs acted mainly as biofilm inhibitors in the case of S. epidermidis strains, while for S. aureus strains, EOs induced either no effect or stimulate biofilm production. In order to elucidate the obtained experimental results, machine learning (ML) algorithms were applied to the EOs’ chemical compositions and the determined associated anti-biofilm potencies. Statistically robust ML models were developed, and their analysis in term of feature importance and partial dependence plots led to indicating those chemical components mainly responsible for biofilm production, inhibition or stimulation for each studied strain, respectively.

Graphical Abstract

1. Introduction

A biofilm is a microbially derived sessile community characterized by cells irreversibly attached to a substrate or interface or to each other, embedded in a self-produced matrix of extracellular polymeric substances, which exhibits an altered phenotype with regard to growth, gene expression and protein production [1]. Biofilm resistance to antimicrobials [2] is a complex phenomenon, driven not only by genetic mutation induced resistance, but also by means of increased microbial cell density that supports resistance by means of horizontal gene transfer across cells [3]. Indeed, other mechanisms are involved, such as: (i) low penetration of antimicrobial agents due to the barrier function exerted by the biofilm matrix, (ii) presence of cells exhibiting a high multidrug tolerance, (iii) reduced susceptibility to antibiotics as a consequence of stress adaptive responses or changes in the chemical biofilm microenvironment [4]. The strategies adopted to treat these challenging infections are rapidly changing due to the increasing understanding of biofilm structure and functions. Nonetheless, the prevention of biofilm formation and the treatment of existing biofilms is currently a difficult challenge; therefore, the discovery of new multi-targeted or combinatorial therapies is increasingly urgent [5].
The development of anti-biofilm agents is therefore considered of major interest and represents an important strategy since non-biocidal molecules to avoid the rapid appearance of resistant mutants are highly valuable. Among bacteria, staphylococci are prevalent causes of biofilm-associated infections [6]. In particular, Staphylococcus aureus (S. aureus) is an opportunistic pathogen that can cause serious diseases in humans, ranging from skin and soft tissue infections to invasive infections of the bloodstream, heart, lungs and other organs [7]. In 2013, Nicholson et al. reported that 30% of U.S. population was colonized by S. aureus while 1.5% was found to be a carrier of methicillin-resistant S. aureus (MRSA), a major cause of healthcare-related infections responsible for a significant proportion of nosocomial infections worldwide. Recently in the U.S. deaths from MRSA infections have exceeded those from many other infectious diseases, including HIV/AIDS [8]. Staphylococcus epidermidis (S. epidermidis), conventionally considered a commensal of human skin, can cause significant problems when breaching the epithelial barrier, especially during biofilm-associated infection of indwelling medical devices [9,10]. Most diseases caused by S. epidermidis exhibit a chronic profile and occur as device-related infections (such as intravascular catheter or prosthetic joint infections) and/or their complications [10]. In view of the above scenario the scientific community is seeking for new agents endowed of anti-biofilm capabilities to fight S. aureus and S epidermidis infections. Recently, several reports indicated in vitro efficacy of non-biocidal essential oils (EOs) as promising treatment to reduce bacterial biofilm production and prevent the inducing of drug resistance [11]. In different applications, EOs have been found of some efficacy in reducing biofilm production of either S. aureus standard strains or MRSA [12,13,14,15,16,17]. In other reports, EOs and some of their purified chemical components have also been proved to inhibit S epidermidis biofilm production [18,19,20].
Recently machine learning (ML) has been proved as a tool able to deeply investigate the modulatory role of EOs’ chemical components on Pseudomonas aeruginosa biofilm production [21,22,23,24]. In particular, 89 EOs extracted in different periods and times of extractions from three different plants were analyzed: 13 EOs (RSEOs) from Ridolfia segetum Moris (RS); 32 EOs (FVEOs) from Foeniculum vulgare Miller (FV) and 44 EOs (CGEOs) from Calamintha nepeta (L.) Savi subsp. glandulosa (Req.) Ball, (CG) [22]. In line with that study and with the objective of investigating EOs’ ability to also reduce bacterial biofilm production in other bacteria, herein is reported an extensive study of the 89 EOs samples as potential antibacterial and anti-biofilm agents against S. aureus ATCC 6538P, S. aureus ATCC 25923, S. epidermidis RP62A and S. epidermidis O-47. To this purpose, like previously reported [22], ML algorithms were applied to the EOs’ chemical compositions and the determined associated anti-biofilm potencies, with the purpose of shedding light on those components likely mainly responsible for either positive or negative modulation of biofilm production.

2. Results

2.1. Antimicrobials Activity of EOs

Some FVEOs and CGEOs samples showed MICs at the highest used concentration (Table 1 and Table 2). No antimicrobial activity on staphylococci was recorded for any of the RSEOs, except for the R30 sample that showed a MIC value of 25 mg/mL (Table 3). Only seven out 89 EOs displayed MIC values of 6.25 mg/mL against the two S. aureus strains (Table 1, Table 2 and Table 3). As control, MIC was also evaluated for conventional a antibiotic, ofloxacin, belonging to the fluoroquinolone family.

2.2. Biofilm Production Modulation by EOs at Selected Fixed Concentrations

Preliminarily, the same representative EOs (2 RSEOs, 3 CGEOs and 3 FVEOs) among the reported 89 used on P. aeruginosa [22] were selected, to evaluate the anti-biofilm potency at different concentrations starting from 25 mg/mL, using scalar dilutions (data not shown).
The obtained preliminary data analyzed in terms of biofilm production modulation and reproducibility led to the selection of two representatives concentrations (3.125 mg/mL and 0.0488 mg/mL). The first concentration was in the range of milligrams while the second one was in the range of micrograms. All 89 EOs were then tested at the two selected concentrations and the biofilm production was measured relatively to untreated bacteria (Figure 1, Figure 2 and Figure 3).
At either selected concentrations EOs modulated the biofilm production with unpredictable results for each strain. These results anticipated that many EOs may act mainly as biofilm inhibitors in the case of RP62A and O-47 strains, while for 6538P and 25923 EOs can either induce no effect or stimulate biofilm production (Table 4). In Table 4, the number of EOs able to inhibit (<100%, <80% and <50%, respectively) or stimulate (≥100%, ≥120%, ≥150% and ≥200%, respectively) biofilm formation is reported. It is worthy to note that on S. epidermidis strains about 30 EOs inhibited more than 50% of biofilm growth even at lowest concentration, while almost none of them showed an activity on S. aureus strains.

2.3. Quantitative Analysis of Selected EOs against Different Strains of S. epidermidis

Representative EOs selected among those able to reduce more than 70% of biofilm formation were further analyzed to evaluate a dose-dependent effect against S. epidermidis RP62A and O-47 (Figure 4, Figure 5 and Figure 6). The inhibition by RSEOs was confirmed at lower concentrations on both strains despite their different biofilm matrix composition and the inhibition of biofilm formation was clearly not dose-dependent (Figure 4). Analogous results were obtained with FVEOs samples (Figure 5).
Differently, CGEOs revealed a dose dependent biofilm inhibition being more pronounced on the strongest biofilm producer S. epidermidis O-47 than on S. epidermidis RP62A (Figure 6).

2.4. Application of Machine Learning Algorithms

2.4.1. PCA Analysis of Datasets

Extraction of the first 3 PCs afforded to a cumulative explained variance of almost 90% (PC1:61.18%; PC1 + PC2: 75.02%; PC1 + PC2 + PC3: 82.92%). As reported [22], the first two principal components indicate at least three clusters (Supplementary Material Figure SM-1) and correctly identified the three plants derived EOs, although no definite separation between RSEOs and FVEOs was observable. PCA related loading plots indicated estragole, o-cymene and pulegone as the chemical components mainly related to high PCs values (Supplementary Material Figure SM-1).

2.4.2. Binary Classification Models

General Results

Analogously as in ML application to Pseudomonas aeruginosa (PA) [22] direct application of linear classification methods using algorithms such as Logistic Regression (LR) and Linear Support Vector Machines (SVM) [25] did not lead to satisfying classifiers (data not shown). At the same time non-linear algorithms like random forest (RF) [26], non-linear support vector machine (SVM) [27] and gradient boosting (GB) [28] also led to insufficiently robust models (data not shown). Therefore, a mixed approach was used and taking the idea from the principal component regression (PCR) as an evolution of multiple linear regression (MLR) a number of PCs were used in place of the original variables (EOs chemical component percentages) as input for the sklearn LR implementation (PCLR). As an initial test, the PCLR was run on the PA dataset leading to highly overlapping results with those obtained with the GB application (data not shown). Nevertheless, as biofilm production assay profiled EOs as either inhibitors of activators (Table 4) accordingly, classification models were tentatively built for all four strains considering either biofilm production inhibition or activation for biofilm percentages observed at the two above introduced concentration levels of 48.8 μg/mL and 3.125 mg/mL. To this aim, initially the optimal biofilm production percentage cutoff for the binary classification was explored by systematically either decreasing it from a starting 80% to 60% or increasing from 120% to 140% for the inhibition or activation models, respectively, being the ranges arbitrarily chosen on the basis of Table 4 filled data. The models’ accuracy was monitored by the MCC value obtained by leave-one-out cross-validation. Following this protocol, for the inhibition models EOs samples characterized by higher values than the best performing cutoff of biofilm production percentage were classified as inactive, while those with lower values were considered active. On the contrary, regarding the biofilm production enhancer models, EOs samples characterized by higher percentages than the cutoff value were classified as active, while those with lower values were considered non-active. Regarding the 6538P inhibition training set, the very low active/inactive ratio at biofilm inhibition below 80% prevented any optimization. Therefore, the grid search analysis to the starting sixteen training sets (four strains by two series of models by two concentrations) afforded to seven optimized models for either concentrations (Table 5). Inspection of optimized models on both hyperparameters and cutoff values revealed for 25923/inhibition, RP62A/activation and O-47/activation sets composed of high unbalanced ratios of actives over non-actives and were hence not further analyzed. Comparing developed models for the two used EOs concentrations revealed 3.125 mg/mL level to lead to more reliable and robust models (Table 6). Based on the above preliminary data, subsequent results and analyses were only carried out on RP62A/inhibition, O-47/inhibition, 6538P/activation and 259237/activation models derived for biofilm modulation recorded at 3.125 mg/mL. This is in full agreement with the fact that EOs samples acted prevalently as reducer of biofilm production for RP62A and O-47 strains, while for 6538P and 25923 the biofilm production was mainly enhanced (Table 4).
To assess either models’ fitness and robustness, their lack of chance correlation was checked by Y-scrambling procedure whose 100 runs of cross-validated scrambled set led to average, standard deviation, maximum and minimum values for AccuracyY-S, MCCY-S, Precision-RecallY-S and ROC-AUCY-S ROC-AUC coefficients always lower than non-cross-validated and cross-validated ones, therefore assessing validity of all final models.

Binary Classification Model for 6538P Biofilm Production Activation

The 6538P/activation/3.125mg/mL optimized derived model was maximum at a cutoff of 133%, using 9 PCs, characterized by a 27:62 (0.44) proportion between actives and non-actives and high values of Accuracy (0.832), MCC (0.667), Precision-Recall (0.772) and ROC-AUC (0.824) coefficients which persisted to be quite good in cross-validation (AccuracyCV = 0.805, MCCCV = 0.613, Precision-RecallCV = 0.698 and ROC-AUCCV = 0.743) (Table 6 and Table 7). Feature importance and partial dependence pointed out compounds 3-octanol, d-limonene and pulegone as more important for biofilm production enhancement modulation (Figure 7, Supplementary Material Figure SM-1 and Table 8).

Binary Classification Model for RP62A biofilm production inhibition

The grid search on the EOs’ chemical composition and their associated RP62A biofilm production inhibitory potencies at 3.125 mg/mL, identified 62% biofilm residual production as the best cutoff value with only 5 PCs and a actives over non-actives ratio of 31:58 (0.53). The final classification model was found characterized by Accuracy, MCC, Precision-Recall and ROC-AUC values of 0.721, 0.455, 0.657 and 0.742, respectively (Table 6). Cross-validation associated coefficients AccuracyCV, MCCCV, Precision-RecallCV and ROC-AUCCV were 0.687, 0.392, 0.584 and 0.683, respectively.
Inspection of model associated EOs’ chemical components importance, the Skater algorithm indicated 3-octanol, phellandral, thymol and d-limonene as those mostly influencing biofilm production inhibition (Figure 8 and Table 7), whose positive control was highlighted by means of partial dependence plots which describe the marginal impact of a feature on model prediction (Supplementary Material Figure SM-2).

Binary Classification Model for O-47 Biofilm Production Inhibition

The O-47/inhibition/3.125mg/mL optimized derived model was also obtained at cutoff of 62%, but with 19 PCs and a 0.51 actives:non-actives ratio (30:59) leading to 0.771, 0.590, 0.682 and 0.753 values for the Accuracy, MCC, Precision-Recall and ROC-AUC coefficients, respectively. Model robustness was assessed by AccuracyCV, MCCCV, Precision-RecallCV and ROC-AUCCV values of 0.738, 0.517, 0.589 and 0.659, correspondingly (Table 6). Y-scrambling application did not revealed the presence of any chance correlation (Table 7). Inspection of feature importance and partial dependence pointed out as more significant for biofilm production inhibition the compounds 3-octanol, o-cymene, d-limonene and β-phellandrene (Figure 9, Supplementary Material Figure SM-3 and Table 8).

Binary Classification Model for 25923 Biofilm Production Activation

The 25923/activation/3.125mg/mL optimized derived model was determined at cutoff of 121%, using 25 PCs, characterized by a 20:69 (0.29) proportion between actives and non-actives. The non cross-validated model was characterized by Accuracy, MCC, Precision-Recall and ROC-AUC coefficients of 0.906, 0.826, 0.956 and 0.961, respectively. Model robustness by cross-validation was characterized by high values of AccuracyCV = 0.763, MCCCV = 0.533, Precision-RecallCV = 0.824 and ROC-AUCCV = 0.834 (Table 6 and Table 7). Feature importance and partial dependence pointed out compounds d-limonene, γ-terpinene, 3-octanol and piperitenone as more important for biofilm production enhancement modulation (Figure 10, Supplementary Material Figure SM-4 and Table 8).

3. Discussion and Conclusions

3.1. EOs Biofilm Bioactivity General Consideration

From the results reported above it could be observed that each EO had a specific effect on biofilm formation, likely depending on its characteristics and unique chemical composition. In particular for S. aureus strains 6538P and 25923 the EOs mainly exhibited an enhancement of biofilm production. Stimulation of bacterial biofilm production by EOs is not surprising as it was previously observed, even by isolated chemical components [29,30,31]. On the other hand and more common [32,33], for S. epidermidis strains RP62A and O-47 an overall inhibition effect on biofilm production was observed by in vitro EOs treatment. Nevertheless, cinnamon EO was reported to stimulate biofilm production on some Staphylococcus epidermidis strains [34].

3.1.1. Bioactivity of RSEOs

The majority of tested RSEO samples did not show inhibitory effects on S. aureus 6538P biofilm formation (a partial inhibitory effect was observed only for R6 essential oil at 3.125 mg/mL, panel A of Figure 1). On the contrary, some RSEO samples (R6, R12, R24, RM4 and RM6) were shown to enhance biofilm formation by up to 140% at 0.0488 mg/mL Differently, several RSEO samples showed a good inhibitory effect on S. epidermids RP62A biofilm production. In particular 4 out of 13 EOs (R6, R24, RM4 and RM6) were able to potently inhibit biofilm formation with a rate of about 80% at either used concentrations (panel B of Figure 1), thus these EOs were selected for further analyses using scalar concentration of each EO starting from 0.0488 mg/mL. An attempt to determine a direct dose dependent effect was not effective (Figure 4).
On O-47 biofilm modulation (panel C of Figure 1) most RSEOs had a slight inhibition effect up to 40% (60% residual of biofilm production) at 0.0488 mg/mL. On the contrary, at the higher concentration RSEOs enhanced biofilm production by up to 130% for most samples. Only RM6 showed a remarkable biofilm production up to 160% at 3.125 mg/mL. For strain 25923 a profile similar to that of RP62A was observed (panel D of Figure 1), therefore no further investigation were pursued on the R6, R24, RM4 and RM6 samples despite the high inhibitory biofilm potency with residual biofilm production ranging 20–30%.

3.1.2. Bioactivity of FVEOs

Among all tested EOs, those from FV comprise among the most active samples able to inhibit biofilm production (Figure 2). In particular only mild effects (positive or negative modulation of biofilm production) were observed on 6538P strain (panel A of Figure 2) with a few exception at both selected concentrations, including FA24, FS2, FO3 and FO6 that increased biofilm production by about 40–60%, and FOM1 that inhibited biofilm production by about 50% at 3.125 mg/mL. On the contrary, some FVEOs proved to be potent antibiofilm agents on S. epidermids RP62A (panel B of Figure 2). In particular 5 out 33 FVEOs (FO1, FO3, FO6, FO24 and FOM3) inhibited biofilm formation with a rate of about 80% (panel B of Figure 2) and were selected for further analyses using scalar concentration of each EO starting from 0.0488 mg/mL. Similarly as for the selected potent RSEOs, no direct dose dependent effect was determined (Figure 5). Interestingly, on O-47 biofilm modulation most FVEOs showed a bioactivity profile almost overlapping that for RP62A with FO1, FO3, FO6, FO24 and FOM3 samples able to reduce biofilm formation of about 50–70% (panel C of Figure 2). While on both S. epidermidis strains RP62A and O-47, FVEOs displayed some peculiar samples with interesting biofilm inhibitory potencies in case of 25923 strain FVEOs displayed an overall bioactivity profile similar to that for RSEOs against O-47 (compare panel C of Figure 1 with panel D of Figure 2).

3.1.3. Bioactivity of CGEOs

EOs from CG are the most modulating biofilm producers either in positive (activators) or in negative (inhibitors). In particular in the case of strain 6538P all samples at either concentrations can be classified as neutral or biofilm promoters (panel A of Figure 3) with a strong inclination to increase biofilm production by up to 500% (CAM1). Other strong biofilm inducers (percentages over 300%) are CAM3, CAM5, CS1, CS3, CS6 and CS24 samples. Many other CGEOs, although to a lesser extent, induced a doubling or even tripling of biofilm production. On the contrary most of CGEOs displayed an inhibition by over 50% of biofilm production by RP62A (panel B of Figure 3). Many CGEO samples were further investigated and for 6 of them (CO2, CO6, COM5, CS2, CS6 and CSM5) a definite dose dependent relation was observed (Figure 6). Regarding biofilm modulation for O-47 CGEOs in this case presented, at either tested concentrations, a mixed scenario in which some samples induced an enhanced biofilm production up to 250-350% (CAM5, CS1, CS3, CS6, CSM1 and CSM3) and 15 different samples showed high inhibition potencies (percentages of residual biofilm lower than 40–50%).

3.2. Machine Learning Classification Models

Application of the PCA coupled with logistic regression led to the formulation of 4 robust models that were characterized by quite good Accuracy, MCC, Precision-Recall and ROC-AUC values (Table 7). Model agnostic feature importance and partial dependence plots were used to find the marginal effect that each EO chemical component has on the predicted outcome of the binary classification models built on the 3.125 mg/mL response variables. Feature importance is a measure of the prediction error of the model after the feature’s values are permuted and highlights the absolute importance of each chemical constituent while partial dependence plots show whether the relationship between the bioactivity and the chemical component is linear, monotonous or more complex.

3.2.1. Biofilm Activation ML Model on 6538P

Inspection of feature importance for model derived on 6538P biofilm percentage production and EOs’ chemical compositions revealed 3-octanol, d-limonene and pulegone as the chemical components more associated to bacterial biofilm production (Figure 7 and Table 8). Further investigation of their partial dependence plots (Supplementary Material Figure SM-2) indicated those three chemicals as all positively correlated with biofilm enhancement.

3.2.2. Biofilm Activation ML Model on 25923

Similarly as for 6538P, also for the 25923 strain a ML model was built to correlate biofilm production enhancement with EOs’ chemical composition. Again, analysis of feature importance found as more important d-limonene, γ-terpinene, 3-octanol and piperitenone (Figure 8). Differently as found for 6538P the main component were not all positively associated to biofilm enhancement production, but d-limonene, γ-terpinene, 3-octanol were suggested to negatively modulate the increase of biofilm, while piperitenone was found positively correlated (Table 8 and Supplementary Material Figure SM-3)

3.2.3. Biofilm Inhibition ML Model on RP62A

Differently from the previous model, feature importance associated to the EOs biofilm inhibition production on RP62A strain highlighted 3-octanol, phellandral, thymol, d-limonene as chemical compounds most important on modulating biofilm reduction (Figure 9 and Table 8). Partial dependence plots for 3-octanol, phellandral, thymol, d-limonene associated the four components as all positively able to inhibit biofilm production (Supplementary Material Figure SM-4).

3.2.4. Biofilm Inhibition ML Model on O-47

Regarding ML model derived on the biofilm inhibition capability of EOs the compounds more responsible for biofilm production modulation were found to be 3-octanol, o-cymene, d-limonene and β-phellandrene (Figure 10 and Table 8). Differently from above RP62A analogous inhibition model only 3-octanol and d-limonene were found positively associated with EOs’ inhibitory ability by partial dependence plots (Supplementary Material Figure SM-5). On the contrary o-cymene and β-phellandrene were associated to a negative action on the inhibition. This could be speculated as a sort of anti-synergic effect that could balance EOs’ potencies.

3.2.5. General Consideration on ML Models

According to the four classification models, two compounds, namely 3-octanol and d-limonene, can be considered as those that most influence biofilm production (Table 8). In particular, d-limonene positively correlated either in inhibiting or in enhancement of biofilm production in three out of the four models while has a negative modulation on the ML model built on the biofilm enhancement of EOs’ on 25923 strain. These data indicate some controversial mechanism associated to d-limonene, it could be speculated that being this compound a highly apolar monoterpene its role could not be indirectly associated to biofilm modulation by altering the bacterial wall [35] allowing other compounds, likely oxygenated ones to enter the cell acting in altering some biochemical mechanism that could end in stimulation or inhibition of biofilm production. Nevertheless, on this topic the data available in the literature is controversial: Natcha and Caoili [36] reported that d-limonene is effective in inhibiting the growth of S. epidermidis RP62A when combined with the antibiotic rifampicin, likely due to d-limonene interference with biofilm formation. The effect of d-limonene in inhibiting bacterial biofilm formation was also proved against species of the genus Streptococcus [37] for which minimal biofilm inhibitory concentration (MBIC) of 400 μg/mL was determined. In a very recent study d-limonene was also reported as a biofilm inhibitor, although less efficient than an EO containing d-limonene [38]. On the contrary Kerekes et al. assayed a series of EOs and a list of chemical components against food-related micro-organisms and found d-limonene was almost deprived of any ability to inhibit biofilm production. In a study from Espina et al. d-limonene at 2000 μL/L was reported to reduce the production of biofilm mass in S. aureus USA300 by 90% after 8 h of incubation, but increase it by 30% after 40 h of incubation [39]. EOs containing d-limonene and the isolated component were found to stimulate biofilm production on Listeria monocytogenes and antibiotic-resistant Enterococcus faecalis strains [29,30,31]. A similar profile and speculation on 3-octanol could also be deduced. 3-Octanol is a molecule resembling normal octanol, a compound commonly used to evaluate compound membrane permeability and lipophilicity through the determination of the logP parameter often used in ADME and QSAR studies. Unfortunately no data are available on the influence of 3-octanol on biofilm production, except for a single report in which the 8-carbon molecules 1-octen-3-ol, 3-octanol and 3-octanone specifically induced conidiation in Trichoderma species colonies placed in the dark [40]. Considering the possible cell wall permeation role of both d-limonene and 3-octanol for strains 6538P and 25923 pulegone, γ-terpinene and piperitenone, on the basis of the ML elaboration, could be the main components responsible for the modulation of EOs’ augmented biofilm production, while in the case RP62A and O-47 phellandral, thymol, o-cymene and β-phellandrene are mainly responsible for positively (phellandral and thymol) or negatively (o-cymene and β-phellandrene) modulating EOs’ biofilm inhibition. Unfortunately no specific data data are available on these isolated components and the herein discussion although based on robust ML calculation are not experimentally based. It is worthy to note that the four bacterial strains tested here produced biofilms with different characteristics. First 6538P and 25923 belong to S. aureus species, while RP62A and O47 belong to S. epidermidis species. 25923 is classified as a strong biofilm producer, and 6538P is a medium/strong biofilm producer according to Cafiso and coworkers [41]. Proteins are the major component in the biofilm matrix of 6538P, while in 25923 the polysaccharides have a predominant role. As regards the S. epidermidis strains, they are both strong biofilm makers and produce a biofilm mainly composed by polysaccharides. Moreover O-47 is a naturally occurring agr mutant [42]. As previously reported [43], agr-negative genotype enhanced biofilm formation on polymer surfaces by an increased expression of the surface protein AtlE, a bifunctional adhesin/autolysin abundant in the cell wall of S. epidermidis. The amount of AtlE present in the cell envelop is one of the reported differences between RP62A and O-47 [43]. The overexpression of AtlE could induce significant changes in the hydrophobicity of the bacterial surface [44]; this effect could explain the different action of EOs on these two strains. Furthermore, the classification models were developed on the same EOs tested on P. aeruginosa biofilm production [22]. In that case, investigation of the most important components by means of feature importance and partial dependence plots indicated estragole and phellandral as the chemical components mostly related to biofilm inhibition of P. aeruginosa, while d-limonene, pulegone and chrysanthenone seem to be related to its biofilm production. Although the use of feature importance and partial dependence plots shed some light on the possible role of some EOs’ components little is yet known on the role of the whole EOs mixture synergisms and anti-synergisms. Further studies on isolated EOs’ chemical components and on their simple mixture are currently under evaluation to develop more refined ML models able to disclose more details on the EOs’ mechanism of action.

4. Materials and Methods

4.1. Essential oil and Chemical Composition Analysis

EOs and their chemical compositions were available from previously reported studies [22,45,46]. Briefly, EOs were obtained by direct fractionated steam distillation and analyzed by a gas chromatographic/mass spectrometric (GC/MS) protocol [47,48].

4.2. Bacterial Strains and Culture Conditions

Bacterial strains used in this work (Table 9) were grown in Brain Heart Infusion broth (BHI, Oxoid, UK). Biofilm formation was assessed in static conditions. Planktonic cultures were grown in flasks under vigorous agitation (180 rpm) at 37 °C. In particular, S. aureus ATCC 6538P (6538P) and S. aureus ATCC 25923 (25923) are reference strains for antimicrobial testing; S. epidermidis RP62A (RP62A) is a reference strain isolated from infected catheter, while S. epidermidis O-47 (O-47) is a clinical isolate strong biofilm producer strain characterized by a genomic mutation in agr locus [42].

4.3. Determination of Minimal Inhibitory Concentration (MIC)

The MIC was determined as the lowest concentration at which the observable bacterial growth was inhibited. MICs were determined according to the guidelines of Clinical Laboratory Standards Institute (CLSI [49]). Each EO was added directly from mother stock and solutions were prepared by two-fold serial dilutions. Mother stock solutions were obtained by solubilizing each EO in DMSO at a final concentration of 1 g/mL. Appropriate dilution (106 cfu/mL) of bacterial culture in exponential phase was used. Ten concentrations were used within the 25–0.045 mg/mL range. Experiments were performed in quadruplicate.

4.4. Biofilm Production Assay

The quantification of biofilm production was based on microtiter plate biofilm assay (MTP): an opportune dilution of bacterial culture in exponential growth phase was added into wells of a sterile 96-well flat-bottomed polystyrene plate in absence and in presence of each EO. Quantification of in vitro biofilm production was based on previously reported methodology [50]. Briefly, the wells of a sterile 96-well flat-bottomed polystyrene plate were filled with 100 µL of the appropriate medium. 1/100 dilution of overnight bacterial cultures was added into each well (about 0.5 OD 600nm). As control, the first row contained bacteria grown in 100 μL of BHI (untreated bacteria). In the second row was added BHI supplemented with each EO at concentrations of 3.125 mg/mL and 0.0488 mg/mL, respectively. The plates were incubated aerobically for 18 h at 37 °C. Biofilm formation was measured using crystal violet staining. After treatment, planktonic cells were gently removed; each well was washed three times with double-distilled water and patted dry with a piece of paper towel in an inverted position. To quantify biofilm formation, each well was stained with 0.1% crystal violet and incubated for 15 min at room temperature, rinsed twice with double-distilled water, and thoroughly dried. The dye bound to adherent cells was solubilized with 20% (v/v) glacial acetic acid and 80% (v/v) ethanol. After 30 min of incubation at room temperature, OD590 was measured to quantify the total biomass of biofilm formed in each well. Each data point is composed of 4 independent experiments, each performed at least in 3-replicates. EOs altering biofilm formation of selected strains were then tested as reported below. Briefly, the wells of a sterile 96-well flat-bottomed polystyrene plate were filled with 100 µL of the appropriate medium. 1/100 dilution of overnight bacterial cultures was added into each well (about 0.5 OD 600nm). As control, the first row contained bacteria grown in 100 μL of BHI (untreated bacteria). Furthermore, BHI broth was added to remaining wells starting from the third row. In the second row was added BHI supplemented with each EO at a concentration of 0.0488 mg/mL. Starting from this lane, samples were serially diluted (1:2 dilutions). The plates were incubated aerobically for 18 h at 37 °C. Biofilm formation was measured using crystal violet staining, as previously reported.

4.5. Statistical Analysis of Biological Evaluation

Data reported were statistically validated using Student’s t-test comparing mean absorbance of treated and untreated samples. The significance of differences between mean absorbance values was calculated using a two-tailed Student’s t-test. A p value of <0.05 was considered significant.

4.6. Machine Learning Binary Classification

4.6.1. General Methods

All calculations were performed using the Python (version 3.6, https://www.python.org/) programming language [51] by executing in-house code in the Jupyter Notebook platform (version 5.7) [52]. The datasets were imported and loaded into a Pandas [53] dataframe and pre-processed to obtain four independent data matrices consisting of 89 rows (essential oil samples) and 54 columns (chemical components). Two dependent target vectors containing 89 biofilm production percentage observations at 48 μg/mL and 3.125 mg/mL were defined. Machine learning algorithms used in this study were implemented using the sklearn library (version 0.20) [54]. Unsupervised dimensionality reduction was performed with Principal component analysis (PCA) [55] while L2 regularized logistic regression was used for the supervised learning analysis. The scores and loadings relatives to the first two principal components (PCs) were graphically inspected on plots generated using the matplotlib library (version 3.0) [56]. To build the classification models, 30 PCs were extracted for each dataset. Cross-validation was used to search for the optimal inhibition/activation percentage cut-off values in order to define active and inactive samples. The optimal cut-off values were used to obtain the hyper-parameters optimized classification models. The Hyper-parameters optimization was achieved through a Bayesian optimization [57] of the number of PCs to be used as features and the regularization parameter of the L2-Logistic Regression (inverse of regularization strength in the sklearn implementation). For each dependent target vector two types of models were built: one to define EOs ability to inhibit biofilm production and another to describe biofilm production enhancement. Percentage ranges of 60–80% and 120–140% biofilm productions were chosen for inhibition and activation models, respectively. Finally, the most appropriate cut-offs for binary classification of biofilm inhibitors/not-inhibitors or biofilm enhancers/not-enhancers EOs were determined from a supervised learning analysis.
The binary classification models were numerically and graphically evaluated by accuracy (ACC), Matthews correlation coefficient (MCC), receiver operating characteristic (ROC) and precision-recall (PR) curves. Finally, the importance of EOs chemical components was evaluated individually through the “feature importance” and “partial dependence” plots [28] as implemented in the Skater python library [58,59]. Feature importance is a generic term for the degree to which a predictive model relies on a particular feature. Skater feature importance implementation is based on an information theoretic criteria, measuring the entropy in the change of predictions, given a perturbation of a given feature [58].

4.6.2. Classification Models’ Validation

Validation of each classification model was carried out by leave-one-out cross-validation and taking into account the accuracy (ACC), the precision or positive predictive value (PPV), the recall or sensitivity or true positive rate (TPR), specificity or true negative rate (TNR), receiver operating characteristic (ROC) curve and the Matthews correlation coefficient (MCC) (see Supplementary Information) [22,60]. Y-scrambling [61,62] was ultimately applied to check any lack of chance correlation and assess coefficients robustness.

Supplementary Materials

The following are available online. Figure SM-1: PCA first 2 PCs graphical plots. Scores plot (panel A) indicate the presence of at least three clusters (circled in panel A). Loading plots (panel B) highlights that estragole, o-cymene and pulegone could be the most important chemical constituents among all the tested EOs, Figure SM-2: 3-octanol, d-limonene and pulegone partial dependence plots for the activation model on 6538P biofilm production, Figure SM-3: d-limonene, γ-terpinene, 3-octanol and piperitenone partial dependence plots for the inhibition model on 25923 biofilm production, Figure SM-4: 3-octanol, phellandral, thymol and d-limonene partial dependence plots for the inhibition model on RP62A biofilm, Figure SM-5: 3-octanol, o-cymene, d-limonene and β-phellandrene partial dependence plots for the inhibition model on O-47 biofilm.

Author Contributions

Conceptualization, S.M., L.S. and R.R.; Data curation, M.B.; Formal analysis, R.B.; Funding acquisition, S.M.; Investigation, A.P., M.A., R.P., M.S., S.G. and G.V.; Supervision, S.M.; Validation, L.S.; Writing—original draft, R.R.; Writing—review & editing, R.R., S.M. and L.S.

Funding

This work was supported by the University of Ferrara (Bando FAR 2018) to SM.

Acknowledgments

R.R. and A.P. thanks Alchemical Dynamics s.r.l. for the technical support.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Donlan, R.M.; Costerton, J.W. Biofilms: Survival mechanisms of clinically relevant microorganisms. Clin. Microbiol. Rev. 2002, 15, 167–193. [Google Scholar] [CrossRef] [PubMed]
  2. Da Costa Krewer, C.; Santos Amanso, E.; Veneroni Gouveia, G.; de Lima Souza, R.; da Costa, M.M.; Aparecido Mota, R. Resistance to antimicrobials and biofilm formation in Staphylococcus spp. isolated from bovine mastitis in the Northeast of Brazil. Trop Anim. Health Prod. 2015, 47, 511–518. [Google Scholar] [CrossRef] [PubMed]
  3. Furuya, E.Y.; Lowy, F.D. Antimicrobial-resistant bacteria in the community setting. Nat. Rev. Microbiol. 2006, 4, 36–45. [Google Scholar] [CrossRef] [PubMed]
  4. Satpathy, S.; Sen, S.K.; Pattanaik, S.; Raut, S. Review on bacterial biofilm: An universal cause of contamination. Biocatal. Agric. Biotech. 2016, 7, 56–66. [Google Scholar] [CrossRef]
  5. Koo, H.; Allan, R.N.; Howlin, R.P.; Stoodley, P.; Hall-Stoodley, L. Targeting microbial biofilms: Current and prospective therapeutic strategies. Nat. Rev. Microbiol. 2017, 15, 740–755. [Google Scholar] [CrossRef] [PubMed]
  6. Otto, M. Staphylococcal biofilms. Curr. Top Microbiol. Immunol. 2008, 322, 207–228. [Google Scholar] [PubMed]
  7. Li, Z. A review of Staphylococcus aureus and the emergence of drug-resistant problem. Adv. Microbiol. 2018, 8, 65–76. [Google Scholar] [CrossRef]
  8. Nicholson, T.L.; Shore, S.M.; Smith, T.C.; Frana, T.S. Livestock-associated methicillin-resistant Staphylococcus aureus (LA-MRSA) isolates of swine origin form robust biofilms. PLoS ONE 2013, 8, e73376. [Google Scholar] [CrossRef]
  9. Dohar, J.E.; Hebda, P.A.; Veeh, R.; Awad, M.; Costerton, J.W.; Hayes, J.; Ehrlich, G.D. Mucosal biofilm formation on middle-ear mucosa in a nonhuman primate model of chronic suppurative otitis media. Laryngoscope 2005, 115, 1469–1472. [Google Scholar] [CrossRef] [PubMed]
  10. Rogers, K.L.; Fey, P.D.; Rupp, M.E. Coagulase-negative staphylococcal infections. Infect. Dis. Clin. North Am. 2009, 23, 73–98. [Google Scholar] [CrossRef] [PubMed]
  11. Wang, Z.; Bai, H.; Lu, C.; Hou, C.; Qiu, Y.; Zhang, P.; Duan, J.; Mu, H. Light controllable chitosan micelles with ROS generation and essential oil release for the treatment of bacterial biofilm. Carbohydr. Polym. 2019, 205, 533–539. [Google Scholar] [CrossRef] [PubMed]
  12. Zhao, X.; Liu, Z.; Liu, Z.; Meng, R.; Shi, C.; Chen, X.; Bu, X.; Guo, N. Phenotype and RNA-seq-based transcriptome profiling of Staphylococcus aureus biofilms in response to tea tree oil. Microb. Pathog. 2018, 123, 304–313. [Google Scholar] [CrossRef] [PubMed]
  13. Vazquez-Sanchez, D.; Galvao, J.A.; Mazine, M.R.; Gloria, E.M.; Oetterer, M. Control of Staphylococcus aureus biofilms by the application of single and combined treatments based in plant essential oils. Int. J. Food Microbiol. 2018, 286, 128–138. [Google Scholar] [CrossRef] [PubMed]
  14. Vaillancourt, K.; LeBel, G.; Yi, L.; Grenier, D. In vitro antibacterial activity of plant essential oils against Staphylococcus hyicus and Staphylococcus aureus, the causative agents of exudative epidermitis in pigs. Arch. Microbiol. 2018, 200, 1001–1007. [Google Scholar] [CrossRef] [PubMed]
  15. Scaffaro, R.; Lopresti, F.; D’Arrigo, M.; Marino, A.; Nostro, A. Efficacy of poly(lactic acid)/carvacrol electrospun membranes against Staphylococcus aureus and Candida albicans in single and mixed cultures. Appl. Microbiol. Biotechnol. 2018, 102, 4171–4181. [Google Scholar] [CrossRef] [PubMed]
  16. Merghni, A.; Noumi, E.; Hadded, O.; Dridi, N.; Panwar, H.; Ceylan, O.; Mastouri, M.; Snoussi, M. Assessment of the antibiofilm and antiquorum sensing activities of Eucalyptus globulus essential oil and its main component 1,8-cineole against methicillin-resistant Staphylococcus aureus strains. Microb. Pathog. 2018, 118, 74–80. [Google Scholar] [CrossRef] [PubMed]
  17. Kot, B.; Wierzchowska, K.; Gruzewska, A.; Lohinau, D. The effects of selected phytochemicals on biofilm formed by five methicillin-resistant Staphylococcus aureus. Nat. Prod. Res. 2018, 32, 1299–1302. [Google Scholar] [CrossRef] [PubMed]
  18. Chovanova, R.; Mezovska, J.; Vaverkova, S.; Mikulasova, M. The inhibition the Tet(K) efflux pump of tetracycline resistant Staphylococcus epidermidis by essential oils from three Salvia species. Lett. Appl. Microbiol. 2015, 61, 58–62. [Google Scholar] [CrossRef] [PubMed]
  19. Chovanova, R.; Mikulasova, M.; Vaverkova, S. In vitro antibacterial and antibiotic resistance modifying effect of bioactive plant extracts on methicillin-resistant Staphylococcus epidermidis. Int. J. Microbiol. 2013, 2013, 7. [Google Scholar] [CrossRef] [PubMed]
  20. Nostro, A.; Sudano Roccaro, A.; Bisignano, G.; Marino, A.; Cannatelli, M.A.; Pizzimenti, F.C.; Cioni, P.L.; Procopio, F.; Blanco, A.R. Effects of oregano, carvacrol and thymol on Staphylococcus aureus and Staphylococcus epidermidis biofilms. J. Med. Microbiol. 2007, 56, 519–523. [Google Scholar] [CrossRef] [PubMed]
  21. Taghadomi-Saberi, S.; Mas Garcia, S.; Allah Masoumi, A.; Sadeghi, M.; Marco, S. Classification of bitter orange essential oils according to fruit ripening stage by untargeted chemical profiling and machine learning. Sensors 2018, 18, 14. [Google Scholar] [CrossRef] [PubMed]
  22. Artini, M.; Patsilinakos, A.; Papa, R.; Bozovic, M.; Sabatino, M.; Garzoli, S.; Vrenna, G.; Tilotta, M.; Pepi, F.; Ragno, R.; et al. Antimicrobial and antibiofilm activity and machine learning classification analysis of essential oils from different mediterranean plants against pseudomonas aeruginosa. Molecules 2018, 23, 13. [Google Scholar] [CrossRef] [PubMed]
  23. Riahi, S.; Pourbasheer, E.; Ganjali, M.R.; Norouzi, P. Investigation of different linear and nonlinear chemometric methods for modeling of retention index of essential oil components: Concerns to support vector machine. J. Hazard Mater. 2009, 166, 853–859. [Google Scholar] [CrossRef] [PubMed]
  24. Drevinskas, T.; Mickienė, R.; Maruška, A.; Stankevičius, M.; Tiso, N.; Šalomskas, A.; Lelešius, R.; Karpovaitė, A.; Ragažinskienė, O. Confirmation of the antiviral properties of medicinal plants via chemical analysis, machine learning methods and antiviral tests: A methodological approach. Anal. Methods 2018, 10, 1875–1885. [Google Scholar] [CrossRef]
  25. Cortes, C.; Vapnik, V. Support-vector networks. Mach. Learn. 1995, 20, 273–297. [Google Scholar] [CrossRef]
  26. Ho, T.K. The random subspace method for constructing decision forests. IEEE Trans. Pattern Anal. 1998, 20, 832–844. [Google Scholar]
  27. Vert, J.-P. Kernel methods in genomics and computational biology. arXiv, 2005; arXiv:q-bio/0510032. [Google Scholar]
  28. Friedman, J.H. Greedy function approximation: A gradient boosting machine. Ann. Stat. 2001, 29, 1189–1232. [Google Scholar] [CrossRef]
  29. Sandasi, M.; Leonard, C.M.; Viljoen, A.M. The in vitro antibiofilm activity of selected culinary herbs and medicinal plants against Listeria monocytogenes. Lett. Appl. Microbiol. 2010, 50, 30–35. [Google Scholar] [CrossRef] [PubMed]
  30. Sandasi, M.; Leonard, C.M.; Viljoen, A.M. The effect of five common essential oil components on Listeria monocytogenes biofilms. Food Control 2008, 19, 1070–1075. [Google Scholar] [CrossRef]
  31. Negreiros, M.O.; Pawlowski, A.; Zini, C.A.; Soares, G.L.; Motta, A.S.; Frazzon, A.P. Antimicrobial and antibiofilm activity of Baccharis psiadioides essential oil against antibiotic-resistant Enterococcus faecalis strains. Pharm. Biol. 2016, 54, 3272–3279. [Google Scholar] [CrossRef] [PubMed]
  32. Szczepanski, S.; Lipski, A. Essential oils show specific inhibiting effects on bacterial biofilm formation. Food Control. 2014, 36, 224–229. [Google Scholar] [CrossRef]
  33. Kannappan, A.; Sivaranjani, M.; Srinivasan, R.; Rathna, J.; Pandian, S.K.; Ravi, A.V. Inhibitory efficacy of geraniol on biofilm formation and development of adaptive resistance in Staphylococcus epidermidis RP62A. J. Med. Microbiol. 2017, 66, 1506–1515. [Google Scholar] [CrossRef] [PubMed]
  34. Nuryastuti, T.; van der Mei, H.C.; Busscher, H.J.; Iravati, S.; Aman, A.T.; Krom, B.P. Effect of cinnamon oil on icaA expression and biofilm formation by Staphylococcus epidermidis. Appl. Environ. Microbiol. 2009, 75, 6850–6855. [Google Scholar] [CrossRef] [PubMed]
  35. Chueca, B.; Pagan, R.; Garcia-Gonzalo, D. Differential mechanism of Escherichia coli inactivation by (+)-limonene as a function of cell physiological state and drug’s concentration. PLoS ONE 2014, 9, e94072. [Google Scholar] [CrossRef] [PubMed]
  36. Rummaneethorn, N.; Caoili, C.M. d-limonene Increases Efficacy of Rifampicin as an Inhibitor of In Vitro Growth of Opportunistic Staphylococcus epidermidis RP62A. In Proceedings of the National Conference On Undergraduate Research (NCUR), Asheville, NC, USA, 7–9 April 2016; University of North Carolina at Asheville: Asheville, NC, USA. [Google Scholar]
  37. Subramenium, G.A.; Vijayakumar, K.; Pandian, S.K. Limonene inhibits streptococcal biofilm formation by targeting surface-associated virulence factors. J. Med. Microbiol. 2015, 64, 879–890. [Google Scholar] [CrossRef] [PubMed]
  38. Cerioli, M.F.; Moliva, M.V.; Cariddi, L.N.; Reinoso, E.B. Effect of the essential oil of minthostachys verticillata (Griseb.) epling and limonene on biofilm production in pathogens causing bovine mastitis. Front. Vet. Sci. 2018, 5, 146. [Google Scholar] [CrossRef] [PubMed]
  39. Espina, L.; Pagan, R.; Lopez, D.; Garcia-Gonzalo, D. Individual constituents from essential oils inhibit biofilm mass production by multi-drug resistant Staphylococcus aureus. Molecules 2015, 20, 11357–11372. [Google Scholar] [CrossRef] [PubMed]
  40. Majumdar, S.; Mondal, S. Perspectives on quorum sensing in Fungi. bioRxiv 2015, 019034. [Google Scholar] [CrossRef]
  41. Cafiso, V.; Bertuccio, T.; Santagati, M.; Demelio, V.; Spina, D.; Nicoletti, G.; Stefani, S. agr-Genotyping and transcriptional analysis of biofilm-producing Staphylococcus aureus. FEMS Immunol. Med. Microbiol. 2007, 51, 220–227. [Google Scholar] [CrossRef] [PubMed]
  42. Heilmann, C.; Gerke, C.; Perdreau-Remington, F.; Gotz, F. Characterization of Tn917 insertion mutants of Staphylococcus epidermidis affected in biofilm formation. Infect. Immun. 1996, 64, 277–282. [Google Scholar] [PubMed]
  43. Vuong, C.; Gerke, C.; Somerville, G.A.; Fischer, E.R.; Otto, M. Quorum-sensing control of biofilm factors in Staphylococcus epidermidis. J. Infect. Dis. 2003, 188, 706–718. [Google Scholar] [CrossRef] [PubMed]
  44. Buttner, H.; Mack, D.; Rohde, H. Structural basis of Staphylococcus epidermidis biofilm formation: Mechanisms and molecular interactions. Front. Cell Infect. Microbiol. 2015, 5, 14. [Google Scholar] [PubMed]
  45. Garzoli, S.; Bozovic, M.; Baldisserotto, A.; Sabatino, M.; Cesa, S.; Pepi, F.; Vicentini, C.B.; Manfredini, S.; Ragno, R. Essential oil extraction, chemical analysis and anti-Candida activity of Foeniculum vulgare Miller—new approaches. Nat. Prod. Res. 2017, 32, 1–6. [Google Scholar] [CrossRef] [PubMed]
  46. Bozovic, M.; Garzoli, S.; Sabatino, M.; Pepi, F.; Baldisserotto, A.; Andreotti, E.; Romagnoli, C.; Mai, A.; Manfredini, S.; Ragno, R. Essential oil extraction, chemical analysis and anti-Candida Activity of Calamintha nepeta (L.) Savi subsp. glandulosa (Req.) ball-new approaches. Molecules 2017, 22, 13. [Google Scholar] [CrossRef] [PubMed]
  47. Bozovic, M.; Navarra, A.; Garzoli, S.; Pepi, F.; Ragno, R. Esential oils extraction: A 24-h steam distillation systematic methodology. Nat. Prod. Res. 2017, 31, 2387–2396. [Google Scholar] [CrossRef] [PubMed]
  48. Garzoli, S.; Pirolli, A.; Vavala, E.; Di Sotto, A.; Sartorelli, G.; Bozovic, M.; Angiolella, L.; Mazzanti, G.; Pepi, F.; Ragno, R. Multidisciplinary approach to determine the optimal time and period for extracting the essential oil from mentha suaveolens Ehrh. Molecules 2015, 20, 9640–9655. [Google Scholar] [CrossRef] [PubMed]
  49. Humphries, R.M.; Ambler, J.; Mitchell, S.L.; Castanheira, M.; Dingle, T.; Hindler, J.A.; Koeth, L.; Sei, K.; Development, C.M.; Standardization Working Group of the Subcommittee on Antimicrobial Susceptibility Testing. CLSI methods development and standardization working group best practices for evaluation of antimicrobial susceptibility tests. J. Clin. Microbiol. 2018, 56, 10. [Google Scholar] [CrossRef] [PubMed]
  50. Papa, R.; Selan, L.; Parrilli, E.; Tilotta, M.; Sannino, F.; Feller, G.; Tutino, M.L.; Artini, M. Anti-biofilm activities from marine cold adapted bacteria against Staphylococci and Pseudomonas aeruginosa. Front. Microbiol. 2015, 6, 10. [Google Scholar] [CrossRef] [PubMed]
  51. Perkel, J.M. Pick up Python. Nature 2015, 518, 125–126. [Google Scholar] [CrossRef] [PubMed]
  52. Kluyver, T.; Ragan-Kelley, B.; Pérez, F.; Granger, B.; Bussonnier, M.; Frederic, J.; Kelley, K.; Hamrick, J.; Grout, J.; Corlay, S.; et al. Jupyter Notebooks—A publishing format for reproducible computational workflows. In Positioning and Power in Academic Publishing: Players, Agents and Agendas; Loizides, F., Scmidt, B., Eds.; IOS Press: Amsterdam, The Netherlands, 2016; pp. 87–90. [Google Scholar] [CrossRef]
  53. McKinney, W. Data Structures for Statistical Computing in Python. In Proceedings of the Python for Scientific Computing, Astin, TX, USA, 28 June–3 July 2010; pp. 51–56. Available online: https://conference.scipy.org/scipy2010/ (accessed on 2 March 2019).
  54. Pedregosa, F.; Varoquaux, G.; Gramfort, A.; Michel, V.; Thirion, B.; Grisel, O.; Blondel, M.; Prettenhofer, P.; Weiss, R.; Dubourg, V.; et al. Scikit-learn: Machine learning in Python. J. Mach. Learn. Res. 2011, 12, 2825–2830. [Google Scholar]
  55. Tipping, M.E.; Bishop, C.M. Probabilistic principal component analysis. J. R. Stat. Soc. B 1999, 61, 611–622. [Google Scholar] [CrossRef]
  56. Hunter, J.D. Matplotlib: A 2D graphics environment. Comput. Sci. Eng. 2007, 9, 90–95. [Google Scholar] [CrossRef]
  57. Snoek, J.; Larochelle, H.; Adams, R.P. Practical Bayesian optimization of machine learning algorithms. In Proceedings of the 25th International Conference on Neural Information Processing Systems—Volume 2; Curran Associates Inc.: Lake Tahoe, NV, USA, 2012; pp. 2951–2959. [Google Scholar]
  58. Choudhary, P.; Kramer, A.; datascience.com team. Datascienceinc/skater: Enable interpretability via rule extraction(BRL). Zenodo, 2018. Available online: https://zenodo.org/record/1198885#.XHsqergRXIU (accessed on 2 March 2019).
  59. Wei, P.; Lu, Z.; Song, J. Variable importance analysis: A comprehensive review. Reliab. Eng. Syst. Saf. 2015, 142, 399–432. [Google Scholar] [CrossRef]
  60. Baldi, P.; Brunak, S.; Chauvin, Y.; Andersen, C.A.F.; Nielsen, H. Assessing the accuracy of prediction algorithms for classification: An overview. Bioinformatics 2000, 16, 412–424. [Google Scholar] [CrossRef] [PubMed]
  61. Klopman, G.; Kalos, A.N. Causality in Structure Activity Studies. J. Comput. Chem. 1985, 6, 492–506. [Google Scholar] [CrossRef]
  62. Rucker, C.; Rucker, G.; Meringer, M. Y-Randomization and its variants in QSPR/QSAR. J. Chem. Inf. Model. 2007, 47, 2345–2357. [Google Scholar] [CrossRef] [PubMed]
Sample Availability: Samples of the compounds are not available from the authors.
Figure 1. Percentages of biofilm production after treatment at two concentrations (3.125 mg/mL and 0.0488 mg/mL) for RSEOs against the four strains S. aureus 6538P (A) and 25923 (B), S. epidermidis RP62A (C) and O-47 (D, respectively). In the ordinate axis are reported the percentage of bacterial biofilm production. Data are reported as percentage of residual biofilm after the treatment in comparison with the untreated one. Each data point is composed of 4 independent experiments each performed with at least three replicates.
Figure 1. Percentages of biofilm production after treatment at two concentrations (3.125 mg/mL and 0.0488 mg/mL) for RSEOs against the four strains S. aureus 6538P (A) and 25923 (B), S. epidermidis RP62A (C) and O-47 (D, respectively). In the ordinate axis are reported the percentage of bacterial biofilm production. Data are reported as percentage of residual biofilm after the treatment in comparison with the untreated one. Each data point is composed of 4 independent experiments each performed with at least three replicates.
Molecules 24 00890 g001aMolecules 24 00890 g001b
Figure 2. Percentages of biofilm production after treatment at two concentrations (3.125 mg/mL and 0.0488 mg/mL) for FVEOs against the four strains S. aureus 6538P (A) and 25923 (B), S. epidermidis RP62A (C) and O-47 (D, respectively). In the ordinate axis are is reported the percentage of bacterial biofilm production. Data are reported as percentage of residual biofilm after the treatment in comparison with the untreated one. Each data point is composed of 4 independent experiments each performed with at least three replicates.
Figure 2. Percentages of biofilm production after treatment at two concentrations (3.125 mg/mL and 0.0488 mg/mL) for FVEOs against the four strains S. aureus 6538P (A) and 25923 (B), S. epidermidis RP62A (C) and O-47 (D, respectively). In the ordinate axis are is reported the percentage of bacterial biofilm production. Data are reported as percentage of residual biofilm after the treatment in comparison with the untreated one. Each data point is composed of 4 independent experiments each performed with at least three replicates.
Molecules 24 00890 g002aMolecules 24 00890 g002b
Figure 3. Percentages of biofilm production after treatment at two concentrations (3.125 mg/mL and 0.0488 mg/mL) for CGEOs against the four strains S. aureus 6538P (A) and 25923 (B), S. epidermidis RP62A (C) and O-47 (D), respectively). In the ordinate axis are is reported the percentage of bacterial biofilm production. The abscissa axis is centered at 100% biofilm production. Data are reported as percentage of residual biofilm after the treatment in comparison with the untreated one. Each data point is composed of 4 independent experiments each performed with at least three replicates.
Figure 3. Percentages of biofilm production after treatment at two concentrations (3.125 mg/mL and 0.0488 mg/mL) for CGEOs against the four strains S. aureus 6538P (A) and 25923 (B), S. epidermidis RP62A (C) and O-47 (D), respectively). In the ordinate axis are is reported the percentage of bacterial biofilm production. The abscissa axis is centered at 100% biofilm production. Data are reported as percentage of residual biofilm after the treatment in comparison with the untreated one. Each data point is composed of 4 independent experiments each performed with at least three replicates.
Molecules 24 00890 g003aMolecules 24 00890 g003b
Figure 4. Antibiofilm effect of selected RSEOs on RP62A and on O-47 strains. In the ordinate axis is reported the percentage of bacterial biofilm production. Data are reported as percentage of residual biofilm after the treatment in comparison with the untreated one. Each data point is composed of 4 independent experiments each performed with at least in three replicates.
Figure 4. Antibiofilm effect of selected RSEOs on RP62A and on O-47 strains. In the ordinate axis is reported the percentage of bacterial biofilm production. Data are reported as percentage of residual biofilm after the treatment in comparison with the untreated one. Each data point is composed of 4 independent experiments each performed with at least in three replicates.
Molecules 24 00890 g004
Figure 5. Antibiofilm effect of selected FVEOs on RP62A and on O-47 strains. In the ordinate axis is reported the percentage of bacterial biofilm production. Data are reported as percentage of residual biofilm after the treatment in comparison with the untreated one. Each data point is composed of 4 independent experiments each performed with at least three replicates.
Figure 5. Antibiofilm effect of selected FVEOs on RP62A and on O-47 strains. In the ordinate axis is reported the percentage of bacterial biofilm production. Data are reported as percentage of residual biofilm after the treatment in comparison with the untreated one. Each data point is composed of 4 independent experiments each performed with at least three replicates.
Molecules 24 00890 g005
Figure 6. Antibiofilm effect of selected CGEOs on RP62A and on O-47 strains. In the ordinate axis is reported the percentage of bacterial biofilm production. Data are reported as percentage of residual biofilm after the treatment in comparison with the untreated one. Each data point is composed of 4 independent experiments each performed with at least three replicates.
Figure 6. Antibiofilm effect of selected CGEOs on RP62A and on O-47 strains. In the ordinate axis is reported the percentage of bacterial biofilm production. Data are reported as percentage of residual biofilm after the treatment in comparison with the untreated one. Each data point is composed of 4 independent experiments each performed with at least three replicates.
Molecules 24 00890 g006
Figure 7. Feature importance plot for the 6538P/activation model defined at 3.125 mg/mL.
Figure 7. Feature importance plot for the 6538P/activation model defined at 3.125 mg/mL.
Molecules 24 00890 g007
Figure 8. Feature importance plot for the 25923/activation model defined at 3.125 mg/mL.
Figure 8. Feature importance plot for the 25923/activation model defined at 3.125 mg/mL.
Molecules 24 00890 g008
Figure 9. Feature importance plot for the RP62A/inhibition model defined at 3.125 mg/mL.
Figure 9. Feature importance plot for the RP62A/inhibition model defined at 3.125 mg/mL.
Molecules 24 00890 g009
Figure 10. Feature importance plot for the O-47/inhibition model defined at 3.125 mg/mL.
Figure 10. Feature importance plot for the O-47/inhibition model defined at 3.125 mg/mL.
Molecules 24 00890 g010
Table 1. MIC determined for CGEOs on S. spp strains. EOID indicates samples names. Sample names are the same as previously reported [22]. Ofloxacin MIC is also reported as a positive reference drug. All data are expressed in mg/mL.
Table 1. MIC determined for CGEOs on S. spp strains. EOID indicates samples names. Sample names are the same as previously reported [22]. Ofloxacin MIC is also reported as a positive reference drug. All data are expressed in mg/mL.
EOIDS. aureus 6538PS. aureus 25923S. epidermidis RP62AS. epidermidis O-47
CJ1>25>25>25>25
CJ2>25>2550>25
CJ32525>25>25
CJ62525>25>25
CJ1225>25>25>25
CJ24>25>25>25>25
CJM112.525>25>25
CJM2>25>25>25>25
CJM312.52525–12.525–12.5
CJM412.52525>25
CJM5>25>25>25>25
CA1>25>25>25>25
CA2>25>25>25>25
CA325>25>25>25
CA625–12.5>252525
CA1212.512.512.5>25
CA2412.512.5>25>25
CAM1>25>25>25>25
CAM212.5–6.2512.5–6.2512.512.5
CAM3>25>25>25>25
CAM46.256.2512.512.5
CAM5>25>25>25>25
CS1>25>25>25>25
CS2>25>25>25>25
CS3>25>25>25>25
CS6>25>25>25>25
CS1212.5–6.2512.512.512.5
CS246.2512.5–6.25>2512.5
CSM1>25>25>25>25
CSM212.512.512.512.5
CSM3>25>25>25>25
CSM412.5–6.2512.5–6.2512.512.5
CSM5>25>25>25>25
CO1>25>25>25>25
CO2>25>25>25>25
CO325–12.5252525
CO625–12.525–12.525–12.525–12.5
CO1212.512.512.512.5
CO24>25>25>25>25
COM1252525–12.525–12.5
COM225–12.525–12.525–12.525–12.5
COM3252525–12.525–12.5
COM412.5252525
COM5>25>25>25>25
Ofloxacin0.0002–0.00040.0004–0.00080.0002–0.00040.0002–0.0004
Table 2. MIC determined for FVEOs samples on S. spp strains. EOID indicates samples names. Sample names are the same as previously reported [22]. Ofloxacin MIC is also reported as a positive reference drug. All data are expressed in mg/mL.
Table 2. MIC determined for FVEOs samples on S. spp strains. EOID indicates samples names. Sample names are the same as previously reported [22]. Ofloxacin MIC is also reported as a positive reference drug. All data are expressed in mg/mL.
EOIDS. aureus 6538PS. aureus 25923S. epidermidis RP62AS. epidermidis O-47
FA1>25>25>25>25
FA225–12.525–12.5>25>25
FA3>25>25>25>25
FA6>25>25>25>25
FA12>25>25>25>25
FA2412.5–6.25>25>25>25
FAM1>25>25>25>25
FAM2>25>25>25>25
FAM3>25>25>25>25
FAM4>25>25>25>25
FAM5>25>25>2525
FS1>25>25>25>25
FS2>25>25>25>25
FS3>25>25>25>25
FS6>25>25>25>25
FS12>25>25>25>25
FS24>25>25>25>25
FSM1>25>25>25>25
FSM2>25>25>25>25
FSM312.5–6.25>25>25>25
FSM4>25>25>25>25
FSM525–12.512.5>25>25
FO12525>25>25
FO2>25>25>25>25
FO32525>25>25
FO62525>25>25
FO12>25>25>25>25
FO24>25>25>25>25
FOM1>25>25>25>25
FOM2>25>25>25>25
FOM325–12.512.5>25>25
FOM4>25>25>25>25
FOM5>25>25>25>25
Ofloxacin0.0002–0.00040.0004–0.00080.0002–0.00040.0002–0.0004
Table 3. MIC determined for RSEOs samples on S. spp strains. EOID indicates samples names. Sample names are the same as previously reported [22]. Ofloxacin MIC is also reported as a positive reference drug. All data are expressed in mg/mL.
Table 3. MIC determined for RSEOs samples on S. spp strains. EOID indicates samples names. Sample names are the same as previously reported [22]. Ofloxacin MIC is also reported as a positive reference drug. All data are expressed in mg/mL.
EOIDS. aureus 6538PS. aureus 25923S. epidermidis RP62AS. epidermidis O-47
R1>25>25>25>25
R2>25>25>25>25
R3>25>25>25>25
R6>25>25>25>25
R12>25>25>25>25
R2425>25>25>25
R3025252525
RM1>25>25>25>25
RM2>25>25>25>25
RM3>25>25>25>25
RM4>25>25>25>25
RM5>25>25>25>25
RM6>25>25>25>25
Ofloxacin0.0002–0.00040.0004–0.00080.0002–0.00040.0002–0.0004
Table 4. Data analysis of biofilm production modulation by EOs at the two selected concentrations as reported in Figure 1, Figure 2 and Figure 3. In pale green background are depicted biofilm inhibition values, while in pale red background are highlighted biofilm enhanced data.
Table 4. Data analysis of biofilm production modulation by EOs at the two selected concentrations as reported in Figure 1, Figure 2 and Figure 3. In pale green background are depicted biofilm inhibition values, while in pale red background are highlighted biofilm enhanced data.
Conc. μg/mLS. spp StrainsBiofilm Production %Number EOs Samples at Biofilm Production %
MINMAX<50%<80%<100%≥100%≥120%≥150%≥200%
31256538P50.98523.83010315838229
2592326.92697.4514167347149
RP62A13.04209.6926427118811
O-4727.12289.88243561281444
48.86538P62.80459.4605206937177
2592337.91501.0133467221073
RP62A11.79202.5729487415621
O-470.44306.6031487415742
Table 5. Characteristics of the grid search optimized models.
Table 5. Characteristics of the grid search optimized models.
Assayed Conc. (μg/mL)Models’ ParametersBiofilm Inhibition ModelsBiofilm Activation Models
RP62AO-4725923RP62AO-476538P25923
3125PCs 151922912925
Actives 231304642720
Non-actives 358598583856269
cutoff 626263126133139139
48.8PCs 18924158920
Actives 232303743045
Non-actives 357598682855944
Cutoff4636362124138133121
1: number of principal components used in the model; 2: number of EOs as inhibitors or enhancers of bacterial biofilm production: 3: number of EOs as non-inhibitors or not-enhancers of biofilm production; 4: optimal values of bacterial biofilm production percentage for binary classification as inhibitors/non-inhibitors or enhancers/not-enhancers of bacterial biofilm production.
Table 6. Fitted and cross-validated Accuracy, MCC, Precision-Recall and ROC-AUC coefficients for the RP62A/inhibition, O-47/inhibition, 6538P/activation and 259237/activation optimized models at 3.125 mg/mL and 0.0488 μg/mL.
Table 6. Fitted and cross-validated Accuracy, MCC, Precision-Recall and ROC-AUC coefficients for the RP62A/inhibition, O-47/inhibition, 6538P/activation and 259237/activation optimized models at 3.125 mg/mL and 0.0488 μg/mL.
ValidationAssayed Conc. (μg/mL)CoefficientBiofilm Inhibition ModelsBiofilm Activation Models
RP62AO-476538P25923
Fitting3125Accuracy0.7210.7710.8320.906
MCC0.4550.5900.6670.826
Precision-Recall0.6570.6820.7720.956
ROC-AUC0.7420.7530.8240.961
48.8Accuracy0.7220.7800.8060.763
MCC0.4520.6040.6320.533
Precision-Recall0.6590.6810.7570.824
ROC-AUC0.7350.7520.8150.834
Cross-Validation3125AccuracyCV0.6870.7380.8050.784
MCCCV0.3920.5170.6130.568
Precision-RecallCV0.5840.5890.6980.782
ROC-AUCCV0.6830.6590.7430.845
48.8AccuracyCV0.6630.7210.7220.606
MCCCV0.3350.4740.4500.214
Precision-RecallCV0.5770.5910.6680.533
ROC-AUCCV0.6660.6600.7530.599
Table 7. Chance correlation control by Y-scrambling procedure results. Mean, standard deviation (St Dev), maximum (max) and minimum (min) values for AccuracyY-S, MCCY-S, Precision-RecallY-S and ROC-AUCY-S ROC-AUC coefficients for cross-validated 100 runs. Values refer to RP62A/inhibition, O-47/inhibition, 6538P/activation and 259237/activation optimized models at 3.125 mg/mL.
Table 7. Chance correlation control by Y-scrambling procedure results. Mean, standard deviation (St Dev), maximum (max) and minimum (min) values for AccuracyY-S, MCCY-S, Precision-RecallY-S and ROC-AUCY-S ROC-AUC coefficients for cross-validated 100 runs. Values refer to RP62A/inhibition, O-47/inhibition, 6538P/activation and 259237/activation optimized models at 3.125 mg/mL.
Type of ModelStrainCoefficientMeanSt DevMaxMin
Biofilm Inhibition ModelsRP62AAccuracyY-S0.5000.0790.6440.219
MCCY-S0.0000.1590.290−0.567
Precision-RecallY-S0.4960.0630.6430.353
ROC-AUCY-S0.4590.0940.6270.198
O-47AccuracyY-S0.4940.0810.6650.286
MCCY-S−0.0110.1640.342−0.429
Precision-RecallY-S0.5060.0680.6680.377
ROC-AUCY-S0.4740.0890.6450.241
Biofilm Activation Models6538PAccuracyY-S0.4920.0820.6370.249
MCCY-S−0.0170.1690.275−0.522
Precision-RecallY-S0.5070.0710.6610.347
ROC-AUCY-S0.4700.1000.6440.166
25923AccuracyY-S0.5080.0830.6800.292
MCCY-S0.0130.1700.361−0.433
Precision-RecallY-S0.5240.0710.7460.355
ROC-AUCY-S0.4900.1020.6750.179
Table 8. Feature importances for each chemical component as derived by the SKATER algorithm for the RP62A/inhibition, O-47/inhibition, 6538P/activation and 259237/activation optimized models at 3.125 mg/mL. Background of more important chemical component for inhibition are colored in darker green, while in darker red background are highlighted components associated to higher biofilm production enhancement.
Table 8. Feature importances for each chemical component as derived by the SKATER algorithm for the RP62A/inhibition, O-47/inhibition, 6538P/activation and 259237/activation optimized models at 3.125 mg/mL. Background of more important chemical component for inhibition are colored in darker green, while in darker red background are highlighted components associated to higher biofilm production enhancement.
Chemical ComponentBiofilm Inhibition ModelsBiofilm Inhibition Models
RP62AO-476538P25923
2-Hydroxypiperitenone0.230.320.380.44
2,3-Pinanediol1.334.071.471.05
3-Methylcyclohexanone0.160.451.091.62
3-Octanol6.477.1610.804.44
4-Terpineol1.991.813.653.57
α-Phellandrene0.610.303.020.67
α-Pinene2.712.981.000.21
α-Terpineol2.403.322.812.22
Apiol0.170.780.830.97
β-Cymene1.461.352.426.37
β-Linalool1.401.080.937.08
β-Myrcene0.720.050.950.67
β-Ocimene0.561.750.780.06
β-Phellandrene1.625.423.300.69
β-Pinene1.311.150.342.12
β-Terpinene0.170.630.970.15
Borneol0.070.250.450.56
Carvacrol0.380.380.340.26
Caryophyllene1.351.672.784.36
Caryophyllene oxide0.131.230.323.25
Chrysanthenone3.873.582.684.79
Cinerolone1.140.400.641.45
cis-β-Terpineol0.550.550.881.85
cis-Sabinol4.603.404.360.17
Citral0.120.390.650.99
Cryptone0.151.630.090.27
d-Limonene5.666.299.221.56
delta-Cadinene1.620.900.401.17
Estragole3.873.210.042.49
Fenchone3.663.080.061.54
γ-Terpinene0.191.141.511.97
Germacrene D0.820.650.650.14
Isocaryophyllene0.860.580.840.07
Isomenthone1.530.260.711.77
Isopiperitenone2.802.381.612.42
Isopulegone0.010.634.890.85
Limonene1.644.022.620.67
Menthol2.171.301.387.16
Menthone3.003.073.817.51
Methyl isopulegone0.270.140.270.13
Myristicin4.162.480.160.59
o-Cymene2.646.403.441.62
p-Cymen-8-ol0.901.760.890.96
p-Cymene2.990.020.391.14
p-Menth-1(7)-en-2-one5.060.891.510.91
p-Menthene0.430.330.260.35
Phellandral6.153.971.402.08
Piperitenone0.210.032.864.95
Piperitenone oxide2.320.461.680.36
Pulegone4.064.109.124.18
Sabinene0.610.620.910.62
Terpinolene0.551.901.071.84
Thymol5.732.940.080.31
trans-p-Mentha-2,8-dienol0.430.330.260.35
Table 9. Details of the used bacterial strains.
Table 9. Details of the used bacterial strains.
StrainNameTypeIsolation
S. aureus 6538P6538Pclinical isolateATCC collection
S. aureus 2592325923clinical isolateATCC collection
S. epidermidis RP62ARP62Ainfected catheter isolated strainATCC collection
S. epidermidis O-47O-47septic arthritis clinical isolateHeilmann et al., 1996 [29]

Share and Cite

MDPI and ACS Style

Patsilinakos, A.; Artini, M.; Papa, R.; Sabatino, M.; Božović, M.; Garzoli, S.; Vrenna, G.; Buzzi, R.; Manfredini, S.; Selan, L.; et al. Machine Learning Analyses on Data including Essential Oil Chemical Composition and In Vitro Experimental Antibiofilm Activities against Staphylococcus Species. Molecules 2019, 24, 890. https://doi.org/10.3390/molecules24050890

AMA Style

Patsilinakos A, Artini M, Papa R, Sabatino M, Božović M, Garzoli S, Vrenna G, Buzzi R, Manfredini S, Selan L, et al. Machine Learning Analyses on Data including Essential Oil Chemical Composition and In Vitro Experimental Antibiofilm Activities against Staphylococcus Species. Molecules. 2019; 24(5):890. https://doi.org/10.3390/molecules24050890

Chicago/Turabian Style

Patsilinakos, Alexandros, Marco Artini, Rosanna Papa, Manuela Sabatino, Mijat Božović, Stefania Garzoli, Gianluca Vrenna, Raissa Buzzi, Stefano Manfredini, Laura Selan, and et al. 2019. "Machine Learning Analyses on Data including Essential Oil Chemical Composition and In Vitro Experimental Antibiofilm Activities against Staphylococcus Species" Molecules 24, no. 5: 890. https://doi.org/10.3390/molecules24050890

APA Style

Patsilinakos, A., Artini, M., Papa, R., Sabatino, M., Božović, M., Garzoli, S., Vrenna, G., Buzzi, R., Manfredini, S., Selan, L., & Ragno, R. (2019). Machine Learning Analyses on Data including Essential Oil Chemical Composition and In Vitro Experimental Antibiofilm Activities against Staphylococcus Species. Molecules, 24(5), 890. https://doi.org/10.3390/molecules24050890

Article Metrics

Back to TopTop