Antiviral Potentialities of Chemical Characterized Essential Oils of Acacia nilotica Bark and Fruits against Hepatitis A and Herpes Simplex Viruses: In Vitro, In Silico, and Molecular Dynamics Studies

Acacia nilotica (synonym: Vachellia nilotica (L.) P.J.H.Hurter and Mabb.) is considered an important plant of the family Fabaceae that is used in traditional medicine in many countries all over the world. In this work, the antiviral potentialities of the chemically characterized essential oils (EOs) obtained from the bark and fruits of A. nilotica were assessed in vitro against HAV, HSV1, and HSV2. Additionally, the in silico evaluation of the main compounds in both EOs was carried out against the two proteins, 3C protease of HAV and thymidine kinase (TK) of HSV. The chemical profiling of the bark EOs revealed the identification of 32 compounds with an abundance of di- (54.60%) and sesquiterpenes (39.81%). Stachene (48.34%), caryophyllene oxide (19.11%), and spathulenol (4.74%) represented the main identified constituents of bark EO. However, 26 components from fruit EO were assigned, with the majority of mono- (63.32%) and sesquiterpenes (34.91%), where trans-caryophyllene (36.95%), Z-anethole (22.87%), and γ-terpinene (7.35%) represented the majors. The maximum non-toxic concentration (MNTC) of the bark and fruits EOs was found at 500 and 1000 µg/mL, respectively. Using the MTT assay, the bark EO exhibited moderate antiviral activity with effects of 47.26% and 35.98% and a selectivity index (SI) of 2.3 and 1.6 against HAV and HSV1, respectively. However, weak activity was observed via the fruits EO with respective SI values of 3.8, 5.7, and 1.6 against HAV, HSV1, and HSV2. The in silico results exhibited that caryophyllene oxide and spathulenol (the main bark EO constituents) showed the best affinities (ΔG = −5.62, −5.33, −6.90, and −6.76 kcal/mol) for 3C protease and TK, respectively. While caryophyllene (the major fruit EO component) revealed promising binding capabilities against both proteins (ΔG = −5.31, −6.58 kcal/mol, respectively). The molecular dynamics simulation results revealed that caryophyllene oxide has the most positive van der Waals energy interaction with 3C protease and TK with significant binding free energies. Although these findings supported the antiviral potentialities of the EOs, especially bark EO, the in vivo assessment should be tested in the intraoral examination for these EOs and/or their main constituents.

Acacia nilotica (common name: Gum Arabic tree) is an important traditional plant for the treatment of several diseases in several countries around the world [22]. Many components were isolated and identified from the extracts of the different organs of this plant comprising the terpenes, tannins, phenolic acids, and flavonoids [21,22,28,29]. Several reports described the interesting biological potentialities of the different extracts of this plant, such as the treatment of inflammation, free radicals, leishmanial, diabetes, cancers, plasmodial, and other infections, along with molluscicidal activity [30][31][32].
Little studies were carried out concerning the composition and pharmaceutical applications of A. nilotica EOs [24]. Thus, the current work aimed to (i) characterize the chemical profile of the EOs derived from the A. nilotica bark and fruits depending upon GC-MS techniques, (ii) evaluate the antiviral efficiency of these two EOs against the HSV1, HSV2, and HAV viruses, and (iii) study of the molecular docking of the main compounds in each oil on 3C protease of HAV (PDB ID: 1QA7) and thymidine kinase of HSV (PDB ID: 1KI3) proteins.

