Targeting Natural Plant Metabolites for Hunting SARS-CoV-2 Omicron BA.1 Variant Inhibitors: Extraction, Molecular Docking, Molecular Dynamics, and Physicochemical Properties Study

(1) Background: SARS-CoV-2 Omicron BA.1 is the most common variation found in most countries and is responsible for 99% of cases in the United States. To overcome this challenge, there is an urgent need to discover effective inhibitors to prevent the emerging BA.1 variant. Natural products, particularly flavonoids, have had widespread success in reducing COVID-19 prevalence. (2) Methods: In the ongoing study, fifteen compounds were annotated from Echium angustifolium and peach (Prunus persica), which were computationally analyzed using various in silico techniques. Molecular docking calculations were performed for the identified phytochemicals to investigate their efficacy. Molecular dynamics (MD) simulations over 200 ns followed by molecular mechanics Poisson–Boltzmann surface area calculations (MM/PBSA) were performed to estimate the binding energy. Bioactivity was also calculated for the best components in terms of drug likeness and drug score. (3) Results: The data obtained from the molecular docking study demonstrated that five compounds exhibited remarkable potency, with docking scores greater than −9.0 kcal/mol. Among them, compounds 1, 2 and 4 showed higher stability within the active site of Omicron BA.1, with ΔGbinding values of −49.02, −48.07, and −67.47 KJ/mol, respectively. These findings imply that the discovered phytoconstituents are promising in the search for anti-Omicron BA.1 drugs and should be investigated in future in vitro and in vivo research.


Introduction
The World Health Organization (WHO) was informed on 31 December 2019 of a cluster of cases of pneumonia of unknown etiology first detected in Wuhan city, Hubei province, China [1]. A new coronavirus (SARS-CoV-2) with high transmissibility was A molecular docking and molecular dynamics investigation was carried out in the current research to study the efficacy of the identified phytochemicals and to comprehend the nature of the interaction with omicron. In view of achieving a better understanding of the potency of the investigated inhibitors, a drug-likeness analysis was carried out to assess the bioeffects and bioavailability of the best Cmds identified from the molecular docking and dynamics study. The current findings will serve as the foundation for future in vivo and in vitro investigations to ensure the efficacy of the developed anti-Omicron BA.1 drugs.

Plant Material
The aerial parts of Echium angustifolium Mill. were collected from Mersa Matruh (Agiba region), Egypt in April 2020., as well as peach plant aerial parts were collected in March 2020 from Upper Egypt. The plant samples were authenticated by staff members of Plant Taxonomy, Faculty of Science, Assiut University, Egypt. Their voucher specimens (T-21 and V-21, respectively) were deposited.

Extraction and Isolation
Air-dried aerial parts of E. angustifolium were ground to yield 1 kg of plant powder. This powder was extracted at room temperature by maceration with 80% EtOH (3 × 4 L, every 48 h). The obtained extracts were combined and vacuum-concentrated at 45 • C, yielding 73 g of the total extract. This dried extract was then fractionated by suspending in 70% MeOH (400 mL) in a separating funnel and partitioning with n-hexane, yielding 50 g of defatted fraction [29]. This fraction was then fractionated on polyamide 6S column chromatography using the eluent of H 2 O-EtOH (100:0 → 0:100). All obtained fractions were screened using silica gel 60 GF254 NP-TLC plates (Merck, Darmstadt, Germany), yielding a promising sub-fraction. Furthermore, LC-ESI-MS/MS was used to identify the chemical constituents (compound 1, 13, and 15) in the Ea-DfD fraction using Acquity UPLC™ (Waters, Milford, MA, USA) in both ionization modes [36].

