Volatile Components, Antioxidant and Phytotoxic Activity of the Essential Oil of Piper acutifolium Ruiz & Pav. from Peru

Piper acutifolium Ruiz & Pav is known as “matico” and belongs to the Piperaceae family, and in Peru it is traditionally used as an infusion or decoction to ameliorate wound healings or ulcers. In this study, the aim was to investigate the volatile components, the antioxidant profile, and the phytotoxic activity of the essential oil (EO) of P. acutifolium from Peru. To identify the phytoconstituents, the EO was injected into a Gas Chromatography-Mass Spectrometry (GC-MS) to obtain the chemical profile of the volatile components, followed by the antioxidant activity carried out by the reaction with three organic radicals (2,2-diphenyl-1-picrylhydrazyl (DPPH); 2,2′-azinobis-(3-ethylbenzothiazoline)-6- sulfonic acid (ABTS); ferric reducing/antioxidant power (FRAP)). Finally, the phytotoxic capabilities of the EO were tested on two model plants, Lactuca sativa seeds and Allium cepa bulbs. As a result, the analysis identified α-phellandrene as its main volatile chemical at 38.18%, followed by β-myrcene (29.48%) and β-phellandrene (21.88%). Regarding the antioxidant profile, the half inhibitory concentration (IC50) in DPPH was 160.12 ± 0.30 µg/mL, for ABTS it was 138.10 ± 0.06 µg/mL and finally in FRAP it was 450.10 ± 0.05 µg/mL. The phytotoxic activity demonstrated that the EO had high activity at 5% and 10% against L. sativa seed germination, the inhibition of root length, and hypocotyl length. Additionally, in A. cepa bulbs, the inhibition root length was obtained at 10%, both comparable to glyphosate, which was used as a positive control. The molecular docking on 5-enolpyruvylshikimate-3-phosphate synthase (EPSPS) revealed that α-phellandrene had −5.8 kcal/mol, being near to glyphosate at −6.3 kcal/mol. The conclusion shows that the EO of P. acutifolium presented antioxidant and phytotoxic activity and might be useful as a bioherbicide in the future.


Introduction
Essential oils (EOs) are a kind of natural product obtained from the different organs of aromatic plants (seeds, leaves, fruits, stems, roots) by different methods, both conventional and unconventional, such as steam distillation, hydrodistillation, solvent extraction, supercritical fluids, and microwaves, among others. [1]. The chemical composition of EOs varies according to the plant genus, soil type, altitude, plant part extracted, growth conditions, etc. [2]. The composition of EOs generally contains compounds identified as monoterpenes and sesquiterpenes which have demonstrated several biological activities such as insecticidal, antitumor, cytotoxic, antimicrobial, anti-inflammatory, hepatoprotective, gastroprotective, anxiolytic, antidepressant, phytotoxic, and analgesic, etc. [3].
In recent decades, a growing number of phytosanitary pesticides have been administered to various agricultural regions. There is research on the possible use of essential oils as bio-pesticides in sustainable agronomy [4]. The Food and Drug Administration of the USA includes EOs in the category of Generally Recognized as Safe [5]. On the contrary, essential oils rich in monoterpenes have shown toxicity assessed in vitro and in vivo studies, so their use should be strictly controlled by rigorous toxicity studies. The phytotoxic effect of EOs has been associated with the presence of volatile compounds such as carvacrol, camphor, α-pinene, 1,8-cineole, limonene and thymol, which have varied individual phytotoxicity levels [6].
Regarding the use of EOs as a bioherbicide, several reports consider its use to treat invasive plants that contaminate crops. The herbicides can be selective (2,4-dichlorophenoxyacetic acid), non-selective (glyphosate), or work during different stages of plant development [7]. Numerous biochemical processes contribute to its herbicidal activity, including the disruption of the cuticle and damage to young tissues, the inhibition of photosynthesis and mitochondrial respiration, the change in enzymatic and phytohormone regulation, the modification of membrane properties and interactions, and the induction of microtubule disruption, genotoxicity, and reactive oxygen and nitrogen species. [8].
Antioxidants are molecules that can react with free radicals, reducing their oxidizing power, causing oxidative stress and slowing or retarding oxidation. These oxidizing agents are very harmful when there is an imbalance in the oxidation-reduction system and can affect lipids, proteins, and carbohydrates in organisms [9]. Regarding the antioxidants and the phytotoxic activity, an inverse relationship exists because the toxicity induced in weeds by bioherbicides can be directly associated with the generation of reactive oxygenated species to produce damage [10]. However, in recent years, the antioxidant activity of essential oils has been studied to prevent food from oxidation during storage, in food stabilizers, and in active packaging and edible coatings in the food industry [11]. Thus, the investigation of antioxidants in essential oils seems to be focused on the inhibition of lipid peroxidation, free radicals, or chelating metal ions, which are processes involved with several pathological conditions in humans. Furthermore, depending on its extraction method and other pretreatments prior to the distillation process, the antioxidant activity can provide different results [12].
In Peru, Piper acutifolium Ruiz & Pav (Piperaceae Family) is commonly known as matico, and a decoction of its leaves (approximately 50 g per liter of water) is used as a wash for sores and ulcers, as well as for vaginal infections [13]. It is also recommended for gastritis and menstrual disorders, with an oral administration of two cups per day [14]. According to Lognay et al., its essential oil is abundant in β-ocimene, α-copaene, allo-aromadendrene, α-cadinene, δ-cadinene, myristicin, and dillapiole [15]. A study showed that benzoic acid derivates isolated from the dichloromethane extract had activity against Leishmania spp., Trypanosoma cruzi, and Plasmodium falciparum [16]. There are no reported studies based on its antioxidant and its phytotoxic activity. Therefore, it is necessary to characterize the essential oil obtained by gas chromatography coupled to mass spectrometry (GC-MS) and to determine the antioxidant activity in vitro, and the phytotoxic activity in lettuce seeds and onion bulbs, which might be useful as a potential bioherbicide in the future.