The Identification of the Chemical Constituents of Bark and Fruits of A. nilotica EOs
The EOs of the bark and fruits of A. nilotica were separately obtained via three hours of hydro-distillation over the Clevenger apparatus, where they yielded 0.072% and 0.056% (v/w) of the oil, respectively. The yield of the EO from bark was found to be more than the fruits, and both were more varied than those described for the Nigeria ecospecies [24]. The oils were analyzed via GC-MS ( Figure 1). albida) [24], A. tortilis (synonym: Vachellia tortilis (Forssk.) Galasso and Banfi) [25], A. mearnsii De Wild. (synonym: A. mearnsii) [26], and A. cyanophylla (synonym: Acacia saligna (Labill.) H.L.Wendl.) [27].
Acacia nilotica (common name: Gum Arabic tree) is an important traditional plant for the treatment of several diseases in several countries around the world [22]. Many components were isolated and identified from the extracts of the different organs of this plant comprising the terpenes, tannins, phenolic acids, and flavonoids [21,22,28,29]. Several reports described the interesting biological potentialities of the different extracts of this plant, such as the treatment of inflammation, free radicals, leishmanial, diabetes, cancers, plasmodial, and other infections, along with molluscicidal activity [30][31][32].
Little studies were carried out concerning the composition and pharmaceutical applications of A. nilotica EOs [24]. Thus, the current work aimed to (i) characterize the chemical profile of the EOs derived from the A. nilotica bark and fruits depending upon GC-MS techniques, (ii) evaluate the antiviral efficiency of these two EOs against the HSV1, HSV2, and HAV viruses, and (iii) study of the molecular docking of the main compounds in each oil on 3C protease of HAV (PDB ID: 1QA7) and thymidine kinase of HSV (PDB ID: 1KI3) proteins.

The Identification of the Chemical Constituents of Bark and Fruits of A. nilotica EOs
The EOs of the bark and fruits of A. nilotica were separately obtained via three hours of hydro-distillation over the Clevenger apparatus, where they yielded 0.072% and 0.056% (v/w) of the oil, respectively. The yield of the EO from bark was found to be more than the fruits, and both were more varied than those described for the Nigeria ecospecies [24]. The oils were analyzed via GC-MS ( Figure 1).  The compounds' names, retention times, relative concentrations, and Kovats indexes are presented in Table 1. The analysis revealed that the terpenes represented the main constituents of the bark with a relative concentration of 95.25% of the hydrocarbons and oxy- genated forms of mono-, sesquit-, and di-terpenes. Furthermore, the terpenes represented the major components of the fruits EO with a relative concentration of 98.23%, including the mono-and sesquiterpene hydrocarbons along with the oxygenated derivatives of monoterpenes. This preponderance of terpenes agreed with the published data from the Nigerian A. nilotica [24] and other Acacia species such as A. albida [24], A. tortilis [25], A. mearnsii [26], and A. cyanophylla [27]. The chemical characterization of the bark EO revealed that the diterpenes are the major constituents with a relative concentration of 54.60%, comprising of diterpene hydrocarbons (53.54%) and traces of oxygenated diterpene (1.06%). With a few exceptions, such as Euphorbia mauritanica L. [33], Lactuca serriola L. [34], and others [35], the phenomenon of diterpene abundance was rare in plants EOs. Four diterpene hydrocarbons were assigned as overall identified diterpenoids, including stachene (48.34%), trachyloban (2.25%), kaur-16-ene (1.94%), and isokaurene (1.01%). Stachene was reported as a major diterpene hydrocarbon in a number of plants, including Chamaecyparis pisifera (Siebold and Zucc.) Endl., Chamaecyparis obtuse (Siebold and Zucc.) Endl., and Thuja orientalis L. [36].
Monoterpenes were found to be the major characterized components of A. nilotica fruits EO, with a relative concentration of 63.32% of mainly oxygenated compounds (55.51%) and hydrocarbons (7.81%). The abundance of the monoterpenoids in the fruits EO agreed with the previous data from Nigerian A. nilotica and other Acacia species [24,26,27,33]. Z-Anethole (22.87%), 2-caren-10-al (3.51%), and borneol (2.31%) were characterized as the main identified oxygenated monoterpenes, while γ-terpinene (7.35%) was the major monoterpene hydrocarbons while 1,8-cineol (3.52%) was assigned. All these compounds were reported in the EOs of Acacia species as traces and/or totally absent [19,24,26,33]. However, the bark EO was found to have traces of monoterpenes in oxygenated forms.
The bark and fruits EOs were found to include the sesquiterpenes with high relative concentrations (39.81% and 34.91%, respectively). The analysis of the bark EO revealed the presence of a low concentration of sesquiterpene hydrocarbons (8.27%) and major oxygenated sesquiterpenes (31.54%). The α-muurolene (2.42%) and bicyclogermacrene (2.41%) were identified as the major sesquiterpene hydrocarbons, while caryophyllene oxide (19.11%) and spathulenol (4.74%) are the main oxygenated sesquiterpenes of the bark EO. Furthermore, the sesquiterpene hydrocarbons were the only detected sesquiterpenes in the EO of the fruits with a relative concentration of 34.91% and a complete absence of oxygenated compounds. From all the assigned sesquiterpene hydrocarbons in fruits EO, trans-caryophyllene (36.95%), α-elemene (4.69%), α-humulene (4.05%), and β-elemene (3.72%) represented the main constituents. The majority of the sesquiterpenes were previously reported from EOs of the Nigerian A. nilotica [24], and other Acacia plants as A. tortilis [25]. Previous studies of EOs derived from Acacia species revealed the prevalence of muurolene, caryophyllene oxide, and spathulenol, where they are determined in considerable concentrations [24][25][26].
Ultimately, the other non-terpenoid compounds are represented as traces of the total mass of the bark EO (3.49%), while they are absent in fruit EO. In bark EO, n-nonacosane (2.51%) and n-dotriacontane (0.98%) were identified. Hydrocarbons are widely identified compounds in Acacia plants [25,26]. As observed in these findings, the net results revealed a significant variation in the quantity and quality as well as the chemical components of the two plant parts and also in Nigerian A. nilotica [24]. This phenomenon of variation might be attributed to the variation of the plant organ, genotypes, humidity, climate, weather, and environmental conditions [41,44].

