Next Article in Journal
Lycocasine A, a Lycopodium Alkaloid from Lycopodiastrum casuarinoides and Its Acid-Sensing Ion Channel 1a Inhibitory Activity
Next Article in Special Issue
Integrating Explicit and Implicit Fullerene Models into UNRES Force Field for Protein Interaction Studies
Previous Article in Journal
Adsorption of Scandium Ions by Amberlite XAD7HP Polymeric Adsorbent Loaded with Tri-n-Octylphosphine Oxide
Previous Article in Special Issue
N-Methyl- and N-Phenylpiperazine Functionalized Styryl Dyes Inside Cucurbiturils: Theoretical Assessment of the Factors Governing the Host–Guest Recognition
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Molecules
Volume 29
Issue 7
10.3390/molecules29071579
molecules-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

DFT Studies of the Activity and Reactivity of Limonene in Comparison with Selected Monoterpenes

by
Katarzyna Rydel-Ciszek
Department of Physical Chemistry, Faculty of Chemistry, Rzeszów University of Technology, Al. Powstańców Warszawy 6, 35-959 Rzeszów, Poland
Molecules 2024, 29(7), 1579; https://doi.org/10.3390/molecules29071579
Submission received: 5 March 2024 / Revised: 27 March 2024 / Accepted: 31 March 2024 / Published: 1 April 2024
(This article belongs to the Special Issue Molecular Modeling: Advancements and Applications II)
Browse Figures
Versions Notes

Abstract

:
Nowadays, the effective processing of natural monoterpenes that constitute renewable biomass found in post-production waste into products that are starting materials for the synthesis of valuable compounds is a way to ensure independence from non-renewable fossil fuels and can contribute to reducing global carbon dioxide emissions. The presented research aims to determine, based on DFT calculations, the activity and reactivity of limonene, an organic substrate used in previous preparative analyses, in comparison to selected monoterpenes such as cymene, pinene, thymol, and menthol. The influence of the solvent model was also checked, and the bonds most susceptible to reaction were determined in the examined compounds. With regard to EHOMO, it was found that limonene reacts more easily than cymene or menthol but with more difficultly than thymol and pienene. The analysis of the global chemical reactivity descriptors “locates” the reactivity of limonene in the middle of the studied monoterpenes. It was observed that, among the tested compounds, the most reactive compound is thymol, while the least reactive is menthol. The demonstrated results can be a reference point for experimental work carried out using the discussed compounds, to focus research on those with the highest reactivity.

Graphical Abstract

1. Introduction

Limonene (1-isopropenyl-4-methylcyclohexene, C10H16) consists of two isoprene units and has two double bonds, one connecting carbon atoms in the ring, an internal or endocyclic bond, and the other occurring outside the ring, that is, external/exocyclic ones [1,2]. Limonene is a monoterpene that occurs naturally in citrus plants such as lemon, orange, and grapes, or even in olive oil [3] and has a wide range of applications. For example, limonene contained in orange peels can be used to produce biogas (using the leaching process) [4], or in the production of biodegradable materials (limonene then serves as an additive) [5]. It was shown that material enriched with a limonene–cyclodextrin/limonene inclusion complex added to poly(L-lactic acid) compared to poly(L-lactic acid) increases water permeability and absorption. This material has an increased barrier to UV–Vis light, which helps protect food against oxidation and also has antibacterial and antifungal properties [5]. The issue of cyclodextrin/limonene complexes is very developing. The incorporation of pure limonene into polymeric materials is difficult. The difficulties are caused by the temperature difference, since limonene has a low boiling point, while polymer processing reactions take place at much higher temperatures. However, the thermal stability of limonene can be increased by trapping it with the formation of an inclusion complex. The use of such a modification allowed the development of a method for obtaining linear polyethylene. Films containing polyethylene and the limonene–cyclodextrin inclusion complex have antifungal and, of course, antibacterial properties, so they can be used in the production of food storage packaging [6]. Interestingly, cyclodextrin and limonene complexes can also be used to improve the durability and aroma of soft drinks [7]. However, there are cases of biomass management in which the antibacterial properties of limonene are a significant disturbing factor; this applies to reactions in which bacteria are used. In the process of obtaining components for the production of biopolymers, using Escherichia coli, which used sugars from orange peel to obtain polyhydroxyalkanoates, limonene, due to its antibacterial properties, had to be eliminated using superoxides produced from potassium peroxide [8]. Rhodococcus Globrulus bacteria are also able to use components of eucalyptus oil, such as limonene, cymene, and thymol, as a source of carbon and energy. The identified and characterized bacterial cytochrome P450 type CYP108N12 is responsible for the biodegradation of these monoterpenes [9].
Limonene, cymene occur in light fractions of pyrolysis oil in car tires [10]. In order to separate these two monoterpenes, the tire pyrolysis oil epoxidation method before the separation was proposed. As a result of the reaction with hydrogen peroxide and a peroxophosphotungstate compound ([(C18H37)2N(CH3)2]3PW4O20), limonene undergoes epoxidation, while cymene does not react under these conditions. The resulting limonene epoxides can be successfully separated using a simple method—distillation, which also expands the possibilities of using waste materials [10].
In 2003, in accordance with the European directive (2003/15/EC), limonene was classified as an allergen [11]. Pure limonene has no allergic properties; however, in the presence of molecular oxygen in the form of long-term contact with air, its autoxidation is possible, and the limonene hydroperoxides formed as a result of this reaction are responsible for allergic reactions [12]. Hydroperoxides acting as an oxidant can oxidize the functional groups in proteins, for example, the sulfur residues of methionine and cysteine or the phenol group in tyrosine [13]. This is one of the reasons why essential oils, including, among others, limonene, cymene, and pinene, are unstable. It is known that limonene undergoes oxidation and degradation in an acidic environment. In order to limit this process, various stabilizers such as whey protein and an electrostatic whey protein–carboxymethyl cellulose complex were used [14]. Nevertheless, a change in the research environment and tests of toilet waters containing limonene showed its high stability [12]. During 9 months, there was no decrease in the initial concentration of limonene, suggesting that limonene is stable in typical water–alcoholic solutions. Unfortunately, a concentrated solution of limonene is much more easily oxidized and, under the influence of air, it may undergo autoxidation to form hydroperoxides [12]. On the other hand, limonene can also be used to deactivate free radicals [15].
Limonene has found practical application in biorefineries, where post-production waste generated during the production of, e.g., orange juices, is processed into commercially useful compounds and constitutes a suitable raw material for the production of important products used in the flavor and fragrance industry [16] in the textile industry for the production of fibers [17] or in pharmacy/medicine [18,19]. However, limonene oxidation products are much more valuable than the substrate from which they are obtained and play an important role as an ingredient for the synthesis of fragrances or drugs [20,21], in the production of biodegradable polymers [22,23,24,25,26], or as a solvent/reactive diluent in the production of epoxy resins [26]. The main products of the C10H16 oxidation reactions are demonstrated in Scheme 1.
In the oxidation of limonene with molecular oxygen, it can undergo an autoxidation reaction, with the formation of limonene hydroperoxide as an intermediate product of the reaction [27]. Hydroxy- and alkylperoxide radicals can be abstracted to allylic hydrogens forming the ketone and alcohol as products, and the acylperoxyl radicals react with the double bond present in the alkene molecules, leading to the epoxide [28]. Bussi et al. [29] showed that the catalyst plays a crucial role in the initial stage of the reaction, consisting of the activation of the reactants and the decomposition of limonene peroxide with the formation of radicals. Nickel and aluminum hydrotalcites were used as catalysts for limonene oxidation by O2 conducted without additional solvent, resulting in the formation of epoxide, alcohol, and ketone [29]. Another heterogeneous complex tested in the oxidation of limonene with dioxygen carried out under mild conditions was the molybdenum (VI) catalyst MoCl2O2Bipy/TiO2-NT with 2,2′-bipyridine-4,4′-dicarboxylate (Bipy) ligand bounded to nanotubes (NT). In the photooxidation reaction (λ = 360 nm) catalyzed by the dioxo-type molybdenum complex, which proceeded practically without the participation of free radicals, limonene-1,2-oxide (LO) was the main product, but carvone (CVN) and carveol (CVL) were also observed [30]. Similar products were observed in another reaction with dioxygen–photooxidation (λ = 360 nm) using TiO2-NT dioxo-Mo (VI) complexes with ligands (L) such as Schiff base, bipyridine, terpyridine. The activity of MoVIO2(L)/TiO2-NT depending on the ligand tested increases in order: Schiff base < bipyridine < terpyridine, where rich in electron ligands, act as “a bridge” for the electron transfer reaction [31]. In turn, by introducing another ligand of the type 2-aminothiazole-4-carboxylic acid and examining the photooxidation reaction of various monoterpenes, pinene was shown to be characterized by greater reactivity than limonene [32]. Using the complexes of iron(II)/(III) [33] and manganese(II) with 2,2′-bipyridine (bpy) [34] formed in situ in the presence of O2, in addition to epoxide, ketone, or alcohol, perillaldehyde and perillyl alcohole were also obtained. Limonene oxide can be selectively produced by the oxidation of limonene with dioxygen and bimetallic complex of ZnCo2O4, the reaction requires isobutyraldehyde as a mediator [35], as well as the use of silylated TiO2 P25 and solar radiation [36].
Oxidants, such as hydrogen peroxide or t-butyl hydroperoxide, are also applied in the oxidation reaction of limonene, and their use was tested in the case of using limonene both as a solvent—the reaction medium and as its substrate, similarly to the use of dioxygen [29]. To ensure hydrophobicity between the limonene as organic solvent and the aqueous phase associated with the addition of an oxidant, the complexes [MoO2(SA(T)P)]2, [MoO2(SATP)]2 [SA(T)P–salicylideneamino(thio)phenolate] for t-Bu-OOH [37] or Ti-salicyldimine with octadecyltrimethoxysilane for HOOH [38] were examined. The data summarized in Table 1 also provide information regarding the use of other oxidants in the limonene oxidation process. Compared to products obtained using dioxygen as an oxidant, diepoxide (DLO), 8,9-LO, or polymer were additionally formed. For HOOH, the following complexes were used: cobalt sandwich-type polyoxometalates [39], tungstophosphates [2], polyoxotungstates [11], Schiff base complexes with Co(II) and Cu(II) and the same compounds but immobilized in zeolite-Y [40], manganese(II) acetylacetonate on MCM41 [41], Al2O3 [42], the ions of non-transition metal [1], methyltrioxorhenium with different ligands [43], activated carbon where the active phase was the magnetite Fe3O4 [44] or MoO2 [45], complexes of VO and copper(II) with Schiff base ligands entrapped in the supercages of zeolite-Y [46], homogeneous and heterogeneous VO and iron(II) with Schiff base ligands [47], γ-Fe2O3/SiO2-NHFeP prepared from nanospheres and 5,10,15,20-tetrakis (pentafluorphenylporphyrin) iron(III) [48], heterogeneous Mn(III), Fe(III), and Co(III) porphyrin-based complexes immobilized on zeolite [49] or others complexes based on zeolite-Y [50,51,52] in which enclosing the catalyst in the porous structure of the support prevents the dimerization of the complexes, ensuring their catalytic activity. Catalysts used with t-butyl hydroperoxide as the oxidant are also zeolites, e.g., zeolite-Y with entrapped VO with Schiff base ligands [50], organic hybrid materials [26,53], Ti-MCM-41, and Ti-MWW compounds [54], iron(II) [55], molybdenum(II) complexes [56], salen-like Jacobsen’s catalysts with manganese(III) [57] or carbon-based complexes with cobalt(II) acetylacetonate [58]. Jacobsen’s compounds with manganese(III) [59,60] with Mn(II), Ni(II), Co(II) [61], or VO(Salten) anchored on SBA-15 (Salten–3-[N,N′-bis-3(salicylidenamino)ethyltriamine]) [62] were also used with other oxidants such as KHSO5 (used as ozone), iodosylbenzene, sodium hypochlorite, or urea hydroperoxide.
Interestingly, some of the mentioned catalysts are also used in the oxidation of other monoterpenes. For example, using analogous catalysts, oxidants, and solvents as given in Table 1, the pinene oxidation reaction was carried out for the conditions presented in Table 1, entries: 1 [30], 2 [32], 16 [40], 23 [47], 25 [52], 26 [51], and the products of these reactions were mainly ketone (verbenone) and alcohol (verbenol). Similarly, the review monograph on biomass management focusses, among others, on the use of both limonene and pinene [64]. In the literature review presented above, in addition to limonene, pinene [14], cymene [9,10], and thymol [9] were mentioned. These natural compounds, similar to limonene, are obtained from plants [65,66,67] and are monoterpenes known for their unique aromatic, therapeutic properties [65,67,68,69]. In turn, due to structural similarity, thymol is often discussed with menthol [70,71,72], a monoterpene resembling hydrogenated thymol, found, e.g., in mint [73], which also has numerous applications [74,75,76]. Interestingly, menthol can be formed from limonene through enzymatic reactions during the monoterpene biosynthetic pathway in peppermint [77]. Therefore, it seems interesting to conduct research on a larger group of monoterpenes, trying to correlate their theoretical activity with data obtained from experiments. In contrast to homogeneous catalysts, their heterogeneous counterparts are recovered, which makes them much more gentle on the environment [47]. Some heterogeneous catalysts are rinsed with large amounts of water and acetone after each catalytic cycle for reuse [46]. Rinsing with water is also carried out to remove undesirable ions, like chloride, and this reaction may contribute to coordinated H2O molecules in space, e.g., zeolite. Furthermore, it has been shown that such complexes (with coordinated water molecules) provide a higher conversion of limonene [46]. On the other hand, the presence of water affects the products of the limonene oxidation reaction and contributes to the ring opening reactions, as a result of which a diol (LDIOL) is obtained from epoxide (LO) [37]. Thus, it seems advisable to investigate the influence of solvents, including water, on the activity of monoterpenes. This can be performed using computational chemistry methods assuming appropriate PMC solvent models. Additionally, although the thesis, that all reactions can be carried out by selecting appropriate catalysts and “additives” as co-catalysts, presented in the review is true [78]. However, in terms of planning preparatory research and selecting appropriate substrates for various reactions, it is crucial to determine their activity in order to exclude the least reactive compounds in the initial stage of the research, which was used, for example, to separate a mixture of oils [10].
For this reason, using DFT computational methods, monoterpenes were examined to determine their activity and to check the possibility of a potential attack of selected monoterpenes on the empty orbital of the catalyst, which is consistent with the Dewar–Chatt–Duncanson model [30]. Monoterpenes were selected based on the review of the literature presented in the manuscript as compounds that can occur in oils in the presence of limonene, the substrate that is the object of my interest and previous research [79]. Determining the differences in the reactivity of analyzed compounds has many applications; for example, it may contribute to the separation of oil components from their mixture, as was achieved in [10], in which the limonene oxidation product (LO) was successfully separated from cymene by distillation. The influence of the use of the solvent model on selected monoterpenes was also examined, and their bonds that are most susceptible to breaking were determined.

