key: cord-0767421-pnqzvqk8 authors: Suleimen, Yerlan M.; Jose, Rani A.; Suleimen, Raigul N.; Arenz, Christoph; Ishmuratova, Margarita Y.; Toppet, Suzanne; Dehaen, Wim; Alsfouk, Bshra A.; Elkaeed, Eslam B.; Eissa, Ibrahim H.; Metwaly, Ahmed M. title: Jusanin, a New Flavonoid from Artemisia commutata with an In Silico Inhibitory Potential against the SARS-CoV-2 Main Protease date: 2022-03-01 journal: Molecules DOI: 10.3390/molecules27051636 sha: d9c40aa419b80145a8a14e280a82968db47dca31 doc_id: 767421 cord_uid: pnqzvqk8 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 M(pro)-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). Natural products served for thousands of years as the major source for the basic needs of treatment and food in human life [1, 2] . Scientists investigated plants [3, 4] and microbes [5, 6] to deeply understand the causatives of these powers. Various metabolites have been isolated and exhibited promising activities as alkaloids [7] , flavonoids [8, 9] , isochromenes [10] , α-pyrones [11] , diterpenes [12] , sesquiterpenes [13] [14] [15] , steroids [16] , and saponins [17, 18] . 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 genus Artemisia was and still is an interesting target for the scientist to study. The isolation and identification of sesquiterpene lactones have been recorded severally in a lot of Artemisia species such as A. tschernieviana and A. sublessingiana [20] . Furthermore, the presence of flavonoids was also reported in A. albida [21] and A. santolinifolia [22] . The isolation of promising secondary metabolites was reported before as epi-ashantin from A. sieversiana [23] , anhydroaustricin from A. albida [24] , matricarin from A. austriaca [25] , cirsineol and cubreva lactone in A. umbrosa [26] . Recently, the in silico anti-SARS-CoV-2 papain-like protease activities of two derivatives of 2-Phenoxychromone of Artemisia spp. has been reported [27] . Additionally, chemical and biological properties of essential oils of diverse Artemisia species were discussed in various records as A. kasakorum [28] A. lercheana, A. sieversiana [29] , A. species [30] , A. umbrosa [31] , five Artemisia species [32] , A. gurganica [33] , A. proceriformis [34] , A. terrae-albae [35] , A. keiskeana [36] , A. littoricola, A. mandshurica [37] , and A. commutata [38] . 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). 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 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 C18H16O8 was decided by the HR-ESI-MS experiment, (+ve mode), that displayed a pseudo molecular ion peak [ Compound 1 was isolated as a yellow powder with a m.p. more than 350 °C. The 1 H NMR spectrum of 1 showed four singlet signals of aromatic protons at H 7.11 (s, H-3), H 6.97 (s, H-8), H 6.58 (s, H-3′), and H 7.45 (s, H-6′). Furthermore, three signals of methoxy groups were detected at H 3.74 (s), H 3.82 (s), and H 3.94 (s) ppm. The characteristic chelated signal of the OH at C-5 appeared as a sharp singlet at H 13.06 (due to intramolecular hydrogen bonding) (see Table 1 ). The 13 С, HSQC, and HMBC spectra of 1 indicated the presence of 15 carbon atoms in addition to 3 methoxy groups. The chemical structure of 1 was further confirmed by the heteronuclear multiple bond correlation (HMBC) experiment ( Figure 2 ). Two HMBC correlations were detected between the proton at H 7.11 (H, s, H-3) and C-10 (C 105.3) and C-1′ (C 107.2). Moreover, the OH at C-5 was engaged in three HMBC correlations with C-5 (δC 152.6), C-6 (δC 131.7), and C-10 (C 105.3). The ring B was confirmed through the HMBC correlations between the proton at H 6.58 (H, s, H-3′) and C-1′ (C 107.2) as well as C-5′ (C 139.7) in addition to the correlations between the proton at H 7.45 (H, s, H-6′) C-2′ (C 152.4), and C-4′ (C 142.0). The HMBC spectrum indicated the attachment of the methoxy groups at H 3.74 (s), H 3.82 (s), and H 3.94 (s) ppm to the carbon atoms C-6 (C 60.0), C-7 (C 56.7), and C-4′ (C 56.5), respectively. Compound 1 was isolated as a yellow powder with a m.p. more than 350 • C. The 1 H NMR spectrum of 1 showed four singlet signals of aromatic protons at δH 7.11 (s, H-3), δH 6.97 (s, H-8), δH 6.58 (s, H-3 ), and δH 7.45 (s, H-6 ). Furthermore, three signals of methoxy groups were detected at δH 3.74 (s), δH 3.82 (s), and δH 3.94 (s) ppm. The characteristic chelated signal of the OH at C-5 appeared as a sharp singlet at δH 13.06 (due to intramolecular hydrogen bonding) (see Table 1 ). The 13 C, HSQC, and HMBC spectra of 1 indicated the presence of 15 carbon atoms in addition to 3 methoxy groups. The chemical structure of 1 was further confirmed by the heteronuclear multiple bond correlation (HMBC) experiment ( Figure 2 ). Two HMBC correlations were detected between the proton at δH 7.11 (H, s, H-3) and C-10 (δC 105.3) and C-1 (δC 107.2). Moreover, the OH at C-5 was engaged in three HMBC correlations with C-5 (δC 152.6), C-6 (δC 131.7), and C-10 (δC 105.3). The ring B was confirmed through the HMBC correlations between the proton at δH 6.58 (H, s, H-3 ) and C-1 (δC 107.2) as well as C-5 (δC 139.7) in addition to the correlations between the proton at δH 7.45 (H, s, H-6 ) C-2 (δC 152.4), and C-4 (δC 142.0). The HMBC spectrum indicated the attachment of the methoxy groups at δH 3.74 (s), δH 3.82 (s), and δH 3.94 (s) ppm to the carbon atoms C-6 (δC 60.0), C-7 (δC 56.7), and C-4 (δC 56.5), respectively. Further chromatographic techniques led to the isolation of 2 (200 mg), 3 (20 m 4 (50 mg). The spectral data of 2-4 were found to correspond to capillartemisin methyl-3-[S-hydroxyprenyl]-cumarate [49] , and β-sitosterol [50] , respectively ( Fig Figure 3 . Chemical structures of compounds 2-4. Further chromatographic techniques led to the isolation of 2 (200 mg), 3 (20 mg), and 4 (50 mg). The spectral data of 2-4 were found to correspond to capillartemisin A [48] , methyl-3-[S-hydroxyprenyl]-cumarate [49] , and β-sitosterol [50] , respectively ( Figure 3 ). 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 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). The results revealed a high degree of similarity between the Jusanin and X77 (Table 2 and Figure 5 ). Discovery studio software was utilized to explore the following molecular properties in Jusanin and the nine ligands; partition coefficient (ALog p) [52] , molecular weight (M. Wt) [53] , hydrogen bond (H-bond) donors (HBA) [54] , H-bond acceptors (HBD) [55] , number of rotatable bonds [56] , number of rings, and also aromatic rings [57] , as well as molecular fractional polar surface area (MFPSA) [58] . 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] . 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 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) . 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 ). 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 ). 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. 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. 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 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. 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. 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. 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 ). 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%. 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 ). 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. 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. 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. 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. 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 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. 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. 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. 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). 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). The DFT parameters of Jusanin were calculated using Discovery Studio software [78] (See Supplementary Materials). 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). ADMET descriptors of Jusanin were determined using Discovery Studio 4.0 [82, 83] (See Supplementary Materials). The toxicity parameters of Jusanin were calculated using Discovery Studio 4.0 [84] [85] [86] (See Supplementary Materials). 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). 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. Traditional ancient Egyptian medicine: A review The Chinese herbal formulae (Yitangkang) exerts an antidiabetic effect through the regulation of substance metabolism and energy metabolism in type 2 diabetic rats Black ginseng and its saponins: Preparation, phytochemistry and pharmacological effects Chemical constituents of stems and leaves of Tagetespatula L. and its fingerprint Comparative biological evaluation of four endophytic fungi isolated from nigella sativa seeds New α-Pyrone derivatives from the endophytic fungus Embellisia sp Two new antileishmanial diketopiperazine alkaloids from the endophytic fungus Trichosporum sp Biological evaluation and molecular docking study of metabolites from Salvadora Persica L. Growing in Egypt Chemical Constituents of Tagetes patula and Their Neuroprotecting Action Nigrosphaerin A a new isochromene derivative from the endophytic fungus Nigrospora sphaerica Antileukemic α-pyrone derivatives from the endophytic fungus Alternaria phragmospora In Vitro and In Silico Cytotoxic and Antibacterial Activities of a Diterpene from Cousinia alata Schrenk Comprehensive In Silico Screening of the Antiviral Potentialities of a New Humulene Glucoside from Asteriscus hierochunticus against SARS-CoV-2 Antileishmanial Derivatives of Humulene from Asteriscus hierochunticus with in silico Synthesis and Molecular Docking of Some Grossgemin Amino Derivatives as Tubulin Inhibitors Targeting Colchicine Binding Site Isolation, Crystal Structure, and In Silico Aromatase Inhibition Activity of Ergosta-5, 22-dien-3β-ol from the Fungus Gyromitra esculenta Induction of apoptosis in human cancer cells through extrinsic and intrinsic pathways by Balanites aegyptiaca furostanol saponins and saponin-coated silvernanoparticles New combination approaches to combat methicillin-resistant Staphylococcus aureus (MRSA) Constituents of Artemisia tschernieviana and their biological activity Sesquiterpene lactones and flavonoids from Artemisia albida Flavonoids from Artemisia santolinifolia 4-Epiashantin from Artemisia sieversiana Anhydroaustricin from Artemisia albida Constituents of Artemisia austriaca and their biological activity Cirsilineol and cubreuva lactone from Artemisia umbrosa and their biological activity Isolation and In Silico Anti-SARS-CoV-2 Papain-like Protease Potentialities of Two Rare 2-Phenoxychromone Derivatives from Artemisia spp Component composition and antimicrobial activity of essential oil from Artemisia kasakorum Component composition of essential oils of Artemisia lercheana and A. sieversiana of the flora of Kazakhstan. Antimicrobial activity of A. sieversiana essential oil Essential oil from Kazakhstan Artemisia species Composition and Bioactivity of Artemisia umbrosa Essential Oil Chemical diversity and biological activity of the volatiles of five Artemisia species from Far East Russia Constituent Composition and Biological Activity of Essential Oil from Artemisia gurganica Chemical composition and antimicrobial activity of essential oils from Acantholippia deserticola, Artemisia proceriformis, Achillea micrantha and Libanotis buchtormensis against phytopathogenic bacteria and fungi Constituent Composition and Biological Activity of Essential Oil from Artemisia terrae-albae Composition and Biological Activity of Essential Oil from Artemisia keiskeana Constituent Compositions of Essential Oils from Artemisia littoricola and A. mandshurica The essential oil composition of six Artemisia species from Mongolia WHO. WHO Coronavirus (COVID-19) Dashboard. Available online Computational methods in drug discovery Isolation and In Silico Anti-COVID-19 Main Protease (Mpro) Activities of Flavonoids and a Sesquiterpene Lactone from Artemisia sublessingiana Molecular Docking and Dynamics Simulation Revealed the Potential Inhibitory Activity of ACEIs Against SARS-CoV-2 Targeting the hACE2 Receptor In Silico Studies of Some Isoflavonoids as Potential Candidates against COVID-19 Targeting Human ACE2 (hACE2) and Viral Main Protease (Mpro) Comprehensive virtual screening of the antiviral potentialities of marine polycyclic guanidine alkaloids against SARS-CoV-2 (COVID-19) Silico Exploration of Potential Natural Inhibitors against SARS-Cov-2 nsp10 Discovery of potential flavonoid inhibitors against COVID-19 3CL proteinase based on virtual screening strategy Artemisia Commutata Besser//Plantarium. Plants and Lichens of Russia and Neighboring Countries: Open Online Galleries and Plant Identification Guide Capillartemisin A and B, two new choleretic principles from Artemisiae capillaris Herba Prenylated coumarates from Artemisia xanthochroa Biosynthesis of β-sitosterol and stigmasterol in Croton sublyratus proceeds via a mixed origin of isoprene units Rabiu, I. Improved deep learning based method for molecular similarity searching using stack of deep belief networks An evaluation of in-silico methods for predicting solute partition in multiphase complex fluids-A case study of octanol/water partition coefficient An adverse outcome pathway for sensitization of the respiratory tract by low-molecular-weight chemicals: Building evidence to support the utility of in vitro and in silico methods in a regulatory context Effect of hydrogen bond donors and acceptors on CO 2 absorption by deep eutectic solvents In silico studies of diarylpyridine derivatives as novel HIV-1 NNRTIs using docking-based 3D-QSAR, molecular dynamics, and pharmacophore modeling approaches In silico strategies for modeling RNA aptamers and predicting binding sites of their molecular targets Ligand-Based Approach for In-silico Drug Designing Development of an in silico prediction model for chemical-induced urinary tract toxicity by using naïve Bayes classifier Predicting differential ion mobility behaviour in silico using machine learning Molecular hybridization conceded exceptionally potent quinolinyl-oxadiazole hybrids through phenyl linked thiosemicarbazide antileishmanial scaffolds: In silico validation and SAR studies In silico design of novel benzohydroxamate-based compounds as inhibitors of histone deacetylase 6 based on 3D-QSAR, molecular docking, and molecular dynamics simulations Structure based pharmacophore modeling, virtual screening, molecular docking and ADMET approaches for identification of natural anti-cancer agents targeting Embracing the diversity of halogen bonding motifs in fragment-based drug discoveryconstruction of a diversity-optimized halogen-enriched fragment library Calculating the aromaticity of heterocycles Lipophilicity of organic compounds calculated using structural similarity and molecular physicochemical descriptors Reoptimization of MDL keys for use in drug discovery Frontier Orbitals and Organic Chemical Reactions Ab initio HF, DFT and experimental (FT-IR) investigation of vibrational spectroscopy of PN, N-dimethylaminobenzylidenemalononitrile (DBM) On contribution of known atomic partial charges of protein backbone in electrostatic potential density maps Náray-szabó, G. A combined electronegativity equalization and electrostatic potential fit method for the determination of atomic point charges Molecular polarization maps as a tool for studies of intermolecular interactions and chemical reactivity Vibrational and electronic investigations, thermodynamic parameters, HOMO and LUMO analysis on Lornoxicam by density functional theory Synthesis, spectroscopic characterization, molecular docking, and ADMET studies of mannopyranoside esters as antimicrobial agents ADMET modeling approaches in drug discovery Prediction of drug absorption using multivariate statistics Molecular dynamics simulations Expression, Purification, and Comparative Inhibition of Helicobacter pylori Urease by Regio-Selectively Alkylated Benzimidazole 2-Thione Derivatives Discovery of new anticancer thiourea-azetidine hybrids: Design, synthesis, in vitro antiproliferative, SAR, in silico molecular docking against VEGFR-2, ADMET, toxicity, and DFT studies Anticancer activity, spectroscopic and molecular docking of some new synthesized sugar hydrazones, Arylidene and α-Aminophosphonate derivatives Discovery of new quinoxaline-2 (1H)-one-based anticancer agents targeting VEGFR-2 as inhibitors: Design, synthesis, and anti-proliferative evaluation Design, molecular docking, in vitro, and in vivo studies of new quinazolin-4 (3H)-ones as VEGFR-2 inhibitors with potential activity against hepatocellular carcinoma New quinoxaline-2 (1 H)-ones as potential VEGFR-2 inhibitors: Design, synthesis, molecular docking, ADMET profile and anti-proliferative evaluations Discovery of new quinazolin-4 (3H)-ones as VEGFR-2 inhibitors: Design, synthesis, and antiproliferative evaluation Design, synthesis, and anti-proliferative evaluation of new quinazolin-4 (3H)-ones as potential VEGFR-2 inhibitors Design, synthesis, molecular modeling, in vivo studies and anticancer activity evaluation of new phthalazine derivatives as potential DNA intercalators and topoisomerase II inhibitors Discovery and antiproliferative evaluation of new quinoxalines as potential DNA intercalators and topoisomerase II inhibitors A web-based graphical user interface for CHARMM The biomolecular simulation program CHARMM-GUI Input Generator for NAMD, GROMACS, AMBER, OpenMM, and CHARMM/OpenMM Simulations Using the CHARMM36 Additive Force Field Optimization of the additive CHARMM all-atom protein force field targeting improved sampling of the backbone phi, psi and side-chain chi(1) and chi(2) dihedral angles Scalable molecular dynamics with NAMD Comparison of simple potential functions for simulating liquid water Extension of the CHARMM General Force Field to sulfonyl-containing compounds and its utility in biomolecular simulations The authors are grateful to W. De Borggraeve and B. Demarsin (KULeuven, Belgium) for help in recording and interpreting mass-and NMR spectra. The authors declare no conflict of interest.Sample Availability: Samples of the compounds are available from the authors.