In Vitro Antiviral Activity
The antiviral activities of EOs derived from the bark and fruits of A. nilotica were in vitro screened against HAV, HSV1, and HSV2, and the maximum non-toxic concentration (MNTC) used for the screening was determined. The results showed that MNTC of EOs of bark and fruits was 500 and 1000 µg/mL, respectively (Table 2). Using the MTT assay, the percentages of antiviral effects of bark and fruits EOs were determined by comparing the viability of cells treated by the virus only with the viability of cells treated by the virus and MNTC of the samples. The results revealed that bark EO has moderate antiviral effects against HAV with an effect of 47.26% alongside a selectivity index (SI) of 2.3, concerning amantadine as a reference drug with SI at 51.62. In addition, the bark EO demonstrated moderate antiviral abilities against HSV1 with an effect of 35.98% and SI of 1.6 compared with acyclovir as a positive control at SI > 387.63. However, this EO exhibited a very weak antiviral effect against HSV2 with a 9.07% comparison with acyclovir with SI at 12.24.
On the other side, the EO derived from the fruits showed weak potentialities against the three tested viruses. This oil showed weak anti-HAV with an effect % of 9.42% and SI of 3.8 compared with amantadine at SI of 51.62. Likewise, this oil also exhibited weak anti-HSV1 with effect% at 14.26% and SI of 5.7 compared with acyclovir as a standard drug. In addition, the results revealed that fruits EO has the lowest activity against HSV2 with an effect% of 3.99% as well as SI at 1.6. The chemical components of the EOs and especially the main compounds played as the main mediators in their pharmaceutical and biological activities [8]. The present results revealed that the bark EO has moderate to weak antiviral potentialities, especially against the two viruses, HAV and HSV1. The GC-MS profiling of this oil revealed the abundance of terpenes, mainly di-and sesqui-terpenes, that were documented to have antiviral effects against several viruses [45,46]. Diterpenes, as the main components of the bark EO, were reported to display significant antiviral effects against HAV and HSV [45,[47][48][49].
Plants' EOs enriched with diterpenes were stated to demonstrate significant antiviral agents such as the different extracts and EOs of Croton lechleri Müll.Arg. [2]. Mechanistically, the diterpenes could inhibit the viral replication process by blocking DNA polymerase activity, as reported with dolastane compounds [49]. Several diterpene skeletons, including kaurenes, have been shown to have significant anti-HSV, anti-HAV, and other antiviral properties [45,49]. In addition, the sesquiterpenes acted in a significant role as inhibitors of several viruses. Astani and his colleagues discovered that caryophyllene oxide has significant anti-HSV1 activity [47]. Many EOs derived from plants, such as Eryngium alpinum L., Eryngium amethystinum L. [50], Melaleuca ericifolia Sm., Melaleuc leucadendron (L.) L., Melaleuc armillaris (Sol. ex Gaertn.) Sm. and Melaleuc styphelioides Sm. [51], have been shown to have antiviral activities due to the majority of caryophyllene oxide and spathulenol. The two sesquiterpenoids, α-muurolene and bicyclogermacrene, were also reported as major components of antiviral active EOs derived from Glechon spathulata Benth. and Glechon marifolia Benth. [52].
Previous studies hypothesized that the anti-HSV effectiveness greatly depends on the binding affinity of EO components to the surface of viruses, and this affinity might be affected by the polarity of compounds [15]. The anti-HSV effects of thymol-related monoterpenoids were reported to decrease with decreasing polarity [53].
Consequently, the present findings were in full agreement with the previous reports. The singular and/or synergetic effects of the major components, along with the minors of the bark EO, might be the main reasons for this anti-HAV and HSV1 efficiency. These results supported the ability of these compounds to cause abnormalities in the HSV membrane protein functions and structures that consequently decrease the penetration or binding of the virus into the cells [54].
On the other hand, the fruits' EO antiviral effects were in concordance with Astani et al. [47] who deduced the anti-HSV with β-caryophyllene and Z-anethole as singular compounds. The increase in the HSV inhibitory effects of the EOs was basically ascribed to the presence of the polyhydroxylated constituents [54]. With the exception of the mono-oxygenated monoterpene, Z-anethole, the other oxygenated monoterpenes were identified as traces that decrease the contribution and effect of each singular compound.

Molecular Docking Studies
In our attempt to study the mechanistic action of each EO, the molecular operating environment (MOE) of the main three compounds in each EO toward the two proteins, HAV 3C protease and HSV TK, were evaluated. The results of in silico studies of the major compounds identified in both oils could rationalize the in vitro antiviral activities noticed by them. As shown in Table 3, caryophyllene oxide and spathulenol, the main constituents of oil of the bark, showed the best affinities (∆G = −5.62, −5.33 and −6.90, −6.76 kcal/mol) for 3C protease and TK, respectively (Figures 2 and 3). Additionally, trans-caryophyllene, as a major constituent of fruit EO, revealed promising binding capabilities against both proteins (∆G = −5.31, −6.58 kcal/mol, respectively) (Figures 2 and 3). compounds identified in both oils could rationalize the in vitro antiviral activities noticed by them. As shown in Table 3, caryophyllene oxide and spathulenol, the main constituents of oil of the bark, showed the best affinities (ΔG = −5.62, −5.33 and −6.90, −6.76 kcal/mol) for 3C protease and TK, respectively (Figures 2 and 3). Additionally, trans-caryophyllene as a major constituent of fruit EO, revealed promising binding capabilities against both proteins (ΔG = −5.31, −6.58 kcal/mol, respectively) (Figures 2 and 3).    Previous studies have shown that the sesquiterpene caryophyllene and its oxides possess suitable antiviral activities with a high selectivity index [47]. The higher activity of the bark oil than the fruit oil might be attributed to the synergistic effect of caryophyllene oxide and spathulenol, as well as the other components. These data go in the same line with the in vitro results in which the bark EO is more active than the fruits EO.

Molecular Dynamic and System Stability
A molecular dynamic simulation was carried out to predict the inhibition performance and interaction of the caryophyllene oxide with the catalytic active site of both 3C protease of HAV and thymidine kinase of HSV targets as well as their stability throughout the simulation [56,57]. The validation of system stability is essential to trace disrupted motions and avoid artifacts that may develop during the simulation. The recorded average RMSD values for all frames of Apo-HAV and Complex-HAV were 1.53 Å and 1.37 Å, respectively. In addition, average RMSD values of 1.51 Å and 1.59 Å were observed for Apo-HSV and Complex-HSV, respectively ( Figure 4A,C). In general, These results revealed that the caryophyllene oxide-bound protein complex system acquired a relatively more stable conformation than the other studied systems.  Previous studies have shown that the sesquiterpene caryophyllene and its oxides possess suitable antiviral activities with a high selectivity index [47]. The higher activity of the bark oil than the fruit oil might be attributed to the synergistic effect of caryophyllene oxide and spathulenol, as well as the other components. These data go in the same line with the in vitro results in which the bark EO is more active than the fruits EO.

Molecular Dynamic and System Stability
A molecular dynamic simulation was carried out to predict the inhibition performance and interaction of the caryophyllene oxide with the catalytic active site of both 3C protease of HAV and thymidine kinase of HSV targets as well as their stability throughout the simulation [56,57]. The validation of system stability is essential to trace disrupted motions and avoid artifacts that may develop during the simulation. The recorded average RMSD values for all frames of Apo-HAV and Complex-HAV were 1.53 Å and 1.37 Å, respectively. In addition, average RMSD values of 1.51 Å and 1.59 Å were observed for Apo-HSV and Complex-HSV, respectively ( Figure 4A,C). In general, These results revealed that the caryophyllene oxide-bound protein complex system acquired a relatively more stable conformation than the other studied systems.  During MD simulation, assessing protein structural flexibility upon ligand binding is critical for examining residue behavior and its connection with the ligand. The protein residue fluctuations were evaluated using the Root-Mean-Square Fluctuation (RMSF) algorithm to evaluate the effect of inhibitor binding toward the respective targets over 20 ns simulations. The computed average RMSF values were 1.21 and 1.00 for Apo-HAV, and caryophyllene oxide-bound to protein systems, respectively, while 1.17 Å and 1.02 Å were recorded for Apo-HSV and Complex-HSV, respectively. Figure 4B,D show the overall residue fluctuations of individual systems. These values reveal that the caryophyllene oxide-bound protein complex system has a lower residue fluctuation than the other systems.

Binding Interaction Mechanism Based on Binding Free Energy Calculation
A popular method for determining the free binding energies of small molecules to biological macromolecules is the molecular mechanics energy technique (MM/GBSA), which combines the generalized born and surface area continuum solvation, and it may be more trustworthy than docking scores. The MM-GBSA program in AMBER18 was used to calculate the binding free energies by extracting snapshots from the trajectories of the systems. As shown in Table 1, all the reported calculated energy components (except ∆G solv ) gave negative values, indicating favorable interactions. The results indicated that binding free energy (∆G bind ) values −19.35 and −32.04 Kcal/mol were obtained for the interactions of caryophyllene oxide with 3C protease of HAV and thymidine kinase of HSV receptors, respectively.
A close look at the individual contributions of energy reveals that the more positive van der Waals energy components drive caryophyllene oxide interactions with 3C protease of HAV and thymidine kinase of HSV enzyme, resulting in the observed binding free energies (Table 4).

Plant Materials and EOs Extraction
Acacia nilotica bark and fruits were collected in August 2020 from Naje Algamal, Sohag Governorate, Egypt (26 • 33 43.2 N 31 • 40 39.7 E). The specimen is authenticated and deposited in the herbarium of the National Research Center with specimen codes NRC-20xYD-AN-20-609.
The collected samples were dried at air temperature (25 • C ± 3) for 7 days and then crushed into a fine powder. The EOs of the air-dried powdered bark and fruits of A. nilotica (300 g, each) were subjected to hydro-distillation over a Clevenger apparatus for three hours. The oily layers were separated via n-hexane and then dried using Na 2 SO 4 anhydrous (0.5 g). Three EO samples from each part were obtained by repeating the extraction process three times with the same sequence. All oil samples were saved in glass vials at 4 • C until further GC-MS and biological analysis.

Gas Chromatography-Mass Spectrometry (GC-MS) Analysis
The GC-MS technique was used for the analysis of the extracted EO samples, as reported in published studies [33,43]. The analysis was performed at the National Research Center using the GC-MS apparatus joined with the TRACE GC Ultra-Gas Chromatograph and quadrupole mass unit, model Thermo-Scientific™ EC, Waltham, MA, USA. The used GC-MS column dimension and film thicknesses were 0.32 mm × 30 m and 0.25 µm. The transporter gas (He) was used with 1 to 10 as the split ratio and 1.0 mL/min as the flow rate. As usual, the temperature was regulated as follows: 60 • C/1 min and elevated to 240 • C during 4 • C/min. The dilution of all oil samples was performed in 1 µL of n-hexane by the ratio of 1:10 (v/v) and injected where both of injector and detector were adjusted at 210 • C. The electron ionization (EI) at 70 eV and m/z 40-450 as spectral range were used for performing the mass spectral data of oil constituents. Finally, the authentication and identification of the chemical constituents occurred depending upon the Automated Mass spectral Deconvolution and Identification (AMDIS) software, (version: 1.0.0.13), NIST library database, Wiley spectral library collection, and retention indices relative to n-alkanes (C 8 -C 22 ).  5, 125, 250, 500, and 1000 µg/mL) of EOs of bark and fruits of A. nilotica in triplicates for another 48 h. Cells were checked for any physical signs of toxicity, e.g., partial or complete loss of the monolayer, rounding, shrinkage, or cell granulation. Then, the cell monolayers were washed 3 times with 1 L sterile phosphate buffer saline (PBS) and treated with 20 µL MTT solution in PBS (5 mg/mL) for each well. After shaking, the plate was incubated for 4 h in 5% CO 2 at 37 • C. The formed formazan crystals in each well were dissolved with 200 µL of 0.04 M HCl in isopropanol. Optical density was determined at 560 nm after removing the background at 620 nm, where it should be directly correlated with cell quantity. The maximum non-toxic concentration (MNTC) of each extract was determined and was used for further biological studies. The MNTC was calculated from the plot of toxicity percent against sample concentration. The DMSO was tested as a control, and it did not show any activity.

