Phytochemical Analysis, Antioxidant, and Antimicrobial Activities of Ducrosia flabellifolia: A Combined Experimental and Computational Approaches

Ducrosia flabellifolia Boiss. is a rare desert plant known to be a promising source of bioactive compounds. In this paper, we report for the first time the phytochemical composition and biological activities of D. flabellifolia hydroalcoholic extract by using liquid chromatography–electrospray tandem mass spectrometry (ESI-MS/MS) technique. The results obtained showed the richness of the tested extract in phenols, tannins, and flavonoids. Twenty-three phytoconstituents were identified, represented mainly by chlorogenic acid, followed by ferulic acid, caffeic acid, and sinapic acid. The tested hydroalcoholic extract was able to inhibit the growth of all tested bacteria and yeast on agar Petri dishes at 3 mg/disc with mean growth inhibition zone ranging from 8.00 ± 0.00 mm for Enterococcus cloacae (E. cloacae) to 36.33 ± 0.58 mm for Staphylococcus epidermidis. Minimal inhibitory concentration ranged from 12.5 mg/mL to 200 mg/mL and the hydroalcoholic extract from D. flabellifolia exhibited a bacteriostatic and fungistatic character. In addition, D. flabellifolia hydroalcoholic extract possessed a good ability to scavenge different free radicals as compared to standard molecules. Molecular docking studies on the identified phyto-compounds in bacterial, fungal, and human peroxiredoxin 5 receptors were performed to corroborate the in vitro results, which revealed good binding profiles on the examined protein targets. A standard atomistic 100 ns dynamic simulation investigation was used to further evaluate the interaction stability of the promising phytocompounds, and the results showed conformational stability in the binding cavity. The obtained results highlighted the medicinal use of D. flabellifolia as source of bioactive compounds, as antioxidant, antibacterial, and antifungal agent.