2. Results and Discussion

Limonene and the structures of selected monoterpenes such as cymene, pinene, thymol, and menthol were optimized using methods with different hybrid functionals and basis sets, Table S1. However, for the tested compounds, the best correlation with experimental data [10] was provided by the B3LYP/6-311+G(d) level method; therefore, it was used to calculate reactivity descriptors. These structures, along with the numbering of carbon atoms presented in GaussView03 and used to discuss BDE, are shown in Scheme 2.
For the optimal structures of selected monoterpenes generated after the calculations performed in the gaseous phase, calculations were also conducted, assuming that the PCM solvent model works for solvents with different polarities, such as water, acetonitrile, and methanol. For each optimized structure, the energy of the highest occupied molecular orbital (HOMO) and the lowest unoccupied molecular orbital (LUMO) was generated, and these data along with the Egap value are presented in Table 2. Those frontier orbitals are the most important orbitals in terms of reactivity [80]. The HOMO energy determines the molecule’s susceptibility to electrophilic attack and is related to the ionization potential. Analogously, the LUMO energy determines the molecule’s predisposition to nucleophilic attack and is related to electron affinity [81]. In turn, the lower the Egap, the less energy is needed to transfer an electron from the HOMO orbital to the LUMO [80]. Therefore, knowing the difference between the HOMO–LUMO energy, it is also possible to determine which of the tested molecules is characterized by the greatest kinetic stability [82].
Analyzing the Egap values collected in Table 2, it was found that, regardless of the solvent model used, the most reactive monoterpene is thymol, which is characterized by the lowest stability, while menthol will be the least reactive of the group of compounds tested. It was also observed that, depending on the tested monoterpene, Egap (Egap = ELUMO − EHOMO [83,84]) can be used to determine the dependency on the solvent model used. For example, on the basis of the data collected in Figure 1, which visualize the surfaces of the HOMO, the LUMO molecular orbitals with their corresponding energies for different ε for limonene, it can be observed that the HOMO–LUMO energy difference increases with increasing ε and is the largest for water, then acetonitrile, and then methanol.
However, analyzing the data collected in Table 2 for the other monoterpene molecules requires individual analysis. Pinene, like limonene, has the largest Egap for using water as a solvent model and the smallest without assuming a solvent model, while the difference between the ELUMO and EHOMO orbitals for cymene for water as a solvent model is the smallest. However, when using the solvent model, remember that it is a polarized continuum model in which solvents are represented by the dielectric continuum medium, and therefore the analyzed models should be verified by conducting experimental work.
Nevertheless, of the ε value, the HOMO energy values of the tested monoterpenes can be arranged as follows: thymol > pienene > limonene > cymene > menthol, as shown in Figure 2, thymol is the compound with the highest EHOMO (in all PCM models), while menthol, among the selected monoterpenes, is the compound with the lowest value of this energy.
In the reaction of monoterpenes with an oxygen-activated catalyst, in the case of electron-rich olefins, the oxygen atom transfer step is based, according to the Dewar–Chatt–Duncanson model, on the attack of the (nucleophilic) olefin on the empty LUMO orbital of the oxygen-catalyst [30]. Therefore, having the values of the highest occupied molecular orbitals, it can be concluded that the monoterpene with the highest EHOMO will be the most reactive in atom transfer reactions. Based on Figure 2, it follows that thymol and pinene are more reactive than limonene, while cymene and menthol are less reactive. The presented results are consistent with previously presented literature data that showed that pinene was characterized by greater reactivity than limonene [32], while from another article, it follows that cymene compared to limonene does not react [10].
The energy of HOMO and LUMO is used to determine the global chemical reactivity descriptors (GCRD)—Table 3. GCRDs are defined for monoterpene molecules in their singlet ground state with the DFT of Parr, Pearson, and Yang [85]. The ionization potential (I), correlated with –EHOMO, is the minimum energy that is necessary to remove an electron from a monoterpene molecule [84]. In all the cases listed in Table 3, thymol has the lowest ionization potential value, followed by pinene, limonene, and cymene. Menthol has the highest ionization potential. In turn, the electron affinity (A) is determined on the basis of the value of –ELUMO and characterizes the ability to attach an electron, resulting in the formation of a negative ion. Table 3 shows that, for the results of calculations carried out without assuming the solvent model, thymol has the highest A value. The situation is similar when using MeOH as a PCM model—thymol is also characterized by the highest electron affinity value of all the tested compounds. In turn, menthol has the lowest A value in all the tested solvent models. The electronegativity (X) values, i.e., the measure of the tendency to “attract” electrons [84], is defined as X = 0.5∙(I + A) [86,87,88] and increase in series: pinene, limonene, thymol, cymene, and menthol with the highest value. Global hardness (η) according to Parr and Pearson—first-order derivative of the chemical potential with respect to the total number of electrons N, with a constant external potential, or second-order derivative of energy (also with respect to the number of electrons N, at a constant external potential). The global hardness is calculated on the basis of the knowledge of I and A. Global softness (S) is related to η. Global hardness and softness concern the sensitivity of electron–electron interactions; for example, for anions that are characterized by the lowest hardness value and the highest softness, their susceptibility to changing the number of electrons is low [89]. For the analyzed monoterpenes, molecular hardness η = 0.5∙(I – A) and molecular softness S = 0.5/η [86,87,88] were calculated; the η values increase from thymol (with the lowest η value), cymene, pienene, limonene, and menthol (with the highest η value and the lowest S value). The global softness values in the given series decrease. In Table 3, ω (where ω = μ2/(2∙η) [88,90]) is characterized by the electrophilicity index, which describes the global electrophilic nature of molecules. It expresses the measure of energy reduction resulting from the flow of electrons between a donor and an acceptor. In the case of a reaction, a molecule with a higher ω will react as an electrophile, while a molecule with a lower ω—as a nucleophile [88]. For example, using water as the PCM model, menthol has high values of both the ω and the X descriptor, so it can act as an electrophile. The last column of Table 3 applies to μ, which describes the chemical potential; μ = −0.5∙(I + A) [90]. Chemical potential shows the sensitivity of the system to changes in electrons, a high μ indicates that the molecule has the properties of a strong electron acceptor, while a low μ characterizes strong electron donors [84]. Of the compounds analyzed, menthol has the lowest μ value; a low μ and high ω for a molecule indicate its good electrophilic nature.
Based on the data collected in Table 3, it can be concluded that, in light of the theory of hard and soft acids and bases (HSAB), thymol should react the easiest/fastest because of its lowest η value, which proves its nucleophilic properties and proton acceptor capabilities, while menthol is the least reactive of the tested group of compounds.
Furthermore, the dissociation enthalpies (BDEs) of individual C–H in the molecules of the tested monoterpenes were calculated according to the reaction:
BDE = Hmonoterpene without-H + HHatom − Hmonoterpene.
The energies required for homolytic breakage of a specific bond calculated using two calculation methods are presented in Tables S2–S6. The numbering of carbon atoms in monoterpene molecules is consistent with the numbering shown in Scheme 2. The data in Tables S2–S6 are listed from the lowest BDE values, i.e., from the bonds in monoterpene molecules that are easiest to break. When comparing the results of the BDE calculations, it can be concluded that, regardless of the chosen calculation method, they provide consistent results. The lowest energy needed to remove a hydrogen atom from the analyzed molecules (rounded to whole numbers) is 81 kcal/mol for the limonene molecule, 77 kcal/mol for thymol, 74 kcal/mol for pinene, 83 kcal/mol for cymene, and 91 kcal/mol for menthol. Therefore, pinene and thymol can most easily undergo reactions involving the transfer of a hydrogen atom to the catalyst molecule, while menthol from the tested group of compounds is the most difficult to oxidize. Data from Tables S2 and S4 show that, for limonene and pinene, the oxidations should occur more easily in the allylic position, which is consistent with the literature data [47]. Additionally, this experimental work confirms that pinene is more reactive than limonene, and one of the main products resulting from the oxidation of both limonene and pinene under similar conditions are their allylic derivatives (Table 1, entries: 16 [40], 23 [47], 25 [52], 26 [51]). In turn, in the case of cymene and limonene and using the same reaction conditions, these monoterpenes were found to oxidize at carbon C7 [9,91]. Interestingly, in the case of both monoterpenes, the calculated BDE values for these bonds are comparable (Tables S2 and S3). In the case of BDE calculated for cymene (Table S3, Scheme S1), the C8-H bond is the easiest to break, which is also confirmed by the products obtained in preparative experiments [92]. In thymol and menthol molecules, the breaking of bonds leading to the formation of the corresponding ketones (Tables S5 and S6, Scheme S1) is favored [77,93,94,95].