Antiviral Effect Percent Determination
All viruses, HSV1, HSV2, and HAV, were obtained from the Microbiology Department, Faculty of Medicine for girls, Al-Azhar University. Vero-E6 cells were distributed in 96well tissue culture plates (10 4 /well) and incubated at a humidified 37 • C incubator under 5% CO 2 conditions. After 24 h, the cells were washed with 1 L of PBS. Equal volumes (1:1 v/v) of MNTC of tested samples and the virus (HSV1, HSV2, and HAV) suspension were incubated for one hour. In triplicate, 100 L of viral or sample suspension was propagated on cells and incubated for 24 h at 37 • C in 5% CO 2 . A total of 20 µL MTT solutions in PBS (5 mg/mL) was added to each well and incubated for 1-5 h to allow MTT to be metabolized. As previously mentioned, the resultant was treated, and the optical density was determined. The antiviral effect percentage was calculated by dividing the viability of cells treated by the virus only/the viability of cells treated by the virus and the sample. Acyclovir was used as a reference drug for the viruses (HSV1 and HSV2), while amantadine was used as a positive control against the HAV virus. The calculation of the SI (selectivity index) values of the tested EO samples was performed from the equation: SI = CC 50 /IC 50 ; where: CC 50 : 50% cytotoxic concentration, and IC 50 : 50% effective concentration [6,11].