Chemical Characterization of the Essential Oil of P. acutifolium
The obtained EO of P. acutifolium presented a yellow color (golden) with an extraction yield of 0.067%, a relative density of 0.9834 ± 0.03 g/mL at 20 • C, and a refractive index of 1.5000 at 20 • C. The profile chemical identified by GC-MS is reported in Table 1.

Antioxidant Activity of the Essential Oil of P. acutifolium
According to Table 2, the results show the antioxidant activity using three organic radicals (DPPH, ABTS and FRAP). As presented, the best method to demonstrate antioxidant Molecules 2023, 28, 3348 5 of 17 activity was in ABTS with an IC 50 of 138.10 ± 0.06 µg/mL and a TEAC of 150.56 ± 8.27 µmol TE/g of EO, followed by DPPH and FRAP. Studies of antioxidant activity using the evaluated methods reveal that in non-polar plant extracts, the correlation between DPPH and ABTS is limited, and decreased in comparison to polar extracts [25]. Other Piper species showed higher levels of antioxidant activity, such as P. auritum (IC 50 = 14.8 ± 0.5 µg/mL), P. aduncum (IC 50 = 30.1 ± 1.8 µg/mL), and P. umbellatum (IC 50 = 32.3 ± 1.3 µg/mL) against DPPH radicals [12,21,26]. Bagheri et al. [12] demonstrated that the EO obtained by hydro distillation had better antioxidant activity than supercritical CO 2 extraction in P. nigrum against DPPH, with an IC 50 of 103.28 µg/mL and 316.27 µg/mL, respectively. Alminderej et al., [27] showed that the essential oil of P. cubeba fruit had an antioxidant activity against DPPH and FRAP of IC 50 = 110.00 ± 0.08 µg/mL, and 106.00 ± 0.11 µg/mL respectively. Valerazo et al., [28] reported the antioxidant activity (IC 50 ) of P. ecuadorense, with values of more than 2.5 mg/mL for DPPH and 1.81 ± 0.09 mg/mL for the ABTS . + radical. Additionally, P. brachypetiolatum had an IC 50 value of 64.8 µg/mL and P. madeiranum had a value of 66.8 µg/mL against the DPPH radical [29]. In this study, the essential oil from P. acutifolium had low antioxidant activity (Table 2), which could be due to the presence of monoterpene hydrocarbons in high percentages as some essential oils with high antioxidant activity contain oxygenated monoterpenes or sesquiterpenes in their composition [30].

