Jusanin, a New Flavonoid from Artemisia commutata with an In Silico Inhibitory Potential against the SARS-CoV-2 Main Protease

A new flavonoid, Jusanin, (1) has been isolated from the aerial parts of Artemisia commutata. The chemical structure of Jusanin has been elucidated using 1D, 2D NMR, and HR-Ms spectroscopic methods to be 5,2′,4′-trihydroxy-6,7,5′-trimethoxyflavone. Being new in nature, the inhibition potential of 1 has been estimated against SARS-CoV-2 using different in silico techniques. Firstly, molecular similarity and fingerprint studies have been conducted for Jusanin against co-crystallized ligands of eight different SARS-CoV-2 essential proteins. The studies indicated the similarity between 1 and X77, the co-crystallized ligand SARS-CoV-2 main protease (PDB ID: 6W63). To confirm the obtained results, a DFT study was carried out and indicated the similarity of (total energy, HOMO, LUMO, gap energy, and dipole moment) between 1 and X77. Accordingly, molecular docking studies of 1 against the target enzyme have been achieved and showed that 1 bonded correctly in the protein’s active site with a binding energy of −19.54 Kcal/mol. Additionally, in silico ADMET in addition to the toxicity evaluation of Jusanin against seven models have been preceded and indicated the general safety and the likeness of Jusanin to be a drug. Finally, molecular dynamics simulation studies were applied to investigate the dynamic behavior of the Mpro-Jusanin complex and confirmed the correct binding at 100 ns. In addition to 1, three other metabolites have been isolated and identified to be сapillartemisin A (2), methyl-3-[S-hydroxyprenyl]-cumarate (3), and β-sitosterol (4).

Artemisia commutata Besser is a perennial herb found in northern and eastern Kazakhstan. The plant is also found in Europe, Mongolia, Altai, western and eastern Siberia, Primorsky Krai, Amur Region, and Sakhalin in the Far East, as well as in the Volga-Kama region of the European part of Russia [19].
The WHO stated on the COVID-19 dashboard on 27 December 2021, that the global confirmed infections had reached 279,114,972 and the number of deaths passed 5,397,580 [39]. These alarming numbers require massive and continuous work in the field of drug discovery to face the COVID-19 pandemic.
The computer-aided (in silico, computational or cheminformatics) drug design and discovery is a widely employed approach. The ability to determine both the physical and chemical properties of a molecule enabled the identification of molecular similarities and hence the prediction of activity. Additionally, molecular docking and dynamic simulations studies can determine the interaction of a specific molecule with a certain protein successfully [40]. Our team presented several reports in which they utilized computational chemistry to find a cure against COVID-19 [41][42][43][44][45].
Interestingly, the inhibitory effect of flavonoids has been reported before. For instance, rutin was predicted to active against COVID-19 3CL proteinase through exhibiting the best binding affinity among 2030 other natural compounds [46].
We herein report the isolation and structure elucidation of the new flavonoid, Jusanin, (jusan is the original kazak name of wormwood) from the aerial parts of A. commutata. Because Jusanin is a new compound, its potential anti-COVID-19 effect was investigated. Additionally, ADMET, toxicity, and DFT properties of Jusanin have been studied. Finally, molecular dynamics simulation studies confirmed the expected affinity of Jusanin. Besides Jusanin, three known metabolites have been isolated and identified as сapillartemisin A (2), methyl-3-[S-hydroxyprenyl]-cumarate (3), and β-sitosterol (4).

Isolation of Jusanin
To study the composition of A. commutata [47], 0.94 kg of the aboveground parts were collected in the East Kazakhstan region (Western Altai Mountains). Raw materials afforded 20 g of chloroform extract.
The chloroform extract was subjected to a further isolation by silica gel column using heptane-ethyl acetate-methanol as a mobile phase in a manner of increasing polarity. Fraction 5 was further purified utilizing Sephadex LH-20 to afford the crystalline substance (1) a yellow color with a m.p. more than 350 • C. Based on the spectral data of UV, 1 H, and 13 C NMR, with two-dimensional experiments, HSQC, Dept, COSY, and HMBC (Table 1), compound 1 was identified to be the new compound 5,2 ,4 -trihydroxy-6,7,5 -trimethoxyflavone, Jusanin ( Figure 1). The molecular formula C 18