Molecular Docking Studies
Molecular Operating Environment (MOE Vs. 2015) docking studies were performed on the catalytic domains of HAV 3C protease (PDB ID: 1QA7) [58] and HSV thymidine kinase (PDB ID: 1KI3) [58] and thymidine kinase of HSV (PDB ID: 1KI3) [59]. The crystal structures of both proteins were retrieved from the protein data bank and processed as previously described [60], as well as the database file (mdb) of the major identified compounds in both EOs, stachene, caryophyllene oxide, spathulenol, trans-caryophyllene, Z-anethole, and γ-terpinene. The docking process was validated by re-docking the co-crystallized ligands in the binding site, which revealed their binding with crucial sub-pockets (Val 28 and Cys 172 and Gln A125, Arg A176, and Glu A83) in both proteins, respectively, at acceptable RMSD (Figures 2 and 3). After processing the study in triplicate (as shown in Table S1), the results of the docking process were presented as the ∆G (kcal/mol) with RMSD values ≤ 2 Å (Table 2). In addition, the interactions of the lowest energy pose with the binding pocket were two and three-dimensionally presented in (Figures 2 and 3).

System Preparation
The crystal structures of the 3c proteinase from the hepatitis A virus receptor and the thymidine kinase from herpes simplex virus type I were retrieved from the protein data bank with codes 1QA7 [58] and 1KI3 [59], respectively. This structure was then prepared for molecular dynamics (MD) studies using UCSF Chimera [61]. Using PROPKA, pH was fixed and optimized to 7.5 (3). Caryophyllene oxide was drawn using ChemBioDraw Ultra 12.1. Altogether, all four prepared systems were subjected to 20 ns MD simulations as described in the simulation section.