Protein Preparation
All in silico calculations were performed using the three-dimensional crystal structure of the SARS-CoV-2 Omicron Variant BA.1 receptor-binding domain (RBD) (PDB ID: 7U0N, resolution: 2.61 Å). The crystal structure of the wild-type SARS-CoV-2 RBD was also investigated (PDB ID: 6M0J, resolution: 2.45 Å). Ions, water molecules, and hetero-atoms were removed from the downloaded viral protein. The H++ web service (http://biophysics.cs.vt.edu/H++, last accessed on 14 May 2022)was employed to study the protonation states of protein amino acids [37,38]. As well, all of the missing hydrogen atoms were added. Physiological parameters such as pH, salinity, internal dielectric, and external dielectric were set to 7.0, 0.15, 10, and 80, respectively, for the H++ computations.

Inhibitor Preparation
Chem3D Pro 12.0 software (version 12.0.2, Cambridge, MA, USA) was exploited in order to sketch and analyze the chemical structures of the fifteen extracted Cmds. The MM2 force field was used to minimize energy in all of the Cmds tested. Such an energy minimization procedure is required prior to molecular docking research to limit the influence of any potential unfavorable torsion angles, bond angles, bond lengths, or undesired non-bonded interactions [39].

Molecular Docking
Molecular docking is widely regarded as the most effective approach in computational drug development for determining the efficacy of drugs under investigation [40]. The binding affinities of these Cmds were investigated using AutoDock Vina (California, CA, USA) [41]. In this investigation, all parameters were left at their normal settings, with the exception of the exhaustiveness parameter, which was set to 200. The BA.1-RBD residues were encased in a docking grid box with XYZ dimensions of 25 × 25 × 25 (Å). Furthermore, the grid spacing was adjusted to 1.0 Å. The resulting nine docked inhibitor poses were examined, and the best one was chosen. The protein-ligand interactions were visualized using BIOVIA Discovery Studio (Waltham, MA, USA) [42]. Table 1 illustrates the mutated residues in BA.1 and the original residues of the wild-type SARS-CoV-2.

Molecular Dynamics Simulations
The YASARA Structure (version 22.5.22, YASARA Biosciences GmbH, Vienna, Austria) protocol was used for molecular dynamics (MD) simulations [43] The main objective of employing MD simulation is to generate deep insight into the protein-ligand complexes' stability. The AMBER14 force field (San Francisco, CA, USA) was the force field used in the MD simulations. Firstly, hydrogen atoms, atomic partial charges, and force field parameters were well studied as the primary structure preparation step. Following system preparation, 5000 step energy minimizations were performed using a combination of steepest and conjugate gradient algorithms to clear any inappropriate geometries or steric clashes. The MD simulations were performed on amino acid residues at the default physiological pH setting (7.4). Water molecules were successfully presented into the system at constant temperature and pressure. To maintain the systems' neutral state, counter ions (Na+ or Cl) at a concentration of 0.9% were added. The Berendsen barostat technique was employed to keep the pressure at 1 atm [44]. The particle-mesh Ewald (PME) method was used to calculate the long-range coulomb forces [45,46]. For non-bonded interactions, the cut-off radius was set to 8 Å. The Langevin thermostat method [47] was used to keep the temperature at 298.15 K through the heating step of the docked protein-ligand complexes from 0 K to 300 K over 10 ps. The periodic boundary conditions were also considered. In each case, the cubic simulation cell was designed to be 20 Å larger than the examined protein-ligand complexes. A consistent simulation speed was maintained for intramolecular processes with a multiple time step of 1.25 fs. The SHAKE method was used to constrain all intermolecular bonds, including hydrogen bonds, at a 2 fs integration step [48]. As a pre-final step, the systems were sufficiently equilibrated for 1 ns. Finally, production processes were completed across simulation periods of 50, 100, 150, and 200 ns. Snapshots of the simulation trajectory were taken every 100 ps, as defined by the root-meansquare deviations (RMSDs) of the solutes from the original structure. The simulation steps were carried out with the help of a pre-installed macro (md_run.mcr) from the YASARA package. The Poisson-Boltzmann method is used in the YASARA program [49]. The surfcost parameter was set at 0.35 in order to calculate the entropic cost of exposing a Å2 to the solvent. The inhibitor's binding energy was calculated using the AMBER14 force field. The binding energies calculated from PBS were divided by a factor of 20 to maintain compatibility with the empirically determined values [50]. Various types of post-MD analyses, such as the root-mean-square deviation (RMSD), radius of gyration (Rg), root-mean-square fluctuation (RMSF), and solvent-accessible surface area (SASA), were all used to study the behavior of the best identified Cmds over the 200 ns MD simulations.

Binding Energy Calculations
The molecular mechanics Poisson-Boltzmann surface area (MM/PBSA) technique was used to calculate the binding free energies of the examined Cmds against BA.1 [51]. The equations shown below were used to calculate the MM/PBSA binding free energy.
The complex, protein, and ligand binding energies are described by ∆G C , ∆G P , and ∆G L , respectively. Furthermore, ∆G Sol , ∆E MM and −T∆S stand for solvation Gibbs energy, gas-phase molecular mechanics change, and conformational entropy, respectively. The term ∆E MM can be calculated by adding the van der Waals and electrostatic interactions. The term ∆G Sol is easily defined as the sum of the polar and non-polar solvation values. The term −T∆S denotes the entropic contribution

Assessment of Drug-Relevant Properties
For the sake of determining physicochemical properties, the Osiris property explorer (https://www.organic-chemistry.org/prog/peo/, last accessed on 21 May 2022) was used. To characterize the extent to which the discovered Cmds were toxic, tumorigenic, mutagenic, irritating, and reproductive qualities were expected. Partition coefficient between n-octanol and water (clogP), water solubility (logS), molecular weight (Mwt), and topological polar surface area (TPSA) were also investigated. All the above parameters were combined to produce a single value known as the drug score (DS), representing the total drug potential. DS can be estimated using the following equation: where (S i ) is: S i = 1 1 + e ap i +b S i represents the contributions estimated directly from cLogP (coctanol/cwater), molecular weight (Mwt), and drug-likeness (pi) via the second equation, which expresses a spline curve. Parameters a and b are (1, −5) for cLogP, (1,5) for logS, (0.012, −6) for Mwt, and (1, 0) for drug-likeness. The ti parameter represents the contributions taken from the four major toxicity risk types (mutagenic, tumorigenic, irritant, and reproductive). The ti values are 1.0, 0.8, and 0.6 for no risk, medium risk, and high risk, respectively.

Identification of Phytoconstituents
Three active components; 1, 13, and 15 were identified from the aerial parts of E. angustifolium. Compound (1) was tentatively identified as echiumin A by comparing its LC-ESI-MS/MS results in the negative ion mode with the previous data [52]. According to our findings, Echiumin A (1) detected a molecular ion peak at m/z 1013. Besides; Twelve active components; 2-12 and 14 were identified from the aerial parts of peach based on literature data;

Molecular Docking
A total of fifteen plant substances were subjected to molecular docking calculations using AutoDock Vina to find out the best active Cmds. The names of the identified Cmds and their docking scores are listed in Tables S1 and S2 and Figure S1. The docking data showed that five Cmds had a docking score of ≥−9.0 kcal/mol kcal/mol against Omicron BA.1. Ligand binding is the fundamental step in the inhibition of protein activity. Accordingly, a comprehensive understanding of the interactions between small Cmds and proteins can serve as the basis for a rational method for drug development.
In Table 2, the docking scores, 2D chemical structures, and binding characteristics for the five most powerful phytochemicals with Omicron BA.1 are all summarized.
For comparison, the best-ranked five compounds were also docked to the identified RBD residues of the wild-type SARS-CoV-2. Figure 2 shows the results of docking analyses in terms of binding mode and binding affinity (kcal/mol).
Contrary to Cmd-1, two additional types of interactions can be noticed, Pi-Alkyl interaction with TYR449 (4.60 Å) and Pi-Pi T-shaped interaction with TYR505 (4.92 Å). It also formed four conventional hydrogen bonds: ARG403 ( For comparison, the best-ranked five compounds were also docked to the identified RBD residues of the wild-type SARS-CoV-2. Figure 2 shows the results of docking analyses in terms of binding mode and binding affinity (kcal/mol).

Figure 2. Depiction of binding modes and docking score for the best-investigated phytochemical
Cmds with the wild-type RBD of SARS-CoV-2.

Molecular Dynamics (MD) Simulations
Molecular dynamics (MD) simulations investigate the stability of receptor-inhibitor complexes, structural details, conformational elasticities and the reliability of receptorinhibitor affinities. Accordingly, MD simulations over 200 ns followed by binding energy estimates for the top five components complexed with the Omicron BA. 1 were performed. The binding free energies were calculated using the MM/PBSA technique. The corresponding binding energies were estimated and shown in Figure 3. Figure 3 shows that in the first 50 ns MD simulations, Cmd-2 and Cmd-4 had binding energies of more than −100.0 KJ/mol. Furthermore, values of binding energies were found to be over −60.0 KJ/mol for Cmd-1 and Cmd-5 and over −50.0 KJ/mol for Cmd-3. In order to increase the reliability of the obtained result, MD simulations were extended to 100 ns and the relevant binding energies were determined. What is noticeable about the estimated binding energies for five Cmds over 100 ns is that Cmd-1 and Cmd-2 still have binding energies of over −60.0 and −100.0 KJ/mol, respectively.
Moreover, computed binding energies over 100 ns for Cmd-3, Cmd-4, and Cmd-5 were slightly lowered by −3.7, −8.0, and −6.6 KJ/mol, respectively. Longer MD simulations of 150 ns and subsequently 200 ns were accomplished for those five potent Cmds in complex with Omicron BA.1 to obtain more reliable results, and the related binding affinities were calculated.
With respect to data from 200 ns MD simulation, Cmd-1 and Cmd-2 obtained binding energies of over −45.0 kJ/mol, whereas Cmd-4 had the highest stability compared to others, with a binding energy value of over −65.0 KJ/mol.

Molecular Dynamics (MD) Simulations
Molecular dynamics (MD) simulations investigate the stability of receptor-inhibitor complexes, structural details, conformational elasticities and the reliability of receptor-inhibitor affinities. Accordingly, MD simulations over 200 ns followed by binding energy estimates for the top five components complexed with the Omicron BA. 1 were performed. The binding free energies were calculated using the MM/PBSA technique. The corresponding binding energies were estimated and shown in Figure 3.  Figure 3 shows that in the first 50 ns MD simulations, Cmd-2 and Cmd-4 had binding energies of more than −100.0 KJ/mol. Furthermore, values of binding energies were found to be over −60.0 KJ/mol for Cmd-1 and Cmd-5 and over −50.0 KJ/mol for Cmd-3. In order to increase the reliability of the obtained result, MD simulations were extended to 100 ns and the relevant binding energies were determined. What is noticeable about the estimated binding energies for five Cmds over 100 ns is that Cmd-1 and Cmd-2 still have binding energies of over −60.0 and −100.0 KJ/mol, respectively.
Moreover, computed binding energies over 100 ns for Cmd-3, Cmd-4, and Cmd-5 were slightly lowered by −3.7, −8.0, and −6.6 KJ/mol, respectively. Longer MD simulations of 150 ns and subsequently 200 ns were accomplished for those five potent Cmds in complex with Omicron BA.1 to obtain more reliable results, and the related binding affinities were calculated.
With respect to data from 200 ns MD simulation, Cmd-1 and Cmd-2 obtained binding energies of over −45.0 kJ/mol, whereas Cmd-4 had the highest stability compared to others, with a binding energy value of over −65.0 KJ/mol.

Post-MD Analyses
To verify the consistency and behavior of the discovered Cmds, their complexes with Omicron BA.1 were structurally and energetically analyzed over 200 ns MD simulations and compared to those of the apoprotein structure. The structural stability of the investigated complexes was assessed by measuring root-mean-square deviation (RMSD), rootmean-square fluctuation (RMSF), the radius of gyration (Rg), and solvent-accessible surface area (SASA).

Post-MD Analyses
To verify the consistency and behavior of the discovered Cmds, their complexes with Omicron BA.1 were structurally and energetically analyzed over 200 ns MD simulations and compared to those of the apoprotein structure. The structural stability of the investigated complexes was assessed by measuring root-mean-square deviation (RMSD), root-meansquare fluctuation (RMSF), the radius of gyration (Rg), and solvent-accessible surface area (SASA).

Root-Mean-Square Deviation (RMSD)
As previously mentioned, the primary goal of the MD simulations is to examine the positional and conformational changes of ligands when they bind to the binding pocket, which provides insight into the binding stability. In order to analyze the structural stability of the best five plant substances in a complex with Omicron BA.1, the root-mean-square deviation (RMSD) of the complete complex backbone atoms was evaluated for comparison purposes. RMSD values indicate the stability of the structures generated during the simulation trajectories [62,63]. RMSD values for the studied phytochemicals in the complex with Omicron BA.1 are plotted in Figure 4. Moreover, variations in the RMSD values for the studied ligands are also plotted in Figure 5.
Lower values of RMSD demonstrate an insignificant difference between the starting and final structures over the MD course.
Throughout the 200 ns MD simulations, the estimated RMSD values for Cmd-2 were found to be 2.81 Å with a difference of 0.51 Å from the apoprotein structure (2.30 Å). As well, Cmd-3, Cmd-4, and Cmd-5 obtained an average RMSD of 3.22, 3.87, and 3.96 Å, respectively. In contrast, Cmd-1 showed higher values of RMSD (avg. = 5.53 Å), which reflects the occurrence of some positional or conformational changes. As shown in Figure 6, the ligand movement for Cmd-2 and Cmd-3 within the binding pocket reflects the RMSD value, which is the lowest compared to the rest of the compounds. which provides insight into the binding stability. In order to analyze the structural stability of the best five plant substances in a complex with Omicron BA.1, the root-meansquare deviation (RMSD) of the complete complex backbone atoms was evaluated for comparison purposes. RMSD values indicate the stability of the structures generated during the simulation trajectories [62,63]. RMSD values for the studied phytochemicals in the complex with Omicron BA.1 are plotted in Figure 4. Moreover, variations in the RMSD values for the studied ligands are also plotted in Figure 5.   which provides insight into the binding stability. In order to analyze the structural stability of the best five plant substances in a complex with Omicron BA.1, the root-meansquare deviation (RMSD) of the complete complex backbone atoms was evaluated for comparison purposes. RMSD values indicate the stability of the structures generated during the simulation trajectories [62,63]. RMSD values for the studied phytochemicals in the complex with Omicron BA.1 are plotted in Figure 4. Moreover, variations in the RMSD values for the studied ligands are also plotted in Figure 5.

Root-Mean-Square Fluctuation (RMSF)
A root-mean-square fluctuation (RMSF) analysis was executed in order to study the flexibility of the Omicron BA.1 residues over the 200 ns MD simulations. This is a numerical measurement comparable to RMSD, except that instead of expressing positional changes across complete structures over time, RMSF calculates individual residue flexibility or how much a specific residue moves (fluctuates) throughout a simulation [64]. RMSF per residue is often displayed vs. residue number and can indicate which amino acids in a protein contribute the most to molecular movement structurally. A high RMSF value indicates greater flexibility, whereas a low RMSF value indicates limited mobility. The RMSF graph for the five studied components and the apoprotein structure is illustrated in Figure 7.
Curr. Issues Mol. Biol. 2022, 44 5040 Lower values of RMSD demonstrate an insignificant difference between the starting and final structures over the MD course.
Throughout the 200 ns MD simulations, the estimated RMSD values for Cmd-2 were found to be 2.81 Å with a difference of 0.51 Å from the apoprotein structure (2.30 Å). As well, Cmd-3, Cmd-4, and Cmd-5 obtained an average RMSD of 3.22, 3.87, and 3.96 Å, respectively. In contrast, Cmd-1 showed higher values of RMSD (avg. = 5.53 Å), which reflects the occurrence of some positional or conformational changes. As shown in Figure  6, the ligand movement for Cmd-2 and Cmd-3 within the binding pocket reflects the RMSD value, which is the lowest compared to the rest of the compounds.

Root-Mean-Square Fluctuation (RMSF)
A root-mean-square fluctuation (RMSF) analysis was executed in order to study the flexibility of the Omicron BA.1 residues over the 200 ns MD simulations. This is a numerical measurement comparable to RMSD, except that instead of expressing positional changes across complete structures over time, RMSF calculates individual residue flexibility or how much a specific residue moves (fluctuates) throughout a simulation [64].
RMSF per residue is often displayed vs. residue number and can indicate which amino acids in a protein contribute the most to molecular movement structurally. A high RMSF value indicates greater flexibility, whereas a low RMSF value indicates limited mobility. The RMSF graph for the five studied components and the apoprotein structure is illustrated in Figure 7.  The averages values of RMSF for apoprotein structure, Cmd-1, Cmd-2, Cmd-3, Cmd-4, and Cmd-5 were 1.50, 1.86, 1.53, 1.54, 1.55, and 1.57 (Å), respectively.

The Radius of Gyration (Rg)
The radius of gyration (Rg) is a parameter that shows how the protein's structural compactness changes during simulation [65]. An Rg analysis reflects the protein's flexibility and stability [65]. Protein has the lowest Rg values and, as a result, the tightest packing and the most stability.

The Radius of Gyration (Rg)
The radius of gyration (Rg) is a parameter that shows how the protein's structural compactness changes during simulation [65]. An Rg analysis reflects the protein's flexibility and stability [65]. Protein has the lowest Rg values and, as a result, the tightest packing and the most stability.
During the 200 ns simulation time, all Cmds demonstrated appropriate Rg behavior, as shown in Figure 8.   Figure 5 further shows that Cmd-3 and Cmd-4 exhibited constant Rg behavior during the simulation, indicating that the protein structure is highly compact.

The Solvent-Accessible Surface area Analysis (SASA)
Solvent-accessible surface area analysis (SASA) was computed to obtain a deeper insight into interactions between the complex and the solvent throughout the 200 ns MD simulation. Figure 9 depicts the SASA vs. simulation time curve for apoprotein and five complexes.
3.4.4. The Solvent-Accessible Surface area Analysis (SASA) Solvent-accessible surface area analysis (SASA) was computed to obtain a deeper insight into interactions between the complex and the solvent throughout the 200 ns MD simulation. Figure 9 depicts the SASA vs. simulation time curve for apoprotein and five complexes. The average SASA values for Cmd-1, Cmd-2, Cmd-3, Cmd-4, and Cmd-5 were 10,317.33, 10,412.28, 10,669.97, 10,437.59 and 10,500.047 (Å 2 ), respectively.

Drug-Relevant Properties
Evaluation of molecular features should be carried out through the inspection of many parameters. For instance, drug solubility (logS) is a crucial element that influences drug transport from the site of administration into the bloodstream. It is also well recognized that low medication solubility can result in poor absorption. The logS value should be greater than −4.0. The second factor is the partition coefficient between n-octanol and water, cLogP, which is used to determine the molecule's hydrophilicity. For the molecule to be well absorbed, values of cLogP should be in the range of −4.0 to 5.0. Figure 10 illustrates the estimated values of logS and cLogP for the five inspected Cmds. All Cmds are within the allowed range, implying good absorption.

Drug-Relevant Properties
Evaluation of molecular features should be carried out through the inspection of many parameters. For instance, drug solubility (logS) is a crucial element that influences drug transport from the site of administration into the bloodstream. It is also well recognized that low medication solubility can result in poor absorption. The logS value should be greater than −4.0. The second factor is the partition coefficient between n-octanol and water, cLogP, which is used to determine the molecule's hydrophilicity. For the molecule to be well absorbed, values of cLogP should be in the range of −4.0 to 5.0. Figure 10 illustrates the estimated values of logS and cLogP for the five inspected Cmds. All Cmds are within the allowed range, implying good absorption.  Numerical values of the drug-likeness parameter should be in the positive range. A positive value indicates that the molecule primarily comprises fragments that are commonly found in commercial medications. The drug score combines all computed values of druglikeness, cLogP, logS, Mwt, and toxicity risks into a single useful value that can be used to assess the compound's overall potential to qualify for drug approval. Figure 11 depicts the values of the estimated drug-likeness and drug score for the five phytoconstituents.
Both Cmd-1 and Cmd-3, which had negative drug-likeness values, had the lowest drug scores when compared to others. Strikingly, Cmd-2 obtained the highest values of drug-likeness and drug score when compared to other Cmds.
Another parameter to consider is toxicity, which indicates whether the molecule has a low risk (G), medium risk (Y), or high risk (R). The toxicity prediction is executed based on a precomputed collection of structural fragments that cause toxicity alerts if detected in the studied structure. Table 3 displays data from the toxicity evaluation for the five plant substances.
As shown in Table 3, Cmd-2 showed low toxicity risks towards the four inspected toxicity parameters. Cmd-1 and Cmd-3 exhibited high reproductive and tumorigenic toxicity risks, respectively. Excluding mutagenicity, Cmd-5 showed low tumorigenic, irritant, and reproductive toxicity risks. Cmd-4 has a medium reproductive side effect. Numerical values of the drug-likeness parameter should be in the positive range. A positive value indicates that the molecule primarily comprises fragments that are commonly found in commercial medications. The drug score combines all computed values of drug-likeness, cLogP, logS, Mwt, and toxicity risks into a single useful value that can be used to assess the compound's overall potential to qualify for drug approval. Figure 11 depicts the values of the estimated drug-likeness and drug score for the five phytoconstituents. Figure 11. Calculated drug-likeness and drug scores of the best studied Cmds.
Both Cmd-1 and Cmd-3, which had negative drug-likeness values, had the lowest drug scores when compared to others. Strikingly, Cmd-2 obtained the highest values of drug-likeness and drug score when compared to other Cmds. Figure 11. Calculated drug-likeness and drug scores of the best studied Cmds. Table 3. Toxicity risks of the five studied Cmds.

Toxicity Risk
Mu Tu Ir Re

Conclusions
In this work, fifteen compounds from Echium angustifolium and peach were successfully identified and screened for possible inhibition of the SARS-CoV-2 Omicron BA.1 variant using various in silico methods. According to a molecular docking study, five phytochemicals showed promising inhibition of Omicron BA.1, with a docking score greater than −9.0 kcal/mol. The dynamic behavior of the five best docked protein-ligand complexes was studied on a time scale of 200 ns. This revealed different effects of the ligands on docking to Omicron BA.1, which were assessed by their respective RMSD, RMSF, Rg and SASA. Moreover, the analysis of molecular properties such as solubility (logS), hydrophilicity (cLogP), toxicity, drug-likeness, and drug score revealed that the studied components had a high chance of becoming drug-like compounds. The results of docking, dynamics and physicochemical property evaluation suggest that the annotated phytoconstituents from Echium angustifolium and peach are promising drug candidates. Further developments, especially in vitro and in vivo studies, are an important step to test these compounds as effective anti-Omicron BA.1 drugs.