Isolation of Jusanin
To study the composition of A. commutata [47], 0.94 kg of the aboveground parts were collected in the East Kazakhstan region (Western Altai Mountains). Raw materials afforded 20 g of chloroform extract.

Molecular Similarity Study
The basic principle to understand the molecular similarity technique is the wellknown link between the bioactivity of a ligand and interactions of that ligand with a protein target. These interactions are based on physical and chemical properties as the ability to donate or accept hydrogen bond interactions as well as the hydrophobic interactions. In consequence, if two molecules have similar chemical structures, they are expected to have similar steric configuration in addition to similar atoms that act as H-bond donors, acceptors, or hydrophobic centers. Accordingly, these two compounds are predicted to exhibit similar bioactivity too [51].
The chemical structure of Jusanin was compared to the chemical structures of eight co-crystallized ligands of eight essential SARS-CoV-2 proteins (Figure 4). This experiment

Molecular Similarity Study
The basic principle to understand the molecular similarity technique is the well-known link between the bioactivity of a ligand and interactions of that ligand with a protein target. These interactions are based on physical and chemical properties as the ability to donate or accept hydrogen bond interactions as well as the hydrophobic interactions. In consequence, if two molecules have similar chemical structures, they are expected to have similar steric configuration in addition to similar atoms that act as H-bond donors, acceptors, or hydrophobic centers. Accordingly, these two compounds are predicted to exhibit similar bioactivity too [51].
The chemical structure of Jusanin was compared to the chemical structures of eight co-crystallized ligands of eight essential SARS-CoV-2 proteins ( Figure 4). This experiment was designed to check if Jusanin has a similarity with one of these ligands and hence it may show a potential inhibitory impact anti-COVID-19 effect.
weight, and number of hydrogen bond donors, the co-crystalized ligand (X77) of the SARS-Cov-2 main protease (PDB ID: 6W63) (red ball) has good similarity with Jusanin (green ball). The minimum distance between the red and green ball is 0.92. and on the other hand, the co-crystallized ligands of the other SARSCoV-2 proteins (blue balls) did not show similarity with Jusanin (green ball).
As shown in Figure 5 which clarifies the correlation between Alog p, molecular weight, and number of hydrogen bond donors, the co-crystalized ligand (X77) of the SARS-Cov-2 main protease (PDB ID: 6W63) (red ball) has good similarity with Jusanin (green ball). The minimum distance between the red and green ball is 0.92. and on the other hand, the co-crystallized ligands of the other SARSCoV-2 proteins (blue balls) did not show similarity with Jusanin (green ball).
The results revealed a high degree of similarity between the Jusanin and X77 (Table 2 and Figure 5).
In order to verify the obtained molecular similarity results, a fingerprint study was carried out using Discovery Studio. The fingerprint technique calculates the 2D molecular features of two different compounds or more in a binary format. The fingerprint study examined the presence or absence of the coming parameters: charge [59], hybridization [60], H-bond acceptor and donor [61], positive and negative ionizable [62], halogen [63], aromatic [64], and ALogP [65].