Plant Material Sampling
In this proposal, D. flabellifolia (Al-Hazaa; Figure 1) plant species were collected from the Hail region (Al-Mu'ayqilat, 27 • 16 41.9" N, 41 • 22 48.0" E) in October 2019. A voucher specimen (AN04) was deposited at the herbarium in the Department of Biology (College of Science, University of Hail, Hail, Kingdom of Saudi Arabia). For the experiment, 4 g of air-dried aerial parts were macerated in 100 mL of methanol-80% at room temperature for 72 h. The filtrate was recuperated by lyophilization, and the yield (expressed in percentage) was calculated using the following Equation (1): where W1 is the weight of extract after the evaporation of solvent and W2 is the dry weight of the sample. The yield of extraction was about 21.56 ± 1.78%.

Study of the Phytochemical Composition
An Agilent Technologies 1260 Infinity liquid chromatography system (Santa Clara, CA, USA) hyphenated to a 6420 Triple Quad mass spectrometer was used for quantitative analyses. Chromatographic separation was carried out on a Poroshell 120 EC-C18 (100 mm × 4.6 mm I.D., 2.7 µm) column (Santa Clara, CA, USA). The previously validated method was used for the analysis of phenolic compounds by LC-ESI-MS/MS [27]. The mobile phase was made up of solvent A (0.1%, v/v formic acid solution) and solvent B (methanol). The gradient profile was set as follows: 0.00 min 2% B eluent, 3.00 min 2% B eluent, 6.00 min 25% B eluent, 10.00 min 50% B eluent, 14.00 min 95% B eluent, 17.00 min 95% B, and 17.50 min 2% B eluent. The column temperature was maintained at 25 • C. The flow rate was 0.4 mL min −1 and the injection volume was 2.0 µL. The tandem mass spectrometer was interfaced with the LC system via an ESI source. The electrospray source of the MS was operated in negative and positive multiple reaction monitoring (MRM) mode and the interface conditions were as follows: capillary voltage of −3.5 kV, gas temperature of 300 • C, and gas flow of 11 L min −1 . The nebulizer pressure was 40 psi. MRM transitions, the optimum collision energies, and retention times for each species are indicated in Supplementary Material S1. In addition, representative LC-ESI-MS/MS chromatograms of phenolic compounds are shown in Supplementary Material S2. Calibration curves and sensitivity properties of the method are also shown in Supplementary Material S3.
In negative and positive multiple reaction monitoring (MRM) mode, the peaks of the analytes were identified by comparing the retention time, together with monitoring ion pairs in an authentic standard solution.
2.4. Screening of the Biological Activities 2.4.1. Antimicrobial Activities The antimicrobial activity of Al-Haza extracts was tested against twelve ESKAPE bacterial strains including Enterococcus faecium, three Staphylococcus species (S. aureus, S. epidermidis, and S. hominis), Klebsiella pneumoniae, Acinetobacter baumannii, Pseudomonas aeruginosa, two Enterobacter species (E. cloacae and E. faecalis), and Escherichia coli. Four types of Candida species were also tested (C. albicans ATCC 20402, C. tropicalis ATCC 1362. C. guillermondii ATCC 6260, and C. utilis ATCC 9255). The effect of the hydroalcoholic extract from D. flabellifolia was estimated using disc diffusion assay [28] by measuring the diameter of the growth inhibition zone tested on a Mueller Hinton agar plate for bacterial strains and Sabouraud dextrose agar for yeast. Sterile discs were impregnated with three different concentrations of the tested extract (1 mg/disc, 1.5 mg/disc, and 3 mg/disc). After incubation, the mean diameter of the growth inhibition zone (mGIZ) was calculated, and the scheme proposed by Parveen et al. [29] was used to interpret the obtained results (no activity: mGIZ = 0; low activity: mGIZ = 1-6 mm; moderate activity: mGIZ = 7-10 mm; high activity: mGIZ = 11-15 mm; and very high activity: GIZ = 16-20 mm). Ampicillin and Amphotericin B were used as controls.
The microdilution assay was used for the determination of MIC (minimal inhibitory concentration) and MBC/MFC values (minimal bactericidal/fungicidal concentration) as previously described by Khalfaoui et al. [30]. Bacterial and fungal cultures were inoculated into the wells of 96-well microtiter plates in the presence of D. flabellifolia methanolic extract at final concentrations varying from 0.039 mg/mL to 100 mg/mL. To interpret the character of the tested extract, we used the ratios (MBC/MIC ratio and MFC/MIC ratio) described by Gatsing et al. [31] and Moroh et al. [32].

Phytochemistry and Antioxidant Activities Screening
Total phenolic content expressed as milligram (mg) of gallic acid per gram of plant extract (mg GAE/g extract) was estimated by using the Folin-Ciocalteu method as previously described by Kumar et al. [33]. Total flavonoids expressed as mg of quercetin equivalents per gram of plant extract (mg QE/g extract) were determined using the AlCl 3 method developed by Benariba et al. [34]. In addition, acidified vanillin method previously described by Broadhurst and Jones [35] was used to estimate the total condensed tannins (expressed mg tannic acid equivalent per gram of plant extract (mg TAE/g extract)).
The ability of the Al-Haza extract against DPPH-H was determined following the same method as Mseddi al. [36]. The method of Koleva et al. [37] for β-Carotene bleaching test and Oyaizu [38] method for the determination of reducing power were used.

Molecular Docking
The 3D structures of LC-MS-detected phytochemicals were retrieved in structural data format (sdf) from the PubChem database. The LigPrep module was used to prepare the entire set of phytochemicals, which involved the addition of hydrogen atoms and suitable charges, as well as the correcting of the valences and the protonation and tautomeric states at pH 7.2 ± 2.0 of the molecules using the Epik tool [39,40]. The X-ray crystal structure of S. aureus type IIA topoisomerase, tyrosyl-tRNA synthetases (TyrRS) of S. aureus, Sap1 of C. albicans, and human peroxiredoxin 5 receptor were retrieved from the PDB server with accession codes of 2XCT, 1JIJ, 2QZW, and 1HD2, respectively, as the receptor for the molecular docking study. The Schrodinger Protein Preparation Wizard has been used for protein preparation and energy minimization, in which crystallographic water molecules are removed, then missing hydrogen and/or side chain atoms are added, and the correct charges and protonation states are given to protein residues at pH 7.0 [41,42]. The protein structure was energy minimized using the OPLS3 force field to alleviate steric conflicts in the protein structure. Following that, the prepared protein was considered for grid generation utilizing the "Receptor Grid Generation" panel where active sites in targeted proteins were defined for grid generation by selecting cognate ligands [43]. Glide's standard precision (SP) mode was used to execute all docking calculations for the prepared ligand molecules, with the default parameters.

Molecular Dynamics (MD) Simulation
MD simulation is regarded as the most important approach for comprehending the nature of biological macromolecules' underlying structure and function. The MD simulations were conducted utilizing the Schrödinger MD simulation program Desmond, which helps us to understand how a ligand-protein complex binds in simulated physiological conditions [44]. MD simulation analysis was conducted for top score phyto-compound Hyperoside in complex with 1JIJ protein. The simulation system was built through the system builder, with the boundaries set by an orthorhombic shape with a diameter of 10 Å × 10 Å × 10 Å and filled with SPC water molecules.
To neutralize the system's charges, sodium (+27) and chloride (−12) ions were supplied as counter ions and salt concentration was maintained at 0.15 M. Following the setup of the system builder, the protein-ligand complex system was minimized using the steepest descent method, which was followed by LBFGS algorithms with a maximum of 2000 iterations [45,46]. Further, the system was equilibrated using the NPT ensemble with a Nose-Hoover chain thermostat to evenly distribute the ions and solvent throughout the protein-ligand complex. During the simulation, the temperature was maintained at 300 K. Furthermore, isotropic position scaling was used to control 1bar barostat pressure [47][48][49][50]. A 100 ns simulation was run, with a total of 1000 frames stored in the simulation system for subsequent analysis using the Desmond program's "Simulation Interactions Diagram" module.

Statistical Analysis
All experiments were performed in triplicate and average values were calculated using the SPSS 25.0 statistical package for Windows. Duncan's multiple-range tests for means with a 95% confidence interval (p ≤ 0.05) were used to calculate the differences in means.

Phytochemical Composition of D. Flabellifolia Hydroalcoholic Extract
Using ESI-MS/MS technique, 23 phytocompounds were identified in the hydroalcoholic extract from D. flabellifolia aerial parts ( The chemical structures of the main identified phytoconstituents are summarized in Figure 2.  Table 1. Table 2 summarizes the results of the quantification of tannins, phenols, and flavonoids. In fact, the hydroalcoholic extract was dominated by phenolic compounds (46.684 ± 0.757 mg GAE/g extract), followed by tannins (6.204 ± 0.401 mg TAE/g extract), and flavonoids (1.816 ± 0.133 mg QE/g extract). The results obtained showed that the tested methanolic extract from D. flabellifolia was able to scavenge the DPPH radical with a low IC 50 value (0.014 ± 0.045 mg/mL) as compared to BHT and AA (0.023 ± 3 × 10 −4 mg/mL and 0.022 ± 53 × 10 -4 mg/mL, respectively). In addition, the ABTS + radical was scavenged with an IC 50 value of about 0.102 ± 0.024 mg/mL of the tested extract. For beta-carotene bleaching inhibiting property, the IC 50 value was estimated at 7.80 ± 0.919 mg/mL, significantly different from BHT and AA (p < 0.05). Table 2. Antioxidant activities of D. flabellifolia methanol/water extract compared to known drugs.

Test Systems D. flabellifolia (Methanol-80%) Butylated Hydroxytoluene Ascorbic Acid
The letters (a,b) indicate a significant difference between the different antioxidant methods according to the Duncan test (p < 0.05).

Antimicrobial Activities Screening
The results of the antimicrobial activity of the hydro-methanolic extract from D. flabellifolia aerial parts showed that the mean diameter of the growth inhibition zone (mGIZ) increased in a concentration-dependent manner (Figure 3). At 3 mg/disc, the mGIZ ranged from 8.00 ± 0.00 mm (E. cloacae) to 36.33 ± 0.58 mm (S. epidermidis). At the same concentration (3 mg/disc), C. utilis ATCC 9255 was the most sensitive strain with mGIZ of about 13.67 ± 0.58 mm. All bacterial strains seem to be more sensitive to the tested extract6 as compared to yeast strains, and the reference drug (ampicillin). Using the scheme proposed by Parveen et al. [29], the tested extract at 3 mg/mL showed high to very high activity against all tested bacterial and fungal strains with mGIZ from 11 to 20 mm. Using the microdilution technique, the minimal inhibitory concentrations values (MICs) ranged from 12.5 to 25 mg/mL for bacterial strains and the minimal bactericidal concentration values (MBCs) ranged from 100 to 200 mg/mL. Using the scheme proposed by Gatsing et al. [31] and Moroh et al. [32], the tested extract showed bacteriostatic activity against almost all tested ESKAPE microorganisms (MBC/MIC ratio higher than 4 and lower than 16), with the exception of E. faecium and E. cloacae (MBC/MIC equal to 4). All these data are summarized in Table 3. Similarly, the tested D. flabellifolia hydroalcoholic extract was able to inhibit the growth of Candida strains in liquid media at 25 mg/mL (Table 4), while a high concentration of the tested extract was needed to completely kill them (from 100 to 200 mg/mL), with MFC/MIC ratios varying from 4 (C. albicans ATCC 20402; fungicidal action) to 8 (Fungistatic action against C. utilis ATCC 9255, C. tropicalis ATCC 1362, and C. guillermondii ATCC 6260). The molecular docking study of identified phyto-compounds was carried out against bacterial fungal and human peroxiredoxin 5 targets. The binding affinities of phytocompounds are shown in the Supplementary Material (S1). The results of the molecular docking study revealed that phyto-compound luteolin 7-glucoside showed the highest binding affinity toward topoisomerase IIA (2XCT) with a docking score −12.562 Kcal/mol. The phyto-compounds apigenin 7-glucoside, quercetin, (−)-epicatechin, eriodictyol, and hyperoside from D. flabellifolia also showed significant binding affinities against topoisomerase IIA (2XCT) with docking scores of −12.514, −10.06, −9.836, −9.582, and −9.159 Kcal/mol, respectively. The best binding pose of luteolin 7-glucoside inside the active site of topoisomerase IIA and the 2D and 3D representations of the interactions with the amino acids inside the binding pocket are presented in Figure 4.
Luteolin 7-glucoside binding interaction reveals that it interacted with Arg1122 and manganese (II) ion via π-cationic and metal coordinates, respectively. Among the assessed phyto-compounds, five compounds, namely, (−)-epicatechin, eriodictyol, 2,5-dihydroxybenzoic acid, chlorogenic acid, luteolin 7-glucoside, and vanillin showed proper and promising antioxidant activity by proper recognition at the binding active site of human peroxiredoxin protein with docking scores of −5.791, −5.255, −5.249, −5.11, −5.06, and −5.02 Kcal/mol, respectively. It was reported that hydrogen bonds and hydrophobic interactions established with the surrounding amino acids are applied to anticipate the binding modes of the novel compounds compared to the antioxidant benzoic acid cognate ligand, at the 1HD2 complex. This active pocket contained conserved amino acid residues such as Thr44, Gly46, Cys47, and Arg127, which play critical roles in docked compound recognition via hydrogen bonding and hydrophobic interactions. The binding interaction shows that two terminal hydroxyl groups connected with asp145 and the147 in the (−)-epicatechin-1HD2 complex, whereas the chromane ring hydroxyl group showed triple hydrogen bonding with Cys47, Thr44, and Arg127 amino acids, suggesting that the identified phyto-compound (−)-epicatechin has antioxidant activity ( Figure 5).  Similarly, in C. albicans Sap1 protein (2QZW), hyperoside showed the highest binding affinity with a docking score of −6.055 Kcal/mol, and phyto-compounds eriodictyol, (−)-epicatechin, and syringic acid showed potential binding affinities against C. albicans Sap1 (2QZW) protein with docking scores of −5.814, −5.461, and −5.042 Kcal/mol, respectively. The residues of Thr6, Asp17, Lys26, Gly103, and Ala104, from Sap1 of C. albicans formed seven significant interactions with the phyto-compound hyperoside with a standard hydrogen binding pattern as shown in Figure 6. Among assessed phyto-compounds, hyperoside also showed a promising docking score (−8.852 Kcal/mol) in the microbial target TyrRS from S. aureus (1JIJ) in which it interacted with active residues, namely, Asp195, Gly193, Gly38, and Asp177 (Figure 7).  Table 5 shows the promising phyto-compounds hydrogen bonding profile in the active site of type IIA topoisomerase, TyrRS, Sap1, and human peroxiredoxin 5 proteins.

MD Simulation
Using molecular dynamics simulation, the docked complex of the phyto-compound hyperoside at the binding site of TyrRS of S. aureus was simulated under biological environments to investigate complex stability and protein flexibility. Hyperoside, a phytocompound, has shown a high affinity for microbial targets, hence it was chosen for the MD simulation studies. The MD trajectories were used to determine RMSD values, root-meansquare fluctuation (RMSF) values, and protein-ligand interactions. Various MD trajectory data analyses for the hyperoside-1JIJ complex are shown in Figure 8. The RMSD figure showed a stable ligand-protein complex throughout the simulation time, as shown by RMSD values ranging from 1.6 to 2.8 Å for protein Cα atoms in the complex with hyperoside ( Figure 8A). RMSD changes remain within 3 Å throughout the simulation period, which is perfectly appropriate for small, globular proteins such as TyrRS. In the case of hyperoside RMSD with respect to protein, it was found that the ligand RMSD ranged from 0.8 to 2.6 Å. Except for a minor fluctuation, the RMSD of hyperoside was found to remain steady throughout the simulation. The maximum ligand RMSD was recorded at 59 and 66 ns, when RMSD reached 2.59 and 2.60 Å, respectively. Table 5. Interacting active site residues of receptors with the best phyto-compounds identified in D. flabellifolia hydroalcoholic extract.

Name of Complex
Interacting Residues Furthermore, the flexibility of the protein system was evaluated during the simulation by computing the RMSF of individual protein amino acid residues. The higher-peaking residues correspond to loop regions (highlighted by white shade) identified by MD trajectories or N-and C-terminal areas. Low RMSF values of binding site residues show that hyperoside binding to TyrRS protein is stable. From Figure 8B  The 2D ligand interaction diagram shows the charged negative amino acid Asp establishing a major hydrogen bond with a hydroxyl group of hyperoside. The amino acids His50, Gly38, Val191, and Thr75 interacted with hyperoside for 57%, 92%, 43%, and 57% of simulation time, respectively. Additionally, there are intramolecular hydrogen bonds in hyperoside between the hydrogen atom of hydroxyl and the carbonyl oxygen of chromen-4-one moiety ( Figure 8C). The binding interactions between hyperoside and active site amino acid residues inside the binding pocket of TyrRS were computed and represented in Figure 8D. Most of the important interactions of hyperoside with the TyrRS protein determined with MD simulations are hydrogen bonds, polar (amino acids mediated hydrogen bonding) interactions, and hydrophobic interactions.

Discussion
In the present study, we report for the first time the identification of bioactive molecules by using the ESI-MS/MS technique in the hydroalcoholic extract from D. flabellifolia aerial parts collected from the Saudi Arabian desert. The main identified phytoconstituents were: protocatechuic acid, chlorogenic acid, 2,5-dihydroxybenzoic acid, 4-hydroxybenzoic acid, caffeic acid, syringic acid, vanillin, sinapic acid, p-coumaric acid, ferulic acid, hyperoside, 2-hydroxycinnamic acid, and pinoresinol.
As far as we know, there are few data available in the literature relating to the chemical compounds from D. flabellifolia extracts. In fact, Talib et al. [14] have demonstrated that D. flabellifolia 95% ethanolic extract collected from Wadi Hassan (Jordan) was a rich source of flavonoids and terpenoids qualitatively estimated by the thin layer chromatography technique. The same team has demonstrated the identification of several flavonoids (flavonols and flavones) by using the HPLC-MS/MS technique in D. flabellifolia ethanolic extract including quercetin, fisetin, kaempferol, luteolin, and apigenin.
Our results have shown that D. flabellifolia hydroalcoholic extract was active against ESKAPE pathogens and Candida species with an increase in the mean diameter of the growth inhibition zone depending on the concentration used. The highest sensitivity for all tested microorganisms was obtained at 3 mg/disc of the tested extract. Using the schemes proposed by Gatsing et al. [31] and Moroh et al. [32], the tested extract exhibited bacteriostatic and fungistatic action against almost all tested microorganisms. No study has hitherto described the antimicrobial activity of D. flabellifolia extracts. While testing Jordanian D. flabellifolia essential oil from dried leaves, Al-Shudiefat et al. [10] reported good activity against C. albicans with a MIC value of about 234 µg/mL; meanwhile, for bacterial strains, the MICs values were about 234 µg/mL against S. aureus, 1870 µg/mL against E. coli, and 1872 µg/mL against P. aeruginosa strain. For D. anethifolia, Elsharkawi and colleagues [11] have demonstrated that the ethyl acetate fraction was active against different Gram-positive bacteria including S. aureus ATCC 25923, S. epidermidis ATCC 49461, Bacillus cereus ATCC 10876, and S. aureus clinical isolate tested at 500 mg/mL. The highest diameter of growth inhibition was recorded against S. epidermidis (14.5 ± 0.5 mm) followed by B. subtilis (14.0 ± 0.00 mm), while no activity was recorded against E. coli ATCC 35218, K. pneumoniae ATCC 700603, K. pneumoniae ATCC 27736, P. aeruginosa, and A. baumannii.
We also reported in this study that D. flabellifolia hydroalcoholic extract was dominated by phenolic compounds, followed by tannins and flavonoids. The same extract showed the ability to scavenge different radicals with low IC 50 values (IC 50 DPPH radical = 0.014 ± 0.045 mg/mL; IC 50 ABTS + radicals = 0.102 ± 0.024 mg/mL; and IC 50 for the beta-carotene test = 7.80 ± 0.919 mg/mL). Previously, Mottaghipisheh et al. [13] reported that the IC 50 value of the free radical scavenging activity of D. anethifolia ethanolic extract was estimated at 122.02 ppm and 354.37 ppm for the ethyl acetate extract. In addition, Elsharkawi and colleagues [11] revealed the presence of high reduction capacity (EC 50 : 0.63 ± 0.03 g/L) and ability to scavenge the free radicals of DPPH with an IC 50 of 0.38 ± 0.02 g/L in ethyl acetate fraction of D. anethifolia.
Overall, the good antimicrobial activity against all tested Gram-positive/Gram-negative and Candida species, and the antioxidant activity of the tested hydroalcoholic extract can be attributed to the presence of many molecules with potential biological activities (Table 6). Table 6. Literature survey showing some biological activities of the identified molecules in D. flabellifolia hydroalcoholic extract.
[ Molecular docking studies were performed in order to explore binding affinity and interaction of identified phyto-compounds against bacterial, fungal, and human peroxiredoxin 5 receptor. A docking score was used to determine the binding affinities of the phyto-compounds to target proteins. A lower docking score indicates a higher affinity. Luteolin 7-glucoside connected with Arg1122 and manganese (II) ion in the S. aureus IIA topoisomerase target via π -cationic and metal coordinates, respectively. This interaction is comparable to that of the co-crystallized standard drug ciprofloxacin. According to the study, the binding site cavity's metal interaction had an impact on the substance's physicochemical properties and antibacterial activity [103,104]. Free radical-scavenging activity properties using the 3-D crystallographic peroxiredoxin 5 (1HD2) were carried out to explore the identified phyto-compounds' recognition in the active site as potential antioxidants. With a docking score of -5 kcal/mol, five phyto-compounds, namely (−)-epicatechin, eriodictyol, 2,5-dihydroxybenzoic acid, chlorogenic acid, luteolin 7-glucoside, and vanillin, showed promising antioxidant activity by appropriate orientation at the binding active site of human peroxiredoxin protein. In fungal target C. albicans Sap1 and bacterial TyrRS protein, the phyto-compound hyperoside showed the highest binding affinity. SB-219383 and its analogs are a family of bacterial TyrRS inhibitors that are both potent and selective. These inhibitors' crystal structures have been solved in complex with the TyrRS from S. aureus, the bacterium that causes the majority of hospital-acquired infections, according to crystal structure. SB-219383 and its analogs interacted with charged negative amino acids Asp195 and Asp80, as well as hydrophobic Tyr170, in the active region of TyrRS from S. aureus [105]. Our docking results for hyperoside match those found in the deposited crystal structure, indicating that it has potential antibacterial activity by inhibiting TyrRS S. aureus. Because of its high affinity for microbial targets, hyperoside was selected for MD simulation studies. The RMSD values of the protein Cα atoms were used to calculate the stability of the protein-ligand complex during dynamics analysis. In fact, Cα atoms in proteins' RMSD is a crucial parameter of the MD simulation trajectory that is used to investigate the back-bone deviation of a single frame created in a dynamic environment. The unfolding of the protein molecule results in a high RMSD value, which suggests compactness. The steady fluctuation in the RMSD value across the simulation time indicates that the protein-ligand complex has equilibrated [106][107][108][109][110]. It is anticipated that the lesser RMSD value during the simulation reveals that the protein-ligand complex is more stable. In contrast, the higher RMSD value indicates the protein-ligand complex is less stable [111][112][113]. The overall RMSD revealed that fluctuations in the range of 1.6 Å to 2.8 Å were within the standard range (1-3 Å) of RMSD, indicating that the protein-ligand complex is stable. Throughout the simulation process, the RMSF value shows the mobility and flexibility of each amino acid in the protein. In this plot, the peaks indicate the fluctuation of each amino acid residue over the entire simulation. It implies that higher RMSF values represent higher residue flexibility, whereas lower RMSF values reflect less residue flexibility and better system stability. If the residues in the active site and main chain fluctuated slightly, it indicated that the conformational change was minimal, implying that the reported lead compound was firmly bound within the cavity of the target protein binding pocket [114,115]. In addition, it was observed that the phyto-compound hyperoside interacted with 31 amino acids of the TyrRS protein during the simulation period, which was highlighted by the green vertical bar. All these intercalated residues have low RMSF values, indicating that the hyperoside binding to TyrRS protein is stable.

Conclusions
In summary, we report for the first time in this study the identification of various phenolic compounds in the hydroalcoholic extract from D. flabellifolia aerial parts growing wild in Saudi Arabia. The tested extract possessed good antimicrobial activities and exhibited notable potency in scavenging free radicals. Docking studies on the identified phyto-compounds in the bacterial, fungal, and peroxiredoxin 5 receptors were carried out to confirm the in vitro results, and they exhibited satisfactory binding profiles. In-depth protein-ligand interaction stability in the dynamic state was evaluated using 100 ns MD simulation studies, which revealed a significant binding affinity of the discovered hyperoside towards the TyrRS protein, implying anti-bacterial effectiveness via TyrRS enzyme inhibition. The obtained results highlighted the possible use of this plant species as a source of molecules with therapeutic effects to be used as antimicrobial and antioxidant agents.

Data Availability Statement:
The data presented in this study are available in the article and supplementary materials.

Conflicts of Interest:
The authors declare no conflict of interest.