3. Materials and Methods

Calculations of thermodynamic parameters were performed in Gaussian 09 and 16 [96] using DFT methods with the Becke 3-parameter hybrid density functional, Lee Young Parr correlation B3LYP [97] and the basis sets 6-31g(d), 6-311+G, Def2SVP, or functionals CAM-B3LYP; B3PW91; ωB97XD [98,99] with 6-311+G(d). The GaussView03 programme was used to model the structures of monoterpenes molecules. Geometry optimization was performed using the B3LYP/6-31G(d), 6-311+G(d) level or Def2SVP method [100,101] using as PCM model water (ε = 78.3553), acetonitrile (ε = 35.688), and methanol (ε = 32.613). The values reported in this paper, combining electronic energies with the enthalpy correction, were used for bond dissociation energy (BDE) calculations. BDE is expressed as a change in the enthalpy of the homolytic cleavage of a selected bond [102] and is one of the basic features of the reactivity of selected compounds.

4. Conclusions

The presented results can contribute to the rational planning and optimization of the experimental work. On the basis of the calculations of the ionization potential, electron affinity, global hardness and softness, the electrophilicity index, electronegativity, and chemical potential, it was found that of the analyzed monoterpenes, thymol will be more reactive than limonene. The most stable and least reactive is menthol, which means that its presence, among others, in essential oils is least exposed to subsequent reactions. Additionally, the energy values of the highest occupied molecular orbital show that EHOMO limonene is “in the middle” of the tested monoterpenes and allow for the ranking of the examined compounds from the most reactive ones, respectively: thymol, pienene, limonene, cymene, and menthol. These calculations are partially confirmed by experimental data [14,32]. The use of the solvent model does not significantly affect the structures of the analyzed monoterpenes. The calculated values of bond dissociation enthalpy also confirm the conclusions obtained from the GCRD analysis that pinene and thymol undergo hydrogen atom transfer reactions more easily than limonene.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/molecules29071579/s1, Table S1: Energy of the HOMO and LUMO orbitals along with ΔEgap levels for limonene and selected monoterpenes, structures optimized using different methods without the PCM model; Table S2: The energies (with and without zero point correction), enthalpies, free energies (G), and bound dissociation enthalpy (BDE) values for the limonene molecules and its radicals were calculated using B3LYP and water as the PCM model; Table S3: The energies (with and without zero point correction), enthalpies, free energies (G), and bound dissociation enthalpy (BDE) values for the cymene molecules and its radicals were calculated using B3LYP and water as the PCM model; Table S4: The energies (with and without zero point correction), enthalpies, free energies (G), and bound dissociation enthalpy (BDE) values for the pinene molecules and its radicals were calculated using B3LYP and water as the PCM model; Table S5: The energies (with and without zero point correction), enthalpies, free energies (G), and bound dissociation enthalpy (BDE) values for the thymol molecules and its radicals were calculated using B3LYP and water as the PCM model; Table S6: The energies (with and without zero point correction), enthalpies, free energies (G), and bound dissociation enthalpy (BDE) values for the menthol molecules and its radicals were calculated using B3LYP and water as the PCM model; Scheme S1: Possible oxidation products of selected monoterpenes.

Funding

Financed by the Minister of Science and Higher Education Republic of Poland within the program “Regional Excellence Initiative”.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The data presented in this study are available in the Supplementary Materials or on request from the corresponding author.

Acknowledgments

This research was supported in part by PLGrid Infrastructure. Computations were performed in the AGH Cyfronet Supercomputer Centre.

Conflicts of Interest

The author declares no conflicts of interest.