Molecular Dynamic (MD) Simulations
The integration of molecular dynamic (MD) simulations into biological systems studies enables exploring the physical motion of atoms and molecules that cannot be easily accessed by any other means. The insight extracted from performing this simulation provides an intricate perspective into the biological systems' dynamical evolution, such as conformational changes and molecule association [61]. The MD simulations of all systems were performed using the GPU version of the PMEMD engine present in the AMBER 18 package [62].
The partial atomic charge of each compound was calculated with ANTECHAMBER's General Amber Force Field (GAFF) technique [63]. The Leap module of the AMBER 18 package implicitly solvated each system within an orthorhombic box of TIP3P water molecules within 10 Å of any box edge. The Leap module was used to neutralize each system by incorporating Na + and Cl − counter ions. A 2000-step initial minimization of each system was carried out in the presence of a 500 kcal/mol applied restraint potential, followed by a 1000-step full minimization using the conjugate gradient algorithm without restraints.
During the MD simulation, each system was gradually heated from 0 K to 300 K over 500 ps, ensuring that all systems had the same amount of atoms and volume. The system's solutes were subjected to a 10 kcal/mol potential harmonic constraint and a 1 ps collision frequency. Following that, each system was heated and equilibrated for 500 ps at a constant temperature of 300 K. To simulate an isobaric-isothermal (NPT) ensemble, the number of atoms and pressure within each system for each production simulation were kept constant, with the system's pressure maintained at 1 bar using the Berendsen barostat [64].
For 20 ns, each system was MD simulated. The SHAKE method was used to constrain the hydrogen bond atoms in each simulation. Each simulation used a 2 fs step size and integrated an SPFP precision model. An isobaric-isothermal ensemble (NPT) with randomized seeding, constant pressure of 1 bar, a pressure-coupling constant of 2 ps, a temperature of 300 K, and a Langevin thermostat with a collision frequency of 1 ps was used in the simulations.