Phytotoxic Activity of the Essential Oil of P. acutifolium on L. sativa Seeds and A. cepa Bulbs
The use of bio-herbicides based on natural products, such as essential oils, is generating much interest in the world with regard to environmental concerns. Glyphosate is a commonly used synthetic herbicide that blocks the enzyme 5-enol-pyruvyl-shikimate-3-phosphate synthase (EPSPS) and catalyzes the sixth step in the shikimic acid pathway, reducing the aromatic amino acids such as phenylalanine, tyrosine, and tryptophan [31]. It is used as a positive control in this study, and its phytotoxic activity is presented in Table 3, Figures 1 and 2, with the inhibition of seed germination, root length, hypocotyl length, and the rate of root length/stem length. In this study, the evaluation of P. acutifolium EO on L. sativa seeds (Table 3) showed phytotoxicity activity at 5% and 10%, comparable to 2% glyphosate. P. acutifolium EO negatively affected the germination of L. sativa seeds (Figure 1), depending on the concentration, while the negative control (0.1% DMSO) did not show any phytotoxic activity. However, germination was inhibited at the highest concentration, being similar to the positive control (2% glyphosate). The root and hypocotyl length of L. sativa was inhibited when exposed to 5% and 10% EO solutions compared with the negative control (Table 3). It is known that some volatile terpenes have revealed phytotoxic activity against L. sativa seeds including α-pinene, γ-terpinene, and p-cymene from Eucalyptus grandis [32]; α-pinene, and β-pinene from Pinus brutia [33], 1,8-cineole, βphellandrene, and α-pinene from Majorana hortensis, o-cymene, and α-pinene from Thymus vulgaris, estragole, limonene, and β-pinene from Carum carvi [34], and β-caryophyllene, (Z)-caryophyllene, and germacrene D from Ailanthus altissima [35]. There was no significant phytotoxic effect between P. acutifolium EO at 5% (0.93 ± 0.12 cm; p = 0.9983) and 10% (0.47 ± 0.06 cm; p = 0.7172) with glyphosate (0.83 ± 0.15 cm) on A. cepa bulbs ( Figure 3); the other concentrations had values of 6.33 ± 0.58 cm at 0.1% (p < 0.0001), 4.5 ± 0.6 cm (p < 0.0001) at 0.5%, and 1.93 ± 0.12 cm (p = 0.0151) at 1% respectively. Studies of A. cepa bulbs indicate that terpenes have a phytotoxic effect, such as β-pinene, δ-carene, and limonene from Heterothalamus psiadioides, and β-caryophyllene and spathulenol from Baccharis patens [35]. Abd-ElGawad et al. concluded that the phytotoxic activity of the EOs are linked to its content of oxygenated terpenoids. However, some sesquiterpenes like β-caryophyllene and derivates have demonstrated phytotoxic action, and in this study only aromadendrene at 0.78% was identified, and this oxygenated sesquiterpene might have a synergism phytotoxic activity on L. sativa and A. cepa bulbs. In addition, monoterpenes and sesquiterpenes extracted from EOs have shown phytotoxic effects causing anatomical and physiological changes in plants, some proposed mechanisms such as the accumulation of lipid globules in the cytoplasm, oxidative stress, the reduction of mitochondria, and the inhibition of DNA synthesis might be involved in its phytotoxicity. of A. cepa bulbs indicate that terpenes have a phytotoxic effect, such as β-pinene, δ-carene, and limonene from Heterothalamus psiadioides, and β-caryophyllene and spathulenol from Baccharis patens [35]. Abd-ElGawad et al. concluded that the phytotoxic activity of the EOs are linked to its content of oxygenated terpenoids. However, some sesquiterpenes like βcaryophyllene and derivates have demonstrated phytotoxic action, and in this study only aromadendrene at 0.78% was identified, and this oxygenated sesquiterpene might have a synergism phytotoxic activity on L. sativa and A. cepa bulbs. In addition, monoterpenes and sesquiterpenes extracted from EOs have shown phytotoxic effects causing anatomical and physiological changes in plants, some proposed mechanisms such as the accumulation of lipid globules in the cytoplasm, oxidative stress, the reduction of mitochondria, and the inhibition of DNA synthesis might be involved in its phytotoxicity.

Molecular Docking Studies
Molecular docking studies were carried out in order to determine the binding affinities of α-phellandrene, β-myrcene, β-phellandrene and glyphosate with EPSP synthase protein. The binding energies (ΔG in kcal/mol) are listed in Table 4. The best ligand poses belonging to the population were found to be the highest in the cluster of RMSD tolerance of 0.5 and the lowest binding energy that were chosen for the docking