References

  1. Mandelli, D.; Kozlov, Y.N.; da Silva, C.A.R.; Carvalho, W.A.; Pescarmona, P.P.; Cella, D.d.A.; de Paiva, P.T.; Shul’pin, G.B. Oxidation of olefins with H2O2 catalyzed by gallium(III) nitrate and aluminum(III) nitrate in solution. J. Mol. Catal. A Chem. 2016, 422, 216–220. [Google Scholar] [CrossRef]
  2. Santos, I.C.M.S.; Gamelas, J.A.F.; Duarte, T.A.G.; Simões, M.M.Q.; Neves, M.G.P.M.S.; Cavaleiro, J.A.S.; Cavaleiro, A.M.V. Catalytic homogeneous oxidation of monoterpenes and cyclooctene with hydrogen peroxide in the presence of sandwich-type tungstophosphates [M4(H2O)2(PW9O34)2]n−, M = CoII, MnII and FeIII. J. Mol. Catal. A Chem. 2017, 426, 593–599. [Google Scholar] [CrossRef]
  3. Al-Dabbas, M.M.; Al-Jaloudi, R.; Abdullah, M.A.; Abughoush, M. Characterization of Olive Oil Volatile Compounds after Elution through Selected Bleaching Materials-Gas Chromatography-Mass Spectrometry Analysis. Molecules 2023, 28, 6444. [Google Scholar] [CrossRef] [PubMed]
  4. Wikandari, R.; Nguyen, H.; Millati, R.; Niklasson, C.; Taherzadeh, M.J. Improvement of Biogas Production from Orange Peel Waste by Leaching of Limonene. BioMed Res. Int. 2015, 2015, 494182. [Google Scholar] [CrossRef] [PubMed]
  5. Dobrzyńska-Mizera, M.; Knitter, M.; Szymanowska, D.; Mallardo, S.; Santagata, G.; Di Lorenzo, M.L. Optical, mechanical, and antimicrobial properties of bio-based composites of poly(L-lactic acid) and D-limonene/β-cyclodextrin inclusion complex. J. App. Polym. Sci. 2022, 139, 52177. [Google Scholar] [CrossRef]
  6. Dobrzyńska-Mizera, M.; Knitter, M.; Piss, M.; Szymanowska, D.; Mallardo, S.; Santagata, G.; Di Lorenzo, M.L. Linear low-density polyethylene modified with d-limonene/β-cyclodextrin inclusion complex: Antimicrobial composite for active food packaging. Polym. Eng. Sci. 2024, 64, 52–61. [Google Scholar] [CrossRef]
  7. Saldanha do Carmo, C.; Pais, R.; Simplício, A.L.; Mateus, M.; Duarte, C.M.M. Improvement of Aroma and Shelf-Life of Non-alcoholic Beverages Through Cyclodextrins-Limonene Inclusion Complexes. Food Bioprocess Technol. 2017, 10, 1297–1309. [Google Scholar] [CrossRef]
  8. Davaritouchaee, M.; Mosleh, I.; Dadmohammadi, Y.; Abbaspourrad, A. One-Step Oxidation of Orange Peel Waste to Carbon Feedstock for Bacterial Production of Polyhydroxybutyrate. Polymers 2023, 15, 697. [Google Scholar] [CrossRef]
  9. Giang, P.D.; Churchman, L.R.; Stok, J.E.; Bell, S.G.; De Voss, J.J. Cymredoxin, a [2Fe-2S] ferredoxin, supports catalytic activity of the p-cymene oxidising P450 enzyme CYP108N12. Arch. Biochem. Biophys. 2023, 737, 109549. [Google Scholar] [CrossRef]
  10. Zhang, Y.; Cheng, S.; Ji, G.; Li, A. Refinement of limonene epoxides from the light distillate of tire pyrolysis oil via catalytic epoxidation. Sep. Purif. Technol. 2023, 319, 124068. [Google Scholar] [CrossRef]
  11. Makrygenni, O.; Vanmairis, L.; Taourit, S.; Launay, F.; Shum Cheong Sing, A.; Proust, A.; Gérard, H.; Villanneau, R. Selective Formation of Epoxylimonene Catalyzed by Phosphonyl/Arsonyl Derivatives of Trivacant Polyoxotungstates at Low Temperature. Eur. J. Inorg. Chem. 2020, 2020, 605–612. [Google Scholar] [CrossRef]
  12. Kern, S.; Granier, T.; Dkhil, H.; Haupt, T.; Ellis, G.; Natsch, A. Stability of limonene and monitoring of a hydroperoxide in fragranced products. Flavour Fragr. J. 2014, 29, 277–286. [Google Scholar] [CrossRef]
  13. Christensson, J.B.; Johansson, S.; Hagvall, L.; Jonsson, C.; Börje, A.; Karlberg, A.T. Limonene hydroperoxide analogues differ in allergenic activity. Contact Dermat. 2008, 59, 344–352. [Google Scholar] [CrossRef] [PubMed]
  14. Kaade, W.; Méndez-Sánchez, C.; Güell, C.; De Lamo-Castellví, S.; Mestres, M.; Ferrando, M. Complexed Biopolymer of Whey Protein and Carboxymethyl Cellulose to Enhance the Chemical Stability of Lemon Oil-in-Water Emulsions. ACS Food Sci. Technol. 2022, 2, 41–48. [Google Scholar] [CrossRef]
  15. Francisco-Márquez, M.; Galano, A. Limonene: A scented and versatile tropospheric free radical deactivator. Int. J. Quantum Chem. 2023, 123, e27103. [Google Scholar] [CrossRef]
  16. Ravichandran, C.; Badgujar, P.C.; Gundev, P.; Upadhyay, A. Review of toxicological assessment of d-limonene, a food and cosmetics additive. Food Chem. Toxicol. Int. J. Publ. Br. Ind. Biol. Res. Assoc. 2018, 120, 668–680. [Google Scholar] [CrossRef] [PubMed]
  17. Lan, W.; Liang, X. Electrospun Polyvinyl Alcohol/d-Limonene Fibers Prepared by Ultrasonic Processing for Antibacterial Active Packaging Material. Molecules 2019, 24, 767. [Google Scholar] [CrossRef] [PubMed]
  18. Shukla, P.; Pant, A.; Pandey, R.K. Limonene Attenuates Oxidative Stress and Extends Longevity in Caenorhabditis elegans. Curr. Sci. 2019, 116, 959–965. [Google Scholar] [CrossRef]
  19. Sales, A.; Pastore, G.M.; Bicas, J.L. Optimization of limonene biotransformation to limonene-1,2-diol by Colletotrichum nymphaeae CBMAI 0864. Process Biochem. 2019, 86, 25–31. [Google Scholar] [CrossRef]
  20. De Oliveira, M.P.; Delolo, F.G.; Villarreal, J.A.A.; dos Santos, E.N.; Gusevskaya, E.V. Hydroformylation and one-pot hydroformylation/epoxy ring cleavage of limonene oxide: A sustainable access to biomass-based multi-functional fragrances. Appl. Catal. A Gen. 2021, 616, 118082. [Google Scholar] [CrossRef]
  21. Aissou, M.; Chemat-Djenni, Z.; Yara-Varón, E.; Fabiano-Tixier, A.-S.; Chemat, F. Limonene as an agro-chemical building block for the synthesis and extraction of bioactive compounds. Comptes Rendus Chim. 2017, 20, 346–358. [Google Scholar] [CrossRef]
  22. Claudino, M.; Mathevet, J.-M.; Jonsson, M.; Johansson, M. Bringing d-limonene to the scene of bio-based thermoset coatings via free-radical thiol–ene chemistry: Macromonomer synthesis, UV-curing and thermo-mechanical characterization. Polym. Chem. 2014, 5, 3245–3260. [Google Scholar] [CrossRef]
  23. Hauenstein, O.; Reiter, M.; Agarwal, S.; Rieger, B.; Greiner, A. Bio-based polycarbonate from limonene oxide and CO2 with high molecular weight, excellent thermal resistance, hardness and transparency. Green Chem. 2016, 18, 760–770. [Google Scholar] [CrossRef]
  24. Li, C.; Sablong, R.J.; van Benthem, R.A.T.M.; Koning, C.E. Unique Base-Initiated Depolymerization of Limonene-Derived Polycarbonates. ACS Macro Lett. 2017, 6, 684–688. [Google Scholar] [CrossRef] [PubMed]
  25. Brito, A.C.; Silvino, C.A. A Brief Review of the Use of Limonene Derivative as a Promising Monomer in the Synthesis of Biodegradable Polymers. Rev. Virtual Química 2021, 13, 1017–1041. [Google Scholar] [CrossRef]
  26. Nunes, M.S.; Gomes, D.M.; Gomes, A.C.; Neves, P.; Mendes, R.F.; Paz, F.A.A.; Lopes, A.D.; Valente, A.A.; Gonçalves, I.S.; Pillinger, M. A 5-(2-Pyridyl)tetrazolate Complex of Molybdenum(VI), Its Structure, and Transformation to a Molybdenum Oxide-Based Hybrid Heterogeneous Catalyst for the Epoxidation of Olefins. Catalysts 2021, 11, 1407. [Google Scholar] [CrossRef]
  27. Kuznetsova, G.M.; Lobanova, T.V.; Rusina, I.F.; Kasaikina, O.T. Kinetic characteristics of initiated oxidation of limonene. Russ. Chem. Bull. 1996, 45, 1586–1591. [Google Scholar] [CrossRef]
  28. Wentzel, B.B.; Alsters, P.L.; Feiters, M.C.; Nolte, R.J.M. Mechanistic Studies on the Mukaiyama Epoxidation. J. Org. Chem. 2004, 69, 3453–3464. [Google Scholar] [CrossRef]
  29. Pena, A.; Veiga, S.; Sapelli, M.; Martínez, N.; Márquez, V.; Dellacassa, E.; Bussi, J. Limonene oxidation by molecular oxygen under solvent-free conditions: The influence of peroxides and catalysts on the reaction rate. React. Kinet. Mech. Catal. 2012, 107, 263–275. [Google Scholar] [CrossRef]
  30. Martinez, Q.H.; Amaya, Á.A.; Paez-Mozo, E.A.; Martinez, O.F.; Valange, S. Photo-assisted O-atom transfer to monoterpenes with molecular oxygen and a dioxoMo(VI) complex immobilized on TiO2 nanotubes. Catal. Today 2021, 375, 441–457. [Google Scholar] [CrossRef]
  31. Martínez, Q.H.; Paez-Mozo, E.A.; Martínez, O.F. Selective Photo-epoxidation of (R)-(+)- and (S)-(−)-Limonene by Chiral and Non-Chiral Dioxo-Mo(VI) Complexes Anchored on TiO2-Nanotubes. Top. Catal. 2021, 64, 36–50. [Google Scholar] [CrossRef]
  32. Martinez Quiñonez, H.; Amaya, Á.A.; Paez-Mozo, E.A.; Martinez Ortega, F. Aminothiazole Ligand-Type Dioxo-Mo(VI) Complex Anchored on TiO2 Nanotubes for Selective Oxidation of Monoterpenes with Light and O2. Top. Catal. 2022, 65, 1088–1101. [Google Scholar] [CrossRef]
  33. Naróg, D.; Szczepanik, A.; Sobkowiak, A. Iron(II, III)-Catalyzed Oxidation of Limonene by Dioxygen. Catal. Lett. 2008, 120, 320–325. [Google Scholar] [CrossRef]
  34. Szczepanik, A.; Sobkowiak, A. Manganese(II)-Induced Oxidation of Limonene by Dioxygen. Catal. Lett. 2008, 126, 261–267. [Google Scholar] [CrossRef]
  35. Liu, J.; Ji, X.; Wang, C.; Wang, L.; Jian, P. Beneficial Synergistic Intermetallic Effect in ZnCo2O4 for Enhancing the Limonene Oxidation Catalysis. Inorg. Chem. 2023, 62, 18750–18757. [Google Scholar] [CrossRef]
  36. Ciriminna, R.; Parrino, F.; De Pasquale, C.; Palmisano, L.; Pagliaro, M. Photocatalytic partial oxidation of limonene to 1,2 limonene oxide. Chem. Commun. 2018, 54, 1008–1011. [Google Scholar] [CrossRef]
  37. Wang, W.; Agustin, D.; Poli, R. Influence of ligand substitution on molybdenum catalysts with tridentate Schiff base ligands for the organic solvent-free oxidation of limonene using aqueous TBHP as oxidant. Mol. Catal. 2017, 443, 52–59. [Google Scholar] [CrossRef]
  38. Yuan, L.S.; Chandren, S.; Efendi, J.O.N.; Ho, C.S.; Nur, H. Hydrophobic effect of silica functionalized with silylated Ti-salicylaldimine complex on limonene oxidation by aqueous hydrogen peroxide. J. Chem. Sci. 2015, 127, 1905–1917. [Google Scholar] [CrossRef]
  39. Kala Raj, N.K.; Puranik, V.G.; Gopinathan, C.; Ramaswamy, A.V. Selective oxidation of limonene over sodium salt of cobalt containing sandwich-type polyoxotungstate [WCo3(H2O)2{W9CoO34}2]10−. Appl. Catal. A Gen. 2003, 256, 265–273. [Google Scholar] [CrossRef]
  40. Godhani, D.R.; Nakum, H.D.; Parmar, D.K.; Mehta, J.P.; Desai, N.C. Zeolite-Y immobilized Metallo-ligand complexes: A novel heterogenous catalysts for selective oxidation. Inorg. Chem. Commun. 2016, 72, 105–116. [Google Scholar] [CrossRef]
  41. Oliveira, P.; Machado, A.; Ramos, A.M.; Fonseca, I.M.; Braz Fernandes, F.M.; Botelho do Rego, A.M.; Vital, J. Anchoring manganese acetylacetonate complex on MCM-41: Catalytic testing on limonene oxidation. Catal. Commun. 2007, 8, 1366–1372. [Google Scholar] [CrossRef]
  42. Bonon, A.J.; Kozlov, Y.N.; Bahú, J.O.; Filho, R.M.; Mandelli, D.; Shul’pin, G.B. Limonene epoxidation with H2O2 promoted by Al2O3: Kinetic study, experimental design. J. Catal. 2014, 319, 71–86. [Google Scholar] [CrossRef]
  43. Michel, T.; Cokoja, M.; Sieber, V.; Kühn, F.E. Selective epoxidation of (+)-limonene employing methyltrioxorhenium as catalyst. J. Mol. Catal. A Chem. 2012, 358, 159–165. [Google Scholar] [CrossRef]
  44. Młodzik, J.; Wróblewska, A.; Makuch, E.; Wróbel, R.J.; Michalkiewicz, B. Fe/EuroPh catalysts for limonene oxidation to 1,2-epoxylimonene, its diol, carveol, carvone and perillyl alcohol. Catal. Today 2016, 268, 111–120. [Google Scholar] [CrossRef]
  45. Wróblewska, A.; Makuch, E.; Młodzik, J.; Koren, Z.C.; Michalkiewicz, B. Oxidation of limonene over molybdenum dioxide-containing nanoporous carbon catalysts as a simple effective method for the utilization of waste orange peels. React. Kinet. Mech. Catal. 2018, 125, 843–858. [Google Scholar] [CrossRef]
  46. Modi, C.K.; Gade, B.G.; Chudasama, J.A.; Parmar, D.K.; Nakum, H.D.; Patel, A.L. Synthesis, spectral investigation and catalytic aspects of entrapped VO(IV) and Cu(II) complexes into the supercages of zeolite-Y. Spectrochim. Acta Part A Mol. Biomol. Spectrosc. 2015, 140, 174–184. [Google Scholar] [CrossRef] [PubMed]
  47. Parmar, D.K.; Butani, P.M.; Thumar, N.J.; Jasani, P.M.; Padaliya, R.V.; Sandhiya, P.R.; Nakum, H.D.; Khan, M.N.; Makwana, D. Oxy-functionalization of olefins with neat and heterogenized binuclear V(IV)O and Fe(II) complexes: Effect of steric hindrance on product selectivity and output in homogeneous and heterogeneous phase. Mol. Catal. 2019, 474, 110424. [Google Scholar] [CrossRef]
  48. Santos, J.S.D.; Faria, A.L.; Amorin, P.M.; Luna, F.M.L.; Caiado, K.L.; Silva, D.O.; Sartoratto, P.P.; Assis, M.D. Iron(III) Porphyrin Covalently Supported onto Magnetic Amino-Functionalized Nanospheres as Catalyst for Hydrocarbon and Herbicide Oxidations. J. Braz. Chem. Soc. 2012, 23, 1411–1420. [Google Scholar] [CrossRef]
  49. Madadi, M.; Rahimi, R. Zeolite-immobilized Mn(III), Fe(III) and Co(III) complexes with 5,10,15,20-tetra(4-methoxyphenyl)porphyrin as heterogeneous catalysts for the epoxidation of (R)-(+)-limonene: Synthesis, characterization and catalytic activity. React. Kinet. Mech. Catal. 2012, 107, 215–229. [Google Scholar] [CrossRef]
  50. Modi, C.K.; Chudasama, J.A.; Nakum, H.D.; Parmar, D.K.; Patel, A.L. Catalytic oxidation of limonene over zeolite-Y entrapped oxovanadium (IV) complexes as heterogeneous catalysts. J. Mol. Catal. A Chem. 2014, 395, 151–161. [Google Scholar] [CrossRef]
  51. Godhani, D.R.; Nakum, H.D.; Parmar, D.K.; Mehta, J.P.; Desai, N.C. Zeolite Y encaged Ru(III) and Fe(III) complexes for oxidation of styrene, cyclohexene, limonene, and α-pinene: An eye-catching impact of H2SO4 on product selectivity. J. Mol. Catal. A Chem. 2017, 426, 223–237. [Google Scholar] [CrossRef]
  52. Mehta, J.P.; Parmar, D.K.; Nakum, H.D.; Godhani, D.R.; Desai, N.C. Enhanced catalytic oxidation of monoterpenes by zeolite-Y entrapped iron complex: Spectral studies and mechanistic vision. J. Porous Mater. 2018, 25, 1649–1658. [Google Scholar] [CrossRef]
  53. Nunes, M.S.; Gomes, D.M.; Gomes, A.C.; Neves, P.; Mendes, R.F.; Paz, F.A.A.; Lopes, A.D.; Pillinger, M.; Valente, A.A.; Gonçalves, I.S. A Molybdenum(VI) Complex of 5-(2-pyridyl-1-oxide)tetrazole: Synthesis, Structure, and Transformation into a MoO3-Based Hybrid Catalyst for the Epoxidation of Bio-Olefins. Catalysts 2023, 13, 565. [Google Scholar]
  54. Gawarecka, A.; Wróblewska, A. Limonene oxidation over Ti-MCM-41 and Ti-MWW catalysts with t-butyl hydroperoxide as the oxidant. React. Kinet. Mech. Catal. 2018, 124, 523–543. [Google Scholar] [CrossRef]
  55. Becerra, J.-A.; González, L.-M.; Villa, A.-L. A bio-inspired heterogeneous catalyst for the transformation of limonene from orange peel waste biomass into value-added products. Catal.Today 2018, 302, 250–260. [Google Scholar] [CrossRef]
  56. Abrantes, M.; Bruno, S.M.; Tomé, C.; Pillinger, M.; Gonçalves, I.S.; Valente, A.A. Epoxidation of DL-limonene using an indenyl molybdenum(II) tricarbonyl complex as catalyst precursor. Catal. Commun. 2011, 15, 64–67. [Google Scholar] [CrossRef]
  57. Oliveira, P.; Machado, A.; Ramos, A.M.; Fonseca, I.; Fernandes, F.M.B.; Rego, A.M.B.d.; Vital, J. MCM-41 anchored manganese salen complexes as catalysts for limonene oxidation. Microporous Mesoporous Mater. 2009, 120, 432–440. [Google Scholar] [CrossRef]
  58. Oliveira, P.; Ramos, A.M.; Fonseca, I.; Botelho do Rego, A.; Vital, J. Oxidation of limonene over carbon anchored transition metal Schiff base complexes: Effect of the linking agent. Catal. Today 2005, 102–103, 67–77. [Google Scholar]
  59. De Fátima Teixeira Gomes, M.; Antunes, O.A.C. Oxidation of limonene catalyzed by MnIII(Salen)Cl-H2O. Catal. Lett. 1996, 38, 133–134. [Google Scholar] [CrossRef]
  60. Cubillos, J.; Vásquez, S.; Montes de Correa, C. Salen manganese (III) complexes as catalysts for R-(+)-limonene oxidation. Appl. Catal. A Gen. 2010, 373, 57–65. [Google Scholar] [CrossRef]
  61. Lima, L.F.; Corraza, M.L.; Cardozo-Filho, L.; Márquez-Alvarez, H.; Antunes, O.A.C. Oxidation of limonene catalyzed by Metal(Salen) complexes. Braz. J. Chem. Eng. 2006, 23, 83–92. [Google Scholar] [CrossRef]
  62. Joseph, T.; Halligudi, S.B. Oxyfunctionalization of limonene using vanadium complex anchored on functionalized SBA-15. J. Mol. Catal. A Chem. 2005, 229, 241–247. [Google Scholar] [CrossRef]
  63. Gonçalves, J.A.; Bueno, A.C.; Gusevskaya, E.V. Palladium-catalyzed oxidation of monoterpenes: Highly selective syntheses of allylic ethers from limonene. J. Mol. Catal. A Chem. 2006, 252, 5–11. [Google Scholar] [CrossRef]
  64. Corma, A.; Iborra, S.; Velty, A. Chemical Routes for the Transformation of Biomass into Chemicals. Chem. Rev. 2007, 107, 2411–2502. [Google Scholar] [CrossRef] [PubMed]
  65. Marchese, A.; Arciola, C.R.; Barbieri, R.; Silva, A.S.; Nabavi, S.F.; Tsetegho Sokeng, A.J.; Izadi, M.; Jafari, N.J.; Suntar, I.; Daglia, M.; et al. Update on Monoterpenes as Antimicrobial Agents: A Particular Focus on p-Cymene. Materials 2017, 10, 947. [Google Scholar] [CrossRef] [PubMed]
  66. Salehi, B.; Upadhyay, S.; Erdogan Orhan, I.; Kumar Jugran, A.; LDJayaweera, S.; ADias, D.; Sharopov, F.; Taheri, Y.; Martins, N.; Baghalpour, N.; et al. Therapeutic Potential of α- and β-Pinene: A Miracle Gift of Nature. Biomolecules 2019, 9, 738. [Google Scholar] [CrossRef] [PubMed]
  67. Gan, C.; Langa, E.; Valenzuela, A.; Ballestero, D.; Pino-Otín, M.R. Synergistic Activity of Thymol with Commercial Antibiotics against Critical and High WHO Priority Pathogenic Bacteria. Plants 2023, 12, 1868. [Google Scholar] [CrossRef] [PubMed]
  68. Kachur, K.; Suntres, Z. The antibacterial properties of phenolic isomers, carvacrol and thymol. Crit. Rev. Food Sci. Nutr. 2020, 60, 3042–3053. [Google Scholar] [CrossRef] [PubMed]
  69. Balahbib, A.; El Omari, N.; Hachlafi, N.E.L.; Lakhdar, F.; El Menyiy, N.; Salhi, N.; Mrabti, H.N.; Bakrim, S.; Zengin, G.; Bouyahya, A. Health beneficial and pharmacological properties of p-cymene. Food Chem. Toxicol. 2021, 153, 112259. [Google Scholar] [CrossRef] [PubMed]
  70. Schaeffer, N.; Abranches, D.O.; Silva, L.P.; Martins, M.A.R.; Carvalho, P.J.; Russina, O.; Triolo, A.; Paccou, L.; Guinet, Y.; Hedoux, A.; et al. Non-Ideality in Thymol + Menthol Type V Deep Eutectic Solvents. ACS Sustain. Chem. Eng. 2021, 9, 2203–2211. [Google Scholar] [CrossRef]
  71. Bergua, F.; Castro, M.; Lafuente, C.; Artal, M. Thymol+l-menthol eutectic mixtures: Thermophysical properties and possible applications as decontaminants. J. Mol. Liq. 2022, 368, 120789. [Google Scholar] [CrossRef]
  72. Benito, C.; Alcalde, R.; Atilhan, M.; Aparicio, S. High-pressure properties of type V Natural Deep Eutectic Solvents: The case of menthol: Thymol. J. Mol. Liq. 2023, 376, 121398. [Google Scholar] [CrossRef]
  73. Kamatou GP, P.; Vermaak, I.; Viljoen, A.M.; Lawrence, B.M. Menthol: A simple monoterpene with remarkable biological properties. Phytochemistry 2013, 96, 15–25. [Google Scholar] [CrossRef] [PubMed]
  74. Galeotti, N.; Di Cesare Mannelli, L.; Mazzanti, G.; Bartolini, A.; Ghelardini, C. Menthol: A natural analgesic compound. Neurosci. Lett. 2002, 322, 145–148. [Google Scholar] [CrossRef] [PubMed]
  75. Ohno, R.; Suzuki, N.; Tsugami, Y.; Nii, T.; Kobayashi, K.; Isobe, N. Menthol application on healthy and inflamed goat udders changes antimicrobial components in milk. Anim. Sci. J. 2023, 94, e13832. [Google Scholar] [CrossRef] [PubMed]
  76. Peel, J.; John, K.; Page, J.; Jeffries, O.; Heffernan, S.M.; Tallent, J.; Waldron, M. Topical application of isolated menthol and combined menthol-capsaicin creams: Exercise tolerance, thermal perception, pain, attentional focus and thermoregulation in the heat. Eur. J. Sport Sci. 2023, 23, 2038–2048. [Google Scholar] [CrossRef] [PubMed]
  77. Davis, E.M.; Ringer, K.L.; McConkey, M.E.; Croteau, R. Monoterpene metabolism. Cloning, expression, and characterization of menthone reductases from peppermint. Plant Physiol. 2005, 137, 873–881. [Google Scholar] [CrossRef] [PubMed]
  78. Shul’pin, G.B. Metal-catalyzed hydrocarbon oxygenations in solutions: The dramatic role of additives: A review. J. Mol. Catal. A Chem. 2002, 189, 39–66. [Google Scholar] [CrossRef]
  79. Rydel-Ciszek, K.; Pacześniak, T.; Chmielarz, P.; Sobkowiak, A. Bio-Inspired Iron Pentadentate Complexes as Dioxygen Activators in the Oxidation of Cyclohexene and Limonene. Molecules 2023, 28, 2240. [Google Scholar] [CrossRef] [PubMed]
  80. Fleming, I. Molecular Orbital Theory. In Molecular Orbitals and Organic Chemical Reactions; Wiley, John Wiley & Sons, Ltd.: Hoboken, NJ, USA, 2009; pp. 1–57. [Google Scholar]
  81. Ahmed, L.; Omer, R.; Qader, I.; Koparir, P. Theoretical Analysis of the Reactivity of Carmustine and Lomustine Drugs. J. Phys. Chem. Funct. Mater. 2022, 5, 84–96. [Google Scholar]
  82. Eryilmaz, S.; Gul, M.; Inkaya, E. Investigation of global reactivity descriptors of some perillaldehyde derivatives in different solvents by DFT method. Indian J. Chem. Technol. 2019, 26, 235–238. [Google Scholar]
  83. Mali, S.N.; Anand, A.; Zaki, M.E.A.; Al-Hussain, S.A.; Jawarkar, R.D.; Pandey, A.; Kuznetsov, A. Theoretical and Anti-Klebsiella pneumoniae Evaluations of Substituted 2,7-dimethylimidazo[1,2-a]pyridine-3-carboxamide and Imidazopyridine Hydrazide Derivatives. Molecules 2023, 28, 2801. [Google Scholar] [CrossRef]
  84. Boukabcha, N.; Benmohammed, A.; Belhachemi, M.H.M.; Goudjil, M.; Yahiaoui, S.; Megrouss, Y.; Djafri, A.; Khelloul, N.; Benyehlou, Z.D.; Djafri, A.; et al. Spectral investigation, TD-DFT study, Hirshfeld surface analysis, NCI-RDG, HOMO-LUMO, chemical reactivity and NLO properties of 1-(4-fluorobenzyl)-5-bromolindolin-2,3-dione. J. Mol. Struct. 2023, 1285, 135492. [Google Scholar] [CrossRef]
  85. Parr, R.G.; Pearson, R.G. Absolute hardness: Companion parameter to absolute electronegativity. J. Am. Chem. Soc. 1983, 105, 7512–7516. [Google Scholar] [CrossRef]
  86. Pearson, R.G. Absolute electronegativity and hardness correlated with molecular orbital theory. Proc. Natl. Acad. Sci. USA 1986, 83, 8440–8441. [Google Scholar] [CrossRef] [PubMed]
  87. Pearson, R.G. Absolute electronegativity and hardness: Applications to organic chemistry. J. Org. Chem. 1989, 54, 1423–1430. [Google Scholar] [CrossRef]
  88. Parthasarathi, R.; Subramanian, V.; Roy, D.R.; Chattaraj, P.K. Electrophilicity index as a possible descriptor of biological activity. Bioorganic Med. Chem. 2004, 12, 5533–5543. [Google Scholar] [CrossRef] [PubMed]
  89. Zaklika, J.; Hładyszowski, J.; Ordon, P. From the Electron Density Gradient to the Quantitative Reactivity Indicators: Local Softness and the Fukui Function. ACS Omega 2022, 7, 7745–7758. [Google Scholar] [CrossRef]
  90. Omer, R.; Koparir, P.; Qader, I.N.; Ahmed, L. Structure reactivity analysis for Phenylalanine and Tyrosine. Cumhur. Sci. J. 2021, 42, 576–585. [Google Scholar] [CrossRef]
  91. Giang, P.D.; Churchman, L.R.; Buczynski, J.B.; Bell, S.G.; Stok, J.E.; De Voss, J.J. CYP108N14: A Monoterpene Monooxygenase from Rhodococcus globerulus. Arch. Biochem. Biophys. 2024, 752, 109852. [Google Scholar] [CrossRef] [PubMed]
  92. Neaţu, F.; Culică, G.; Florea, M.; Parvulescu, V.I.; Cavani, F. Synthesis of Terephthalic Acid by p-Cymene Oxidation using Oxygen: Toward a More Sustainable Production of Bio-Polyethylene Terephthalate. ChemSusChem 2016, 9, 3102–3112. [Google Scholar] [CrossRef]
  93. Jakaria, M.; Cho, D.Y.; Ezazul Haque, M.; Karthivashan, G. Neuropharmacological Potential and Delivery Prospects of Thymoquinone for Neurological Disorders. Oxid. Med. Cell Longev. 2018, 2018, 1209801. [Google Scholar] [CrossRef]
  94. Liu, C.; Gao, Q.; Shang, Z.; Liu, J.; Zhou, S.; Dang, J.; Liu, L.; Lange, I.; Srividya, N.; Lange, B.M.; et al. Functional Characterization and Structural Insights Into Stereoselectivity of Pulegone Reductase in Menthol Biosynthesis. Front. Plant Sci. 2021, 12, 780970. [Google Scholar] [CrossRef] [PubMed]
  95. Kani, İ.; Taşkınlar, İ.; Uzel, Z.; Avan, İ. Catalytic oxidation of thymol and carvacrol with Mn(II)-benzoylbenzoate-bipyridine complex. Polyhedron 2024, 249, 116772. [Google Scholar] [CrossRef]
  96. Frisch, M.J.; Trucks, G.W.; Schlegel, H.B.; Scuseria, G.E.; Robb, M.A.; Cheeseman, J.R.; Scalmani, G.; Barone, V.; Petersson, G.A.; Nakatsuji, H.; et al. Gaussian, version 16, revision B.01; Gaussian, Inc.: Wallingford, CT, USA, 2016. [Google Scholar]
  97. Lee, C.; Yang, W.; Parr, R.G. Development of the Colle-Salvetti correlation-energy formula into a functional of the electron density. Phys. Rev. B 1988, 37, 785–789. [Google Scholar] [CrossRef] [PubMed]
  98. Garza, A.J.; Scuseria, G.E. Predicting Band Gaps with Hybrid Density Functionals. J. Phys. Chem. Lett. 2016, 7, 4165–4170. [Google Scholar] [CrossRef] [PubMed]
  99. Moto Ongagna, J.; Tamafo Fouegue, A.D.; Ateba Amana, B.; Mouzong D’ambassa, G.; Zobo Mfomo, J.; Mbaze Meva’A, L.; Bikele Mama, D. B3LYP, M06 and B3PW91 DFT assignment of nd8 metal-bis-(N-heterocyclic carbene) complexes. J. Mol. Model. 2020, 26, 246. [Google Scholar] [CrossRef] [PubMed]
  100. Grimme, S.; Ehrlich, S.; Goerigk, L. Effect of the damping function in dispersion corrected density functional theory. J. Comput. Chem. 2011, 32, 1456–1465. [Google Scholar] [CrossRef] [PubMed]
  101. Becke, A.D. Density-functional thermochemistry. III. The role of exact exchange. J. Chem. Phys. 1993, 98, 5648–5652. [Google Scholar] [CrossRef]
  102. Trung, N.Q.; Mechler, A.; Hoa, N.T.; Vo, Q.V. Calculating bond dissociation energies of X−H (X = C, N, O, S) bonds of aromatic systems via density functional theory: A detailed comparison of methods. R. Soc. Open Sci. 2022, 9, 220177. [Google Scholar] [CrossRef]