Post-MD Analysis
After saving the trajectories obtained by MD simulations every 1 ps, the trajectories were analyzed using the AMBER18 suite's CPPTRAJ [65] module. The Origin data analysis program and Chimera were used to create all graphs and visualizations.

Thermodynamic Calculation
The Poisson-Boltzmann or generalized born and surface area continuum solvation (MM/PBSA and MM/GBSA) approach has been found to be useful in the estimation of ligand-binding affinities [66,67]. The protein-ligand complex molecular simulations used by MM/GBSA and MM/PBSA compute rigorous statistical-mechanical binding free energy within a defined force field. Binding free energy averaged over 500 snapshots extracted from the entire 50 ns trajectory. The estimation of the change in binding free energy (∆G) for each molecular species (complex, ligand, and receptor) can be represented as follows [68]: G sol = G GB + G SA (4) The terms Egas, Eint, Eele, and Evdw symbolize gas-phase energy, internal energy, Coulomb energy, and van der Waals energy. The Egas was directly assessed from the FF14SB force field terms. Solvation free energy (Gsol) was evaluated from the energy involvement of the polar states (GGB) and non-polar states (G). The non-polar solvation free energy (GSA) was determined from the solvent accessible surface area (SASA) [69,70] using a water probe radius of 1.4 Å. In contrast, solving the GB equation assessed the polar solvation (GGB) contribution. Items S and T symbolize the total entropy of the solute and temperature, respectively.

Statistical Analysis
The presented results were obtained from the mean values of three experimental processes independently performed. CC 50 and IC 50 values were determined from the curve of the dose response along with the regression analysis of the triplicates of the values. A one-way ANOVA followed by multiple Tukey's tests was applied for the comparison of the EOs with p < 0.05 significance in the antiviral assay.

Conclusions
The analysis of the EOs derived from the bark and fruits of Acacia nilotica revealed high relative concentrations of terpenoids in both oils. Stachene, caryophyllene oxide, spathulenol, trans-caryophyllene, Z-anethole, and γ-terpinene represented the main constituents. The EO of the bark exhibited moderate anti-HAV and anti-HSV1, while fruit EO showed weak effects against HAV, HSV1, and HSV2. Caryophyllene oxide and spathulenol exhibited the best affinities against the 3C protease and TK proteins. The molecular dynamics simulation proved the significant van der Waals energy of caryophyllene oxide with 3C protease of HAV and thymidine kinase of HSV enzyme. The present findings revealed the effects of the main constituents of A. nilotica EO. However, in vivo studies should be evaluated for these EOs and/or their major compounds, either in combination or singular, to determine the actual action mechanisms and safety.
Supplementary Materials: The following supporting information can be downloaded at: https:// www.mdpi.com/article/10.3390/plants11212889/s1, Table S1: Three results of docking simulations of the major compounds identified in EOs of both bark and fruits of Acacia nilotica.