Molecular Docking Studies
Molecular docking studies were carried out in order to determine the binding affinities of α-phellandrene, β-myrcene, β-phellandrene and glyphosate with EPSP synthase protein.
The binding energies (∆G in kcal/mol) are listed in Table 4. The best ligand poses belonging to the population were found to be the highest in the cluster of RMSD tolerance of 0.5 and the lowest binding energy that were chosen for the docking study. A total of 95% of the ligand pose population belonging to the lowest binding energy cluster showed that the α-phellandrene occupied the deep core of the binding pocket and was involved in alkyl interactions with Ala174, Arg200, and Ile325 residues (Table 4). Rest residues were involved in non-bonded van Der Waal's interactions, as displayed in Figure 2A. Similar interactions were also observed for beta-myrcene, where alkyl interactions were observed with Lys28, Ala174, and Arg200 ( Figure 2B). β-phellandrene exhibited alkyl interactions with Ala174, Arg200, and Ile325 (Table 3, Figure 2C). On the other site, α-phellandrene exhibited the lowest binding energy ∆G = −5.8 kcal/mol and the highest affinity for the protein. The standard ligand glyphosate was involved in conventional hydrogen bonds with Lys28, Thr101, Ala100, Arrg405, and Arg357 residues and presented a high affinity and a binding energy ∆G = −6.3 kcal/mol (Table 3, Figure 2D).

Molecular Dynamics Simulation
Molecular dynamics and simulation studies were carried out to determine the stability and convergence of 3FJZ-glyphosate, 3FJZ-α-phellandrene, 3FJZ-β-phellandrene, and 3FJZβ-myrcene. The simulation of 200 ns displayed a stable conformation while comparing the root mean square deviation (RMSD) values. The RMSD of Cα-backbone of 3FJZ-glyphosate exhibited a deviation of 3.0 Å, while α-phellandrene bound protein exhibited RMSD 2.5 Å ( Figure 4A). RMSD of 3FJZ-β-phellandrene exhibited 2.7 Å ( Figure 4A), while 3FJZ-βmyrcene exhibited unstable RMSD up to 60 ns and later stabilized at 2.9 Å ( Figure 4A). All of the RMSD values are within the acceptable range, which should be below or near 3 Å. A stable RMSD plot during simulation signifies a good convergence and stable conformations. Therefore, it can be suggested that α-phellandrene and β-phellandrene bound to 3FJZ are quite stable in complex due to the higher affinity of the ligand. The plot for root means square fluctuations (RMSF) displayed large spikes of fluctuation in 3FJZ-glyphosate are observed at 10-16, 105-110 and 152-160 residues, which may be due to the higher flexibility of the residues ( Figure 4B). 1M17 with compound 2 also exhibited fluctuations at 20-54, 125-146, 150, 175-182, 230-240 and 220-230 residue positions ( Figure 4B), while 3FJZ-α-phellandrene exhibited fewer significant fluctuations except at 10-30, 110-120, 150, 175-182, 230-240 and 220-230 residues ( Figure 4B). The flexibility of the residues was observed to be more in the case of α-phellandrene. For β-phellandrene and β-myrcene, similar peaks were observed to that of glyphosate bound complex ( Figure 4B). Most of the residues fluctuated less during the entire 200 ns simulation, indicating the rigid amino acid conformations during the simulation time. Fewer fluctuations indicate a less active protein, since more flexibility allows the protein to become functionally active. Therefore, based on the RMSF plots it can be suggested that the α-phellandrene bound protein is more flexible and becomes active during simulation in ligand-bound conformations. The radius of gyration (Rg) is the measure of the compactness of the protein. Here in this study, the 3FJZ-glyphosate Cα-backbone displayed a lowering of Rg from 19.0 to 18.5 Å ( Figure 4C).
On the other hand, a lowering and stable pattern was also observed for α-phellandrene, β-phellandrene and β-myrcene bound to 3FJZ from 19.0 to 18.5 and 18.6 Å, respectively, ( Figure 4C). The initial lowering of Rg indicates a highly compact orientation of the protein in a ligand-bound state. A number of hydrogen bonds between the protein and the ligand suggest the significant interaction and stability of the complex. The number of hydrogen bonds between 3FJZ-glyphosate showed significant (4) numbers, and with α-phellandrene, β-phellandrene and β-myrcene, a single hydrogen bond was observed on an average ( Figure 4D) throughout the simulation time of 200 ns. Followed by Rg analysis, the pattern was also observed in the solvent accessible surface area (SASA) in both the ligand-bound and unbound states. The unbound state of ligands to the 3FJZ protein ( Figure 4F-H) displayed a high surface area accessible to solvents in all the cases ( Figure 4E-H). The SASA value lowered as compared to the unbound state when bound with a ligand (Figure 4E-H). The study of Rg means that the binding of ligands forces the respective proteins to become more compact.