Molecules 29 01579 sch001
Scheme 1. Possible products of limonene oxidation.
Scheme 1. Possible products of limonene oxidation.
Molecules 29 01579 sch001
Molecules 29 01579 sch002
Scheme 2. Structures of selected monoterpenes.
Scheme 2. Structures of selected monoterpenes.
Molecules 29 01579 sch002
Molecules 29 01579 g001
Figure 1. Energies of the HOMO and LUMO molecular orbitals for the limonene molecule with their visualized surfaces for the corresponding orbitals, B3LYP/6-311+G, assuming different solvent models.
Figure 1. Energies of the HOMO and LUMO molecular orbitals for the limonene molecule with their visualized surfaces for the corresponding orbitals, B3LYP/6-311+G, assuming different solvent models.
Molecules 29 01579 g001
Molecules 29 01579 g002
Figure 2. Energies of the HOMO molecular orbitals for limonene, cymene, pinene, thymol, and menthol, B3LYP/6-311+G, assuming different solvent models.
Figure 2. Energies of the HOMO molecular orbitals for limonene, cymene, pinene, thymol, and menthol, B3LYP/6-311+G, assuming different solvent models.
Molecules 29 01579 g002
Table 1. Complexes, oxidants, solvents and main oxidation products used in limonene oxidation reactions.
Table 1. Complexes, oxidants, solvents and main oxidation products used in limonene oxidation reactions.
EntryCatalystsOxidantsSolventsOxidations ProductsRef.
1MoCl2O2(Bipy)/TiO2-NTO2/(λ = 360 nm)MeCNLO, DLO, CVN[30]
2MoO2(L1–L4)/TiO2-NTO2/(λ = 360 nm)MeCNLO, DLO, CVN[31,32]
3CYP108N12O2Tris-HClPALC, PALD[9]
4[Fe(bpy)2]2+, [Fe(bpy)2]3+air, O2MeCNLO, CVN, CVL, PALD[33]
5[Mn(bpy)2]2+air, O2MeCNLO, CVN, CVL, PALD, PALC [34]
6TiO2-P25O2MeCNLO[36]
7ZnCo2O4, isobutyraldehydeO2MeCNLO[35]
8Pd(OAc)2/PTSA/BQ, Na2PdCl4/PTSA/BQO2MeOH, ethanol, 2-Propanolallylic ethers[63]
9NiAl-HTO2LimoneneLO, CVN, CVL[29]
10[MoO2(SAP)]2, [MoO2(SATP)]2t-Bu-OOHLimoneneLO, LDIOL[37]
11Ti:OTMS H2O2LimoneneCVN, CVL, LO[38]
12Na10[Co5W19O70H4]·44H2Oair, H2O2MeCN, MeOH, acetoneLO, CVN, CVL[39]
13[M4(H2O)2(PW9O34)2]n−, M–CoII, MnII, FeIIIH2O2MeCNLDIOL, CVN, CVL[2]
14[(C18H37)2N(CH3)2]3PW4O20, H2O2Tire pyrolysis oilLO, DLO, LDIOL[10]
15(nBu4N)3[NaHAsW9O33[P(O)R]2]
(R = t-Bu or CH2CH2COOH)
(n-Bu4N)3-[NaHPW9O34[As(O)p-C6H4NH2]2]		H2O2MeCNLO, DLO, LDIOL[11]
16Co(II)-Y, Cu(II)-Y with Schiff base ligandsH2O2MeCNCVN, CVL, LO, LDIOL[40]
17[Mn(acac)2APTS]@MCM-41H2O2Acetone–t-butanolLO, CVL, CVN, polymer[41]
18Al2O3H2O2Ethyl acetateLO, DLO, 8,9-LO[42]
19Ga(NO3)3, Al(NO3)3H2O2Ethyl acetate, THF LO, DLO, LDIOL, 8,9-LO[1]
20MTO:L5-7H2O2CH2Cl2LO, 8,9-LO, DLO, CVL, CVN[43]
21carbon EuroPh with FeH2O2MeOHPALC, CVL, CVN, LO, LDIOL[44]
22[VO(L8)H2O]-Y, [Cu(L8)H2O]-YH2O2MeCNLDIOL, CVL, CVN, LO[46]
23[VO(sal2bz)]2, [VO(sal2bz)]2-Y
[Fe(sal2bz)(H2O)2]2⋅2H2O, [Fe(sal2bz)(H2O)2]2-Y		H2O2MeCNLDIOL, PALC, CVN, CVL [47]
24Mn(III)/Fe(III)/Co(III)/L9/Y/ammonium acetateH2O2MeCNLO, 8,9-LO[49]
25[FeII(L10)2(H2O)2]-YH2O2MeCNCVN, CVL[52]
26RuY, FeY, 3Y–6YH2O2MeCNCVN, CVL, LO, LDIOL [51]
27γ-Fe2O3/SiO2-NHFePm-CPBA, H2O2MeCNLO, CVN, CVL[48]
28MoO2-EuroPhH2O2, t-Bu-OOH MeOHCVN, CVL, LO, PALC[45]
29[VO(VFCH)2]-Y, [VO(VTCH)2]-Y, [VO(SFCH)·H2O]-Y, [VO(STCH)·H2O]-YH2O2, t-Bu-OOHMeCNLO, LDIOL, CVN, CVL[50]
30[MoO3(Hpto)]∙H2Ot-Bu-OOHα,α,α- trifluorotolueneLO, LDIOL, DLO[53]
31[MoO3(Hpytz)]t-Bu-OOHLO, LDIOL, DLO[26]
32Ti-MCM-41, Ti-MWWt-Bu-OOHMeOHLO, CVN, CVL, PALC[54]
33FePcCl16-NH2-SiO2t-Bu-OOHAcetoneCVN, LO, CVL[55]
34cobalt(II)-(acac)-carbon-based catalystst-Bu-OOHAcetone–t-butanolLO, CVN, CVL, polymer[58]
35(η5-C9H7)Mo(CO)3Met-Bu-OOHDecane, t-butanolLO, DLO, LDIOL[56]
36MCM-41-Mn(4-OHsalen), MCM-41Mn(4-OHsalhd), MCM-41 Mn(4-OHsalophen)t-Bu-OOHAcetone–t-butanolLO, CVN, CVL, polymer[57]
37Mn(III)-Jacobsen-type catalystsKHSO5 Acetone–H2ODLO[60]
38Mn(Salen)Cl∙H2OPhIOCH2Cl2LO, CVN, PALD[59]
39M(Salen)Cl∙H2O M = MnII, NiII, CoIIPhIO, NaOClAcetone, MeCN, CH2Cl2, ethyl acetateLO, CVN, CVL[61]
40VO(Salten)-SBA-15UHPMeCNLO, CVN, CVL, carvacrol [62]
Abbreviations: Oxidants: m-CPBA—3-chloro peroxybenzoic acid; t-Bu-OOH—t-butyl hydroperoxide; PhIO—iodosylbenzene; UHP—urea hydroperoxide. Solvents: THF—tetrahydrofurane. Ligands: NT—nanotube; Bipy—2,2′-bipyridine-4,4′-dicarboxylate; L1—Schiff base; L2—bipyridine; L3—terpyridine; L4—2-aminothiazole-4-carboxylic acid; bpy—2,2′-bipyridine; PTSA—p-toluenesulfonic acid; BQ—benzoquinone; HT—hydrotalcites; SA(T)P—salicylidene amino(thio)phenolate; OTMS—octadecyltrimethoxysilane; ATPS—aminopropyltrimethoxysilane; MTO—methyltrioxorhenium; L5—t-butylpyridine; L6—4,4′-dimethyl-2,2′-bipyridine; L7—pyrazole; L8—(Z)-4-(((2hydroxyphenyl)amino)methylene)- 3-methyl-1-phenyl-1H-pyrazol-5-ol; sal2bzH2—(Z)-2-[4′-Methyleneamino-biphenyl-4-ylimino)- methyl]-phenol; L9—5,10,15,20-tetra(4-methoxyphenyl) porphyrin; L10—(Z)-2-((4-hydroxybenzylidene)-amino)benzoic acid; FeY—Na26.1Fe1.3[(AlO2)30(SiO2)162; RuY—Na27.3Ru0.9[(AlO2)30(SiO2)162; 3Y—Na27.4[Fe0.86(1)0.61(AlO2)30(SiO2)162; 4Y—Na27.12[Fe0.96(2)0.59(AlO2)30(SiO2)162; 5Y—Na28.41[Ru0.53(1)0.46(AlO2)30(SiO2)162; 6Y—Na28.62[Ru0.46(2)0.41(AlO2)30(SiO2)162; VTCH—vanillin thiophene-2-carboxylic hydrazine; VFCH—vanillin furoic-2-carboxylic hydrazone; H2STCH—salicylaldehyde thiophene-2-carboxylichydrazone; H2SFCH—salicylaldehyde furoic-2-carboxylic hydrazone; Hpto—5-(2-pyridyl-1-oxide) tetrazole; Hpytz—5-(2-pyridyl)tetrazole; Pc—phthalocyanine; 4-HOsalen)—N′-bis(4-hydroxysalicylaldehyde) ethylenediamine; 4-HOsalhd—N,N’-bis(4-hydroxysalicylaldehyde) cyclohexane-diimine; 4-HOsalophen—N,N′-bis(4-hydroxysalicylaldehyde)phenylenediamine; Salten—3-[N,N′-bis-3(salicylidenamino)ethyltriamine.
Table 2. Energy of the HOMO and LUMO orbitals along with ΔEgap levels for limonene and selected monoterpenes, structures optimized using the method without and with PCM model.
Table 2. Energy of the HOMO and LUMO orbitals along with ΔEgap levels for limonene and selected monoterpenes, structures optimized using the method without and with PCM model.
GasH2OMeCNMeOH
EHOMO, [eV]ELUMO, [eV]Egap, [eV]EHOMO, [eV]ELUMO, [eV]Egap, [eV]EHOMO, [eV]ELUMO, [eV]Egap, [eV]EHOMO, [eV]ELUMO, [eV]Egap, [eV]
limonene −6.41867−0.023676.39500−6.140570.293616.43418−6.144110.288906.43309−6.144920.287906.43282
cymene−6.45105−0.328176.12288−6.22846−0.116476.11199−6.23010−0.117836.11227−6.23037−0.118106.11227
pinene−6.207240.027486.23472−5.911450.331986.24343−5.915530.327086.24261−5.916350.326276.24262
thymol−6.08669−0.392125.69457−5.83036−0.112115.71825−5.83362−0.117015.71661−5.84614−0.127895.71825
menthol−7.27230−0.038377.23393−7.038010.370897.40890−7.042090.3635547.40564−7.042630.362197.40482
Table 3. Calculated GCRD for limonene and selected monoterpenes, structures optimized using the B3LYP/6-311+G method without and with PCM model.
Table 3. Calculated GCRD for limonene and selected monoterpenes, structures optimized using the B3LYP/6-311+G method without and with PCM model.
TerpeneE0
[a.u.]		I
[eV]		A
[eV]		X
[eV]		η
[eV]		S
[eV]		ω
[eV]		μ
[eV]
limonene−390.760169526.4190.0243.2213.1970.1561.623−3.221
cymene−389.593602626.4510.3283.3903.0610.1631.876−3.390
pinene−390.740183436.207−0.0273.0903.1170.1601.531−3.090
thymol−464.835402776.0870.3923.2392.8470.1761.843−3.239
menthol−468.448847837.2720.0383.6553.6170.1381.847−3.655
MeCNlimonene−390.761981366.144−0.2892.9283.2170.1551.332−2.928
cymene−389.595970356.2300.1183.1743.0560.1641.648−3.174
pinene−390.741046125.916−0.3272.7943.1210.1601.251−2.794
thymol−464.841565215.8340.1172.9752.8580.1751.549−2.975
menthol−468.453288997.042−0.3643.3393.7030.1351.506−3.339
H2Olimonene−390.762035106.141−0.2932.9233.2170.1551.328−2.923
cymene−389.596046006.2280.1163.1723.0560.1641.647−3.172
pinene−390.741070985.911−0.3322.7903.1220.1601.247−2.790
thymol−464.841762915.8300.1122.9712.8590.1751.544−2.971
menthol−468.453416937.038−0.3713.3343.7050.1351.500−3.334
MeOHlimonene−390.761972216.145−0.2882.9293.2160.1551.333−2.929
cymene−389.595957516.2300.1183.1743.0560.1641.648−3.174
pinene−390.741041885.916−0.3262.7953.1210.1601.251−2.795
thymol−464.841537395.8460.1282.9872.8590.1751.560−2.987
menthol−468.453267217.043−0.3623.3403.7020.1351.507−3.340
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Rydel-Ciszek, K. DFT Studies of the Activity and Reactivity of Limonene in Comparison with Selected Monoterpenes. Molecules 2024, 29, 1579. https://doi.org/10.3390/molecules29071579

AMA Style

Rydel-Ciszek K. DFT Studies of the Activity and Reactivity of Limonene in Comparison with Selected Monoterpenes. Molecules. 2024; 29(7):1579. https://doi.org/10.3390/molecules29071579

Chicago/Turabian Style

Rydel-Ciszek, Katarzyna. 2024. "DFT Studies of the Activity and Reactivity of Limonene in Comparison with Selected Monoterpenes" Molecules 29, no. 7: 1579. https://doi.org/10.3390/molecules29071579

Article Metrics

Back to TopTop