Structural Fingerprints Study
In order to verify the obtained molecular similarity results, a fingerpri carried out using Discovery Studio. The fingerprint technique calculates the 2 features of two different compounds or more in a binary format. The finge examined the presence or absence of the coming parameters: charge [59], h [60], H-bond acceptor and donor [61], positive and negative ionizable [62], aromatic [64], and ALogP [65].
In structural keys, the chemical structure of a molecule is encoded into string (0′s and 1′s), each bit of which corresponds to a "pre-defined" struc (e.g., substructure or fragment). If the molecule has a pre-defined feature, th corresponding to this feature is set to 1 (ON). Otherwise, it is set to 0 (OFF) [ Fingerprint approach depends on the Tanimoto coefficient. In this wo

Structural Fingerprints Study
In structural keys, the chemical structure of a molecule is encoded into a binary bit string (0 s and 1 s), each bit of which corresponds to a "pre-defined" structural feature (e.g., substructure or fragment). If the molecule has a pre-defined feature, the bit position corresponding to this feature is set to 1 (ON). Otherwise, it is set to 0 (OFF) [66].
Fingerprint approach depends on the Tanimoto coefficient. In this work, there are three parameters controlling the binary bit string. SA refers to the bits number that computed in the co-crystallized ligands of SARSCoV-2 proteins and Jusanin. SB refers to the bits number that computed in the co-crystallized ligands of SARSCoV-2 proteins but not in Jusanin. SC refers to the bits number that computed in Jusanin but not in the co-crystallized ligands of SARSCoV-2 proteins.
The results of the fingerprint study revealed that X77 has the highest similarity value (0.512) with Jusanin. In addition, X77 showed high SA value with acceptable rage of SB and SC values indicating the presence of a high number of bits in both X77 and Jusanin. These findings verified the significant fingerprint similarity of Jusanin with X77 (Table 3).

Flexible Alignment Studies
Three-dimensional flexible alignment of Jusanin with the co-crystallized ligands (X77) was studied. The result of flexible alignment revealed the generally good overlap of Jusanin with the reference molecule (X77). In addition, Jusanin showed the same spatial orientation of X77. In detail, the 2,3-dimethoxyphenol moiety showed the same orientation of both pyridine and 1H-imidazole moieties of X77. In addition, the 4H-pyran-4-one moiety of Jusanin exhibited close orientation to the phenyl ring of X77. Furthermore, the 2methoxybenzene-1,4-diol moiety of Jusanin showed the same orientation of the cyclohexyl moiety of X77 ( Figure 6).

Flexible Alignment Studies
Three-dimensional flexible alignment of Jusanin with the co-crystallized ligands (X77) was studied. The result of flexible alignment revealed the generally good overlap of Jusanin with the reference molecule (X77). In addition, Jusanin showed the same spatial orientation of X77. In detail, the 2,3-dimethoxyphenol moiety showed the same orientation of both pyridine and 1H-imidazole moieties of X77. In addition, the 4H-pyran-4-one moiety of Jusanin exhibited close orientation to the phenyl ring of X77. Furthermore, the 2-methoxybenzene-1,4-diol moiety of Jusanin showed the same orientation of the cyclohexyl moiety of X77 ( Figure 6).

DFT Studies
The DFT study is an in silico method that calculates the molecular orbital analysis in addition to the molecular electrostatic potential maps of a certain compound based on the examined parameters (Table 4) [67,68]. The DFT study determines the degree of reactivity of a compound. Accordingly, the DFT properties of Jusanin and X77 were studied by Discovery Studio software to disclose the similarity between them from this point. The results were summarized in Table 4, Figures 7 and 8. The functional used in this test was PWC of local density approximation (LDA). In addition, the quality was adjusted to be Coarse which uses DN basis set with SCF density convergence of 1.0e-4 as utilized from Accelrys in the DMol3 module of Materials Studio package.

DFT Studies
The DFT study is an in silico method that calculates the molecular orbital analysis in addition to the molecular electrostatic potential maps of a certain compound based on the examined parameters (Table 4) [67,68]. The DFT study determines the degree of reactivity of a compound. Accordingly, the DFT properties of Jusanin and X77 were studied by Discovery Studio software to disclose the similarity between them from this point. The results were summarized in Table 4, Figures 7 and 8. The functional used in this test was PWC of local density approximation (LDA). In addition, the quality was adjusted to be Coarse which uses DN basis set with SCF density convergence of 1.0e-4 as utilized from Accelrys in the DMol3 module of Materials Studio package.

Molecular Orbital (MO) Analysis
The compounds also exhibited almost equal dipole moment values of 2.395 and 2.906. Additionally, the gap energy of Jusanin (0.083 Ha) was less than that of X77 (0.096 Ha) which denotes a higher degree of reactivity. Consequently, Jusanin may serve as a promising candidate for further studies.
As shown in Figure 7, the analysis of EHOMO and ELUMO, which are corresponding to the highest occupied and lowest unoccupied molecular orbital energies, respectively. EHOMO and ELUMO were related to the chemical reactivity of the investigated compound and its stability. The HOMO spatial (the electron transfer zones) distributions of

Molecular Orbital (MO) Analysis
The compounds also exhibited almost equal dipole moment values of 2.39 Additionally, the gap energy of Jusanin (0.083 Ha) was less than that of X77 which denotes a higher degree of reactivity. Consequently, Jusanin may serve ising candidate for further studies. denotes a higher degree of reactivity. Consequently, Jusanin may serve as a promising candidate for further studies.
As shown in Figure 7, the analysis of EHOMO and ELUMO, which are corresponding to the highest occupied and lowest unoccupied molecular orbital energies, respectively. EHOMO and ELUMO were related to the chemical reactivity of the investigated compound and its stability. The HOMO spatial (the electron transfer zones) distributions of X77 are mainly on the 1H-imidazole and the two amide moieties, while its LUMO spatial (the electron acceptor zones) distributions are on the pyridine and phenyl moieties. For compound Jusanin, the HOMO spatial distributions are mainly distributed on 5-hydroxy-6,7-dimethoxy-4H-chromen-4-one moiety, while its LUMO spatial distributions are located on the 2-methoxybenzene-1,4-diol.

Molecular Electrostatic Potential Maps (MEP)
MEP is an in silico method that calculates partial charges [69], electronegativity [70], intermolecular interaction [71], and chemical reactivity [71] to determine the electrostatic potential of a compound in a 3D form [72]. In MEP study, the electronegative atoms (Hbond acceptors) appear in red. On the other side, the electron-poor atoms (H-bond donors) appear in blue. Additionally, the neutral atoms (hydrophobic interactors) appear in a green to yellow color [73].
The MEP of Jusanin and X77 was illustrated in Figure 8A,B, respectively. The ligand, X77, showed four red patches indicating H-bond acceptors and one blue patch (an H-bond donor). Additionally, the high possibility for hydrophobic interactions was indicated by the presence of yellow patches on the aromatic and aliphatic moieties.
For Jusanin, it showed eight red patches (H-bond acceptors) and two blue patches (H-bond donors). Additionally, there is a yellow patch on the aromatic system indicating the possibility of hydrophobic interaction. These findings validate the ability of Jusanin to interact with M pro like X77.

Docking Studies
To validate the results of the previous experiments, the binding interaction of Jusanin was investigated against the SARS-CoV-2 main protease (PDB ID: 6W63) through docking studies. Compound X77, was used as a reference. The binding free energy (∆G) and the correct binding mode were employed as the base to evaluate the biding. First of all, a validation procedure of the docking process was applied via re-docking of X77 against M pro . The resulting RMSD value was 1.7 • A to indicate the validity of the docking process ( Figure 9). X77 are mainly on the 1H-imidazole and the two amide moieties, while its LUMO spatial (the electron acceptor zones) distributions are on the pyridine and phenyl moieties. For compound Jusanin, the HOMO spatial distributions are mainly distributed on 5-hydroxy-6,7-dimethoxy-4H-chromen-4-one moiety, while its LUMO spatial distributions are located on the 2-methoxybenzene-1,4-diol. MEP is an in silico method that calculates partial charges [69], electronegativity [70], intermolecular interaction [71], and chemical reactivity [71] to determine the electrostatic potential of a compound in a 3D form [72]. In MEP study, the electronegative atoms (Hbond acceptors) appear in red. On the other side, the electron-poor atoms (H-bond donors) appear in blue. Additionally, the neutral atoms (hydrophobic interactors) appear in a green to yellow color [73].
The MEP of Jusanin and X77 was illustrated in Figure 8A,B, respectively. The ligand, X77, showed four red patches indicating H-bond acceptors and one blue patch (an H-bond donor). Additionally, the high possibility for hydrophobic interactions was indicated by the presence of yellow patches on the aromatic and aliphatic moieties.
For Jusanin, it showed eight red patches (H-bond acceptors) and two blue patches (H-bond donors). Additionally, there is a yellow patch on the aromatic system indicating the possibility of hydrophobic interaction. These findings validate the ability of Jusanin to interact with M pro like X77.

Docking Studies
To validate the results of the previous experiments, the binding interaction of Jusanin was investigated against the SARS-CoV-2 main protease (PDB ID: 6W63) through docking studies. Compound X77, was used as a reference. The binding free energy (∆G) and the correct binding mode were employed as the base to evaluate the biding. First of all, a validation procedure of the docking process was applied via re-docking of X77 against M pro . The resulting RMSD value was 1.7 °A to indicate the validity of the docking process ( Figure 9).  The binding of X77 exhibited a binding free energy of −22.81 kcal/mol. The pyridine moiety was buried in the first pocket of the receptor to form one H-bond with His163 and one hydrophobic interaction with Leu141. Furthermore, the 1H-imidazole was oriented into the second pocket and interacted hydrophobically with Cys145 and His41. The tertbutylbenzene moiety occupied the third pocket closely to the amino acids Arg188, Met49, and Leu167. The same moiety reacted with two hydrophobic interactions with Cys145 and His41. The cyclohexyl moiety was buried in the fourth pocket forming a hydrophobic interaction with the amino acid Pro168. The amide moiety formed one H-bond interaction with Glu166 ( Figure 10). The binding of X77 exhibited a binding free energy of −22.81 kcal/mol. The pyridine moiety was buried in the first pocket of the receptor to form one H-bond with His163 and one hydrophobic interaction with Leu141. Furthermore, the 1H-imidazole was oriented into the second pocket and interacted hydrophobically with Cys145 and His41. The tertbutylbenzene moiety occupied the third pocket closely to the amino acids Arg188, Met49, and Leu167. The same moiety reacted with two hydrophobic interactions with Cys145 and His41. The cyclohexyl moiety was buried in the fourth pocket forming a hydrophobic interaction with the amino acid Pro168. The amide moiety formed one H-bond interaction with Glu166 ( Figure 10).  The binding mode of the Jusanin showed a binding free energy of -19.54 kcal/mol ( Table 5). The 2,3-dimethoxyphenol moiety was oriented into the first pocket of the receptor to form one H-bond with Cys145. It showed close contact with Asn142, Glu166, Phe140, His163, and Leu141. In addition, the 4H-pyran-4-one was buried in the second pocket forming one hydrophobic interaction with Cys145. It formed one H-bond with Gly143. The 2-methoxybenzene-1,4-diol moiety occupied the third pocket in close contact with Arg188 and Met49. It formed three hydrophobic bonds with Cys44, Met49, and His41. Additionally, it formed one H-bond with Cys44 ( Figure 11).

In Silico ADMET Analysis
The possibility of a molecule being a drug depends on its activity as well as its bioavailability [74]. In consequence, the in silico ADMET parameters of Jusanin were determined by Discovery Studio software. Simeprevir was used as a reference drug. The results were illustrated in Table 5 and Figure 12.
The results revealed that Jusanin has a very low chance to penetrate BBB indicating its high safety margin against CNS. Moreover, Jusanin exhibited good aqueous solubility and intestinal absorption. Finally, Jusanin was predicted to be CYP2D6 non-inhibitor and can bind plasma protein by more than 90%.

In Silico Toxicity Studies
The proposed toxicity of Jusanin was determined in silico utilizing Discovery Studio software. Seven different toxicity models were applied using Simeprevir as a reference drug. The results were summarized in Table 6. FDA rodent carcinogenicity in male mice expected that Jusanin is non-carcinogenic. Additionally, Jusanin showed high carcinogenic potency TD50 value of 13.5782 mg/kg body weight/day, comparing to Simeprevir (0.2803 mg/kg/day). Additionally, Jusanin showed a high rat maximum tolerated dose value of 0.3502 g/kg compared to Simeprevir (0.0030 g/kg body weight). Jusanin showed a rat oral LD50 value of 0.4621 g/kg which was higher than the reference drug LD50 = 0.2088 g/kg. Regarding the rat chronic LOAEL model, Jusanin exhibited a high value of 0.0624 g/kg, higher than Simeprevir (0.0021 g/kg body weight). Finally, Jusanin was predicted to have mild irritancy against the ocular irritancy model and non-irritant against the skin irritancy model. His41. Additionally, it formed one H-bond with Cys44 ( Figure 11).

In silico ADMET Analysis
The possibility of a molecule being a drug depends on its activity as well as its bioavailability [74]. In consequence, the in silico ADMET parameters of Jusanin were determined by Discovery Studio software. Simeprevir was used as a reference drug. The results were illustrated in Table 5 and Figure 12.
The results revealed that Jusanin has a very low chance to penetrate BBB indicating its high safety margin against CNS. Moreover, Jusanin exhibited good aqueous solubility and intestinal absorption. Finally, Jusanin was predicted to be CYP2D6 non-inhibitor and can bind plasma protein by more than 90%. Figure 12. The expected ADMET study. ADMET_AlogP98: lipid-water partition coefficient; AD-MET_PSA_2D: polar molecular surface area. Two-dimensional polar surface area (PSA_2D) each drug is plotted against their computed atom-type partition coefficient (ALogP98). The area encompassed by the ellipse represents good absorption without any violation of the ADMET properties. Based on Egan et al. [75] absorption model the 95% and 99% confidence limit ellipses corresponding to the blood-brain barrier (BBB) and intestinal absorption models are indicated.

In Silico Toxicity Studies
The proposed toxicity of Jusanin was determined in silico utilizing Discovery Studio software. Seven different toxicity models were applied using Simeprevir as a reference drug. The results were summarized in Table 6.
FDA rodent carcinogenicity in male mice expected that Jusanin is non-carcinogenic. Additionally, Jusanin showed high carcinogenic potency TD50 value of 13.5782 mg/kg body weight/day, comparing to Simeprevir (0.2803 mg/kg/day). Additionally, Jusanin showed a high rat maximum tolerated dose value of 0.3502 g/kg compared to Simeprevir Figure 12. The expected ADMET study. ADMET_AlogP98: lipid-water partition coefficient; AD-MET_PSA_2D: polar molecular surface area. Two-dimensional polar surface area (PSA_2D) each drug is plotted against their computed atom-type partition coefficient (ALogP98). The area encompassed by the ellipse represents good absorption without any violation of the ADMET properties. Based on Egan et al. [75] absorption model the 95% and 99% confidence limit ellipses corresponding to the blood-brain barrier (BBB) and intestinal absorption models are indicated.

Molecular Dynamics (MD) Simulations Studies
MD simulation is a powerful accurate in silico method that can explore the exact binding mode, flexibility, and stability of a specific receptor-ligand complex for a certain time [76]. Molecular dynamics (MD) simulations studies were employed to mimic the dynamic nature of M pro -Jusanin interaction investigating the simulation for binding complex stability at 100 ns.
The dynamic changes of atoms in addition to the conformational variations of atoms backbone in the M pro -Jusanin complex were estimated by RMSD to investigate their stability in both apo-and ligand-bonded states. The M pro -Jusanin exhibited low RMSD values (less than 0.6 nm) with no major fluctuations declaring the great stability of the complex. The M pro -Jusanin complex was stable till 70 ns~. Although the M pro -Jusanin complex showed a minor fluctuation later, it started to be more stable again at 85 ns~( Figure 13A).
The flexibility of M pro was estimated in terms of RMSF to identify the regions of M pro that are being fluctuated through the simulation. Figure 13B shows that the binding of M pro with Jusanin did not change the flexibility of M pro . The compactness of the M pro -Jusanin was represented by the radius of gyration (R g ). The decrease in fluctuation throughout the simulation period indicates the higher compactness of the examined system. The R g of the M pro -Jusanin complex reached a stable conformation with the radius of gyration fluctuating around 2.2 nm ( Figure 13C). The interaction between M pro -Jusanin complex and solvents was investigated by solvent accessible surface area (SASA) over the simulation time. Consequently, SASA of the M pro -Jusanin complex was calculated to investigate the degree of the conformational changes that happened after the interaction. Interestingly, as shown in Figure 13D, M pro exhibeted a decrease in the surface area showing a comparatively lower value of SASA than the starting period. H-bonding between the M pro -Jusanin complex as anessential factor to stabilize the complex observed ( Figure 13E). The highest number of conformations of the M pro formed up to three H-bonds with Jusanin.

General Experimental Section
NMR spectra were carried out on a commercial instrument (Bruker Avance 300 a 600 MHz), chemical shifts () are presented in parts per million (ppm) and re-calculat with respect to tetramethylsilane (TMS) ( 1 H) or carbon signals of deuterium solvents ( 13 C Spin-spin coupling constants (J) are given in hertz (Hz). The refinement of 13 C NMR sp tra signals was carried out using Dept, HSQC, HMBC NMR spectra. Mass spectra we recorded on an HP5989A instrument (CI and EI, ionization energy 70 eV) with Apollo 3 data, and on a Kratos MS50TC instrument for accurate calculations (reaching by elect shock (ESI), common solvent mixture: CH2Cl2 -MeOH  NH4OAc) with MASSLYNX syste data. UV spectra were obtained on a Perkin-Elmer Lambda 20 Spectrometer instrume Melting points were determined on Reichert Thermovar. For column chromatograph silica gel 0.06-0.2 mm (Acros) was used as the stationary phase. Silica gel 32-63 mesh w used for flash column chromatography.

Plant Material
To study the component composition of Artemisia commutata Besser, or wormwo substitute (family Asteraceae), the aboveground part of the plant is collected in the easte Kazakhstan region (Western Altai Mountains) on July, phase of blooming-beginning

General Experimental Section
NMR spectra were carried out on a commercial instrument (Bruker Avance 300 and 600 MHz), chemical shifts (δ) are presented in parts per million (ppm) and re-calculated with respect to tetramethylsilane (TMS) ( 1 H) or carbon signals of deuterium solvents ( 13 C). Spin-spin coupling constants (J) are given in hertz (Hz). The refinement of 13 C NMR spectra signals was carried out using Dept, HSQC, HMBC NMR spectra. Mass spectra were recorded on an HP5989A instrument (CI and EI, ionization energy 70 eV) with Apollo 300 data, and on a Kratos MS50TC instrument for accurate calculations (reaching by electric shock (ESI), common solvent mixture: CH 2 Cl 2 − MeOH + NH 4 OAc) with MASSLYNX system data. UV spectra were obtained on a Perkin-Elmer Lambda 20 Spectrometer instrument. Melting points were determined on Reichert Thermovar. For column chromatography, silica gel 0.06-0.2 mm (Acros) was used as the stationary phase. Silica gel 32-63 mesh was used for flash column chromatography.

Plant Material
To study the component composition of Artemisia commutata Besser, or wormwood substitute (family Asteraceae), the aboveground part of the plant is collected in the eastern Kazakhstan region (Western Altai Mountains) on July, phase of blooming-beginning of flowering.
Species is identified by botanists of the Altai Botanical Garden (Rider city, eastern Kazakhstan). The herbarium sample is stored in the International Scientific Research Holding «Phytochemistry» Fund. The herbarium sample code is 2007.10.02.02.03.

Extraction and Isolation
In total, 1.04 kg of raw material was placed in a round-bottomed flask and filled with chloroform and heated to the boiling point of solvents. This operation was repeated three times. The solvent was evaporated on a rotary evaporator under the vacuum of a water-jet pump to obtain an extract weighing 20 g, which was used for preparative chromatographic separation by column chromatography on silica gel using Heptane-Ethylacetate-Methanol followed by Sephadex LH-20 (Pharmacia Fine Chemicals), Uppsala, Sweden.

Molecular Similarity
The molecular similarity of Jusanin against the co-crystallized ligands of SARS-Cov-2 was carried out and calculated using Discovery Studio 4.0 [77] (See Supplementary Materials).

Fingerprint Study
A fingerprint study of Jusanin against the nine co-crystallized ligands of SARS-Cov-2 was carried out and calculated using Discovery Studio 4.0 (See Supplementary Materials).

DFT
The DFT parameters of Jusanin were calculated using Discovery Studio software [78] (See Supplementary Materials).

Docking Studies
The docking investigation was carried out for Jusanin using MOE2014 software. The results of the docking process were then visualized using Discovery Studio 4.0 software [79][80][81] (See Supplementary Materials).

Molecular Dynamics Simulations
The system was prepared using the web-based CHARMM-GUI [87][88][89] interface with the CHARMM36 force field [90]. All the simulations were carried out using the NAMD 2.13 [91] package. The TIP3P explicit solvation model and CHARMM general force field were used [92,93] (See Supplementary Materials).

Conclusions
A new flavonoid, Jusanin, (1) has been isolated from the aerial parts of Artemisia commutata. Jusanin showed a high structural similarity degree with X77, the co-crystallized ligand of the COVID-19 main protease (PDB ID: 6W63), M pro . This result was indicated by molecular similarity, fingerprint, and DFT studies. The molecular docking of 1 against M pro confirmed the correct binding of 1 inside M pro exhibiting a binding energy of −19.54 Kcal/mol. ADMET and toxicity properties of 1 indicated its likeness to be a drug as well as its general safety. The MD simulation studies at 100 ns confirmed the correct binding of the M pro -Jusanin complex. These interesting results open the door to finding a treatment against COVID-19 after in vitro and in vivo studies.
Supplementary Materials: The following data are available online, 1D, 2D NMR and HR-Ms spectra data of Jusanin, detailed toxicity report of Jusanin, in addition to the utilized in silico methods.