Molecular Mechanics Generalized Born Surface Area (MM-GBSA) Calculations
Utilizing the MD simulation trajectory, the binding free energy along with other contributing energy in the form of MM-GBSA were determined for each 3FJZ-glyphosate, 3FJZ-α-phellandrene, 3FJZ-β-phellandrene and 3FJZ-β-myrcene. The results (Table 5) suggested that the maximum contribution to ∆G bind in the stability of the simulated complexes were due to ∆G bind Coulomb, ∆G bind vdW and ∆G bind Lipo, while, ∆G bind Covalent and ∆G bind SolvGB contributed to the instability of the corresponding complexes. The 3FJZ-glyphosate, 3FJZ-α-phellandrene complexes showed comparatively higher binding free energies. These results supported the potential of glyphosate and α-phellandrene molecules with 3FJZ, and showed the efficiency in binding to the selected protein and the ability to form stable protein-ligand complexes. Table 5. Binding free energy components for the 3FJZ-glyphosate, 3FJZ-α-phellandrene, 3FJZ-βphellandrene and 3FJZ-β-myrcene calculated by MM-GBSA.

In-Silico Toxicological Study of P. acutifolium EO
The human hepatotoxicity (HH), carcinogenicity and Ames's toxicity analysis are presented in Table 5. Compounds 1,4,5,6,7,8,9,11,12, and glyphosate had a good safety profile. However, compounds 2, 3 and 10 might be hepatotoxic. In addition, compounds 2, 5, 8 and 9 might be carcinogenic. In a study, the oral administration of β-myrcene (#2) induced nonneoplastic lesions in the kidney of rats [36]; with regard to other compounds, there are no toxicological reports. The results for the environmental toxicity predictions of the volatile compounds of the EO fall within 1000 L/kg < BCF < 5000 L/kg categorized as bioaccumulative by the United States Environmental Protection Agency under the Toxic Substances Control Act [37]. However, glyphosate showed a low value (BCF = 0.151), which indicates that it is not bioaccumulative. Regarding the other parameters, such as IGC 50 , LC 50 FM, and LC 50 DM, the volatile components had higher values than glyphosate, indicating that these compounds isolated in the EO of P. acutifolium might have less toxicity when compared with glyphosate. According to the data in Table 6, some compounds can be considered toxic in silico, but the use of EO as a bioherbicide has an advantage compared to synthetic herbicides; for example, they are volatile and rapidly degrade in the ecosystem [38]. However, data on its toxicity have yet to be reported, and it will depend mainly on its chemical composition and in vivo tests. Table 6. In-silico analysis of toxicity of the volatile compounds of P. acutifolium and glyphosate. List of volatile compounds: 1. α-pinene; 2. β-myrcene; 3. α-phellandrene; 4. o-cymene; 5. β-phellandrene; 6. α-gurjunene; 7. β-caryophyllene; 8. humulene; 9. germacrene D; 10. bicyclogermacrene; 11. δ-cadinene; 12. aromadendrene.

Toxicity
Environmental Toxicity

Plant Material
A quantity of 9 kg of the P. acutifolium leaves were collected from Cajabamba, Cajamarca, Peru, located at 2658 masl (7 • 37 25" S and 78 • 02 39" W)) in December 2022. The material plant was authenticated by biologist Jose Campos (certificate Id. 09092022-JP). Dried and cleaned leaves (yield = 90%) were incorporated in a Clevenger apparatus to obtain the essential oil by steam distillation for 2 h. Finally, the essential oil was separated from the hydrolate by decantation, and then anhydrous Na 2 SO 4 was added to purify the EO. Finally, the EO was kept in a sealed amber vial at 4 • C.

Identification of the Volatile Compounds by Gas Chromatography-Mass Spectrometry (GC-MS)
The volatile components of the essential oil of P. acutifolium leaves were identified using a Gas Chromatograph coupled to Mass Spectrophotometry (GCMS-QP2010 SE, Shimadzu) with a Quadrex capillary column (007 METHYL 5%) (30 m × 0.25 mm, 0.20 µm film thickness, Agilent Technologies). The EO was diluted in chloroform HPLC grade (1000 µL) (1: 100). The chromatographic conditions were set according to a previous analysis protocol [39]. The identification of volatile components was based on a comparison of relative retention indices (RIs) and mass spectra data (NIST-5 library) and published literature. Each Retention Index was calculated in comparison to a homologous series of n-alkanes C9-C25 (C9, 99% BHD purity and C10-C25, 99% Fluka purity). The EO of P. acutifolium (150 µL) at different concentrations (0.1; 0.5; 1.0; 5.0 and 10.0%) reacted with 850 µL FRAP radical's solution, which was elaborated with acetate buffer pH 3.6, and 2,4,6-Tris(2-pyridyl)-s-triazine diluted in HCl, and ferric chloride hexahydrate. The reaction was stopped at 4 min and the absorbances were read at 593 nm. Each test was carried out in triplicate [42].

Phytotoxicity Test on Lactuca sativa Seeds
Germination and growth bioassays were performed in Petri dishes (90 mm in diameter), with Whatman filter paper number 42. Commercial seeds of Lactuca sativa (lettuce) were treated with different solutions of the EO solubilized in 0.1% DMSO at 0.1; 0.5; 1.0; 5 and 10%, respectively. The Petri dishes were then incubated at 22 • C for 5 days in the absence of light to simulate the germination process. Several parameters were evaluated, including the number of germinated seeds, the root length, and the hypocotyl length. A solution of 0.1% DMSO was used as a negative control, and a 2% glyphosate solution (Fuego ® , NeoAgro, Lima, Peru) was used as a positive control. Each Petri dish had twenty seeds, and four replicates were utilized [43].

Phytotoxicity Test on Allium cepa Bulbs
A. cepa phytotoxicity assay was carried out according to Haq et al. [44]. A. cepa bulbs weighing between 70-90 g were obtained from a local market for this study. Three bulbs were cut at the base to generate root growths. Finally, A. cepa bulbs were placed in 50 mL Falcon tubes containing different concentrations of the EO (0.1, 0.5, 1.0, 5.0 and 10.0%). The negative control was 0.1% DMSO, which was used as the diluent, and a 2% glyphosate was used as a positive control. Tubes were kept at room temperature for 5 days. After 5 days, the root lengths were measured in each evaluated concentration. The phytotoxicity assay was carried out in triplicate.

Molecular Docking
This study's purpose is to provide insight into the binding propensities of the main volatile components of P. acutifolium which represented the major percentages in the total composition, and these compounds were α-phellandrene, β-myrcene, and β-phellandrene, with 5-enolpyruvylshikimate-3-phosphate synthase (EPSPS) synthase (PDB ID: 3FJZ) [45]. Additionally, glyphosate was used as the standard to validate the ligand binding. Autodock vs. 4.2.6 was used to perform docking investigations on the molecules in question. The main volatile components and the protein EPSPS were saved in pdbqt format after combining non-polar hydrogens. Molecular docking was analyzed within a grid box dimension of 20.04 × 21 × 20.5 Å. The binding pocket residues were determined to form the adjacent residues at 4Å radius of the co-crystallized ligand [45].
3.9. In-Silico Toxicity of the Volatile Components of P. acutifolium The ADMETlab 2.0 program (http://admetmesh.scbdd.com, accessed on 10 March 2023) was used to estimate the toxicity parameters and environmental toxicity. The parameters evaluated were human hepatotoxicity, the bioconcentration factor expressed as −log10[(mg/L)/(1000*MW)], the IGC 50 : concentration of the test chemical in water (mg/L) that causes 50% growth inhibition to the species Tetrahymena pyriformis after 48 h, the LC 50 FM that refers to the concentration of the test chemical in water (mg/L) that causes 50% of fathead minnow to die after 96 h, and the LC 50 DM: concentration of the test chemical in water (mg/L) that causes 50% of the species Daphnia magna to die after 48 h. These three parameters (IGC 50 , LC 50 FM, and LC 50 DM) were also expressed as −log10[(mg/L)/(1000*MW)]. To interpret the hepatotoxic, Ames' test, and the carcinogenicity test results, they were evaluated based on probability values, where values near to 1 are highly toxic [46].

Molecular Dynamics Simulation
Molecular dynamics simulation studies of 3FJZ-glyphosate, 3FJZ-α-phellandrene, 3FJZ-β-phellandrene and 3FJZ-β-myrcene were carried out in Desmond [47] for a 200 nanoseconds time scale. Docking experiments were the first step in creating protein and ligand complexes for the molecular dynamics simulation. In static settings, molecular docking studies can predict the ligand binding state. Since docking provides a static representation of a molecule's binding pose in a protein's active site, MD simulations typically compute atom motions over time by integrating Newton's classical equation of motion [48]. The ligand binding status in the physiological environment was predicted using simulations. The Protein Preparation Wizard was used to pre-process the proteinligand complex, which included complex optimization and minimization. The System Builder tool was used to create all of the systems. As transferable intermolecular interaction potential 3 points, the solvent model with an orthorhombic box was chosen. The simulation used the OPLS 2005 force field. Counter ions were used to make the models neutral. An amount of 0.15 M salt (NaCl) was added to simulate physiological circumstances. For the entire simulation, the NPT ensemble with 300 K temperature and 1 atm pressure was used. Before the simulation, the models were relaxed [49]. The trajectories were saved for evaluation every 100 ps, and the simulation's stability was confirmed by comparing the RMSD, Rg, hydrogen bonds, and the solvent accessible surface area of the protein and ligand over time.

Binding Free Energy Analysis
The molecular mechanics and generalized Born surface area method (MM-GBSA) was utilized to determine the binding free energy of standard glyphosate, α-phellandrene, β-phellandrene, and β-myrcene complexes, respectively. Using the Python script thermal mmgsba.py and the VSGB solvation model, along with the OPLS5 force field and one-step sample size, the MM-GBSA free energy of binding was measured in the final fifty frames of the simulation trajectory. The individual energy modules of columbic, covalent, hydrogen bond, van Der Waals, self-contact, lipophilic, and solvation of ligand and protein were aggregated using the concept of additivity in order to determine the binding free energy of MM-GBSA (kcal/mol). To determine the ∆G bind , the following equation was used: ∆G bind is binding free energy (kcal/mol), ∆G MM designates free energy differences of the 3FJZ + ligand complex and the total energies of 3FJZ and ligands in isolated form, ∆G Solv is the solvation energy of the ligand-receptor complex and the sum of the solvation energies of the 3FJZ and ligands in the unbound state, and ∆G SA is the surface area energy differences between 3FJZ and ligands.

Statistical Analysis
Graph Pad Prism v4.0 was used to introduce the results obtained in this study. An ANOVA test and Dunnett's post-test were used as statistical tests, and a p value more than 0.001 was considered statistically significant in the phytotoxicity test. The Trolox equivalent antioxidant capacity (TEAC) was expressed as µmol, and the Trolox equivalent/g of EO and linear regression using Microsoft Excel was used to calculate the half inhibitory concentration (IC 50 ).

Conclusions
Despite the fact that some synthetic herbicides have harmful side effects, they are becoming vital for effective weed management. Recent interest and demand for organic fruits, vegetables, dairy products, and beverages throughout the world, especially in industrialized nations, has led to an interest in bioherbicides and the search for products with lower harmful side effects. It is concluded that the EO presented 14 compounds and as the main fingerprint to β-myrcene, which was determined by GC-MS. Regarding the antioxidant activity against DPPH, ABTS and FRAP, the best result was using the ABTS radical and its antioxidant profile compared to other Piper species reported in the literature. Furthermore, the phytotoxic activity revealed that at 5% and 10%, the EO inhibited several parameters in the development of Lettuce seeds and in onion bulbs, such as germination and root length. The glyphosate control used in this study also showed phytotoxic activity, and its main target, known as 5-enolpyruvylshikimate-3-phosphate synthase (EPSPS), was inhibited by the α-phellandrene, being similar to glyphosate. The molecular dynamic showed that α-phellandrene was very stable during 200 ns of evaluation. The in silico toxicological evaluation showed that all compounds are safe for the environment, but are bioaccumulative. As such, the EO of P. acutifolium could be used as a bioherbicide in the future. In addition, an in vivo toxicological evaluation, and an environmental study of the EO is needed to assess its safety.

Data Availability Statement:
The data is available upon reasonable request.

Conflicts of Interest:
The authors declare that they have no conflict of interest.