key: cord-0902889-xn5u5j82 authors: Suleimen, Yerlan M.; Jose, Rani A.; Suleimen, Raigul N.; Arenz, Christoph; Ishmuratova, Margarita; Toppet, Suzanne; Dehaen, Wim; Alsfouk, Aisha A.; Elkaeed, Eslam B.; Eissa, Ibrahim H.; Metwaly, Ahmed M. title: Isolation and In Silico Anti-SARS-CoV-2 Papain-Like Protease Potentialities of Two Rare 2-Phenoxychromone Derivatives from Artemisia spp. date: 2022-02-11 journal: Molecules DOI: 10.3390/molecules27041216 sha: 10b99c306b04265fd609b0c13b4922bc9f97a82a doc_id: 902889 cord_uid: xn5u5j82 Two rare 2-phenoxychromone derivatives, 6-demethoxy-4`-O-capillarsine (1) and tenuflorin C (2), were isolated from the areal parts of Artemisia commutata and A. glauca, respectively, for the first time. Being rare in nature, the inhibition potentialities of 1 and 2 against SARS-CoV-2 was investigated using multistage in silico techniques. At first, molecular similarity and fingerprint studies were conducted for 1 and 2 against co-crystallized ligands of eight different COVID-19 enzymes. The carried-out studies indicated the similarity of 1 and 2 with TTT, the co-crystallized ligand of COVID-19 Papain-Like Protease (PLP), (PDB ID: 3E9S). Therefore, molecular docking studies of 1 and 2 against the PLP were carried out and revealed correct binding inside the active site exhibiting binding energies of −18.86 and −18.37 Kcal/mol, respectively. Further, in silico ADMET in addition to toxicity evaluation of 1 and 2 against seven models indicated the general safety and the likeness of 1 and 2 to be drugs. Lastly, to authenticate the binding and to investigate the thermodynamic characters, molecular dynamics (MD) simulation studies were conducted on 1 and PLP. Natural products have been utilized in different folklores to treat many illnesses since ancient times [1] [2] [3] . The research in the field of natural product chemistry identified various bioactive metabolites that are engaged with the bioactivities such as diterpenes [4] , phenolics [5] , sesquiterpenes [6, 7] , steroids [8] , flavonoids [9, 10] , alkaloids [11] , saponins [12] [13] [14] , isochromenes [15] , benzofuran glucosides [16] , and α-pyrones [17, 18] . The chemical composition as well as the biological activities of plants that belong to the genus Artemisia L. (family Asteraceae) drew the attention of scientists to explore and investigate. For example, anhydroaustricin, epiashantin, and matricarin are biologically active compounds that were isolated from A. albida [19] , A. sieversiana [20] , and A. austriaca [21] , respectively. Besides, cirsineol is an active lactone that was reported in A. umbrosa [22] . Also, the presence of sesquiterpene lactones and flavonoids were reported in A. tschernieviana [23] , A. albida [24] , and A. santolinifolia [25] , respectively. On the other hand, the bioactivities and phyto analysis of essential oils were studied for A. kasakorum [26] , A. lercheana and A. sieversiana [27] , A. umbrosa [28] , A. gurganica [29] , A. proceriformis [30] , A. terrae-albae [31] , A. keiskeana [32] , A. littoricola, A. mandshurica [33] , A. commutata [34] , in addition to five plants belonging to the Artemisia species [35] . A. commutata [36] and A. glauca [37] are perennial herbs that are cultivated in northern and eastern Kazakhstan [38] . The lack of scientific records about the phytochemistry and the biological effects of these two species encouraged us to examine their phytochemical properties aiming at the isolation of promising secondary metabolites. The WHO declared on December 28, 2021 , that the confirmed cases of COVID-19 infections on a global base are 280,119,931. Sadly, 5,403,662 of them are dead [39] . In regard to these bad numbers, massive work is demanded from scientists all over the world to find a cure. The field of computer-aided/based (in silico or computational) drug design and discovery is a fast-growing approach that has shown several successes and advancements during the last decade [40] . The structural informatics explosion allowed the identification of the exact chemical and physical properties of molecules. Consequently, the ability to compare different structures and predict a specific biological activity became much easier [41] . Also, the huge advancement in the field of proteomic enabled scientists to expect the biding and hence the activity of a molecule and target enzyme [42] . Our teamwork reported the utilization of the in silico methods to suggest a treatment that could be effective against COVID-19 in many reports [43] [44] [45] [46] [47] . We here in this manuscript, report the isolation and structure elucidation of two rare 2-phenoxychromone derivatives, 6-demethoxy-4'-O-capillarsine (1) and tenuflorin C (2) from the areal parts of A. commutata and A. glauca, respectively, for the first time. Because of the lack of research on the isolated compounds and the severe need to find a cure against COVID-19, the potential inhibitory effects of 1 and 2 against COVID-19 were investigated depending on a multiphase in silico screening method. The potentialities of 1 and 2 against COVID-19 Papain-Like Protease were confirmed according to molecular similarity, structure fingerprint, molecular docking, and molecular dynamics simulations studies. Additionally, ADMET and toxicity properties of 1 and 2 were studied to predict their likeness to be used as drugs. The areal parts of A. commutata, (940 g) and A. glauca (1.02 kg) were collected from the western part of the Altai Mountains. The dried plants were extracted with chloroform. The chloroform extract of A. commutata was subjected to chromatographic purifications steps to afford the amorphous powder of 1. Compound 1 was identified to be 6-demethoxy-4'-O-capillarsine ( Figure 1 ). Identification of 1 was done depending on 1 H and 13 C NMR (Table S1, in the Supplementary Materials), comparing the published data. Compound 1 was previously isolated from the leaves of Mimosa tenuiflora [48] , A. capillaris [49] , and A. rupestris [50] . The chloroform extract of A. glauca was subjected to chromatographic purification steps to give the amorphous powder of 2. Compound 2 was identified to be 5,7,4 -trihydroxy, 3 -methoxy-2-phenoxychromone ( Figure 1 ). Compound 2 was identified depending on 1 H and 13 C NMR (Table S2 , in the Supplementary Materials) compared to the published data. Compound 2 was previously isolated from Mimosa tenuiflora [48] . previously isolated from the leaves of Mimosa tenuiflora [48] , A. capillaris [49] , and A. rupestris [50] . The chloroform extract of A. glauca was subjected to chromatographic purification steps to give the amorphous powder of 2. Compound 2 was identified to be 5,7,4′trihydroxy, 3′-methoxy-2-phenoxychromone ( Figure 1 ). Compound 2 was identified depending on 1 H and 13 C NMR (Table S2 , in the supporting data) compared to the published data. Compound 2 was previously isolated from Mimosa tenuiflora [48] . . To understand the principle of molecular similarity, we have to understand that the bioactivity of any compound is a result of some well-known interactions with a specific protein target. These protein-ligand interactions depend on specific physical and chemical interactions such as hydrogen bonding and hydrophobic interactions. Consequently, the similarity in the chemical structures of two compounds will give a similarity in the number and position of hydrogen bond donors, acceptors, and hydrophobic centers in addition to the steric configuration. This similarity could cause a noticeable similarity in the bioactivity too [51] . To evaluate the expected anti-COVID-19 potentialities of compounds 1 and 2, the chemical structures of them were compared with the chemical structures of eight different co-crystallized ligands of eight crucial COVID-19 proteins ( Figure 2 ) in a structural similarity experiment. Utilizing Discovery Studio software, the following molecular properties were measured to 6-demethoxy-4`-O-capillarsine (1) and tenuflorin C (2) as well as the eight ligands; partition coefficient (ALog p) [52] , molecular weight (M. W) [53] , Hbond donors (HBD) [54] and acceptors (HBA) [55] , rotatable bonds [56] , aromatic rings [57] , heterocyclic rings [58] , along with the molecular fractional polar surface area (MFPSA) [59] . The results revealed the great similarity between compounds 1 and 2 and the co-crystallized ligands (TTT) of SARS-CoV-PLP (PDB ID: 3E9S) (Table 1, Figure 3A ,B). To understand the principle of molecular similarity, we have to understand that the bioactivity of any compound is a result of some well-known interactions with a specific protein target. These protein-ligand interactions depend on specific physical and chemical interactions such as hydrogen bonding and hydrophobic interactions. Consequently, the similarity in the chemical structures of two compounds will give a similarity in the number and position of hydrogen bond donors, acceptors, and hydrophobic centers in addition to the steric configuration. This similarity could cause a noticeable similarity in the bioactivity too [51] . To evaluate the expected anti-COVID-19 potentialities of compounds 1 and 2, the chemical structures of them were compared with the chemical structures of eight different co-crystallized ligands of eight crucial COVID-19 proteins ( Figure 2 ) in a structural similarity experiment. Utilizing Discovery Studio software, the following molecular properties were measured to 6-demethoxy-4'-O-capillarsine (1) and tenuflorin C (2) as well as the eight ligands; partition coefficient (ALog p) [52] , molecular weight (M. W) [53] , H-bond donors (HBD) [54] and acceptors (HBA) [55] , rotatable bonds [56] , aromatic rings [57] , heterocyclic rings [58] , along with the molecular fractional polar surface area (MFPSA) [59] . The results revealed the great similarity between compounds 1 and 2 and the co-crystallized ligands (TTT) of SARS-CoV-PLP (PDB ID: 3E9S) (Table 1, Figure 3A ,B). The molecular fingerprint is an in silico method that analyzes the similarities of two molecules or more. Molecular fingerprint explores the property profiles of a certain compound, the technique computes these properties in the forms of bits vectors examining the existence, absence, and frequencies it in the reference and target compounds. The molecular fingerprint is an in silico method that analyzes the similarities of two molecules or more. Molecular fingerprint explores the property profiles of a certain compound, the technique computes these properties in the forms of bits vectors examining the existence, absence, and frequencies it in the reference and target compounds. The fingerprint approach examined the existence or absence of the these parameters: charge [60] , hybridization [61] , H-bond acceptor and donor [62] , negative and positive ionizables [63] , halogen [64] , aromatic [65] , and ALogP [66] . The experiment was achieved depending on Discovery Studio. The results verified the similarity of 6-demethoxy-4'-Ocapillarsine (1), tenuflorin C (2), and TTT (Table 2 ). The reported main pharmacophoric features of SARS-CoV-2 Papain-Like protease inhibitors (PLPIs) are two hydrophobic systems separated be a linker [67] . As shown in Figure 4 , the co-crystallized ligand (TTT) has naphthyl and phenyl groups as hydrophobic centers. The two hydrophobic groups are separated by hydrophilic linker (amide group). In addition, the terminal phenyl group is substituted by amino group as a hydrogen bonding group. Similarly, compounds 1 and 2 have coumarin moiety as the first hydrophobic center in addition to phenyl group as the second hydrophobic moiety. In each compound, the two hydrophobic centers were separated by oxygen atom as a hydrophilic linker. This linker facilitates these compounds to have the same configuration of TTT. Additionally, the two phenyl rings were substituted by hydrophilic groups as which can form hydrogen bonds with the target receptors as TTT. From medicinal chemistry point of view, this investigation revealed that there is a great similarity between the two tested molecules and the co-crystallized ligand (TTT). centers. The two hydrophobic groups are separated by hydrophilic linker (amide group). In addition, the terminal phenyl group is substituted by amino group as a hydrogen bonding group. Similarly, compounds 1 and 2 have coumarin moiety as the first hydrophobic center in addition to phenyl group as the second hydrophobic moiety. In each compound, the two hydrophobic centers were separated by oxygen atom as a hydrophilic linker. This linker facilitates these compounds to have the same configuration of TTT. Additionally, the two phenyl rings were substituted by hydrophilic groups as which can form hydrogen bonds with the target receptors as TTT. From medicinal chemistry point of view, this investigation revealed that there is a great similarity between the two tested molecules and the co-crystallized ligand (TTT). 3D-Flexible alignment of the compounds 1 and 2 with TTT was presented in Figure 5 . From the figure, it is possible to observe that, in general, the structure of compounds 1 and 2 have a good overlap with the reference molecule (TTT). In addition, the two tested compounds showed the same spatial orientation of TTT due to the presence of the flexible linker at the center of these compounds. 3D-Flexible alignment of the compounds 1 and 2 with TTT was presented in Figure 5 . From the figure, it is possible to observe that, in general, the structure of compounds 1 and 2 have a good overlap with the reference molecule (TTT). In addition, the two tested compounds showed the same spatial orientation of TTT due to the presence of the flexible linker at the center of these compounds. To examine the obtained results from structural similarity and fingerprint studies, molecular docking experiments were achieved for 6-demethoxy-4`-O-capillarsine (1) and tenuflorin C (2) against PLP (PDB ID: 3E9S) using TTT, the co-crystallized ligand, as a reference. We considered the binding free energies (∆G) besides the binding modes as the bases of evaluation. Firstly, a validation process has proceeded via the re-docking of TTT against PLP. To examine the obtained results from structural similarity and fingerprint studies, molecular docking experiments were achieved for 6-demethoxy-4'-O-capillarsine (1) and tenuflorin C (2) against PLP (PDB ID: 3E9S) using TTT, the co-crystallized ligand, as a reference. We considered the binding free energies (∆G) besides the binding modes as the bases of evaluation. Firstly, a validation process has proceeded via the re-docking of TTT against PLP. The obtained RMSD value between the two poses was 2.1 • A confirming the validity of the docking method ( Figure 6 ). To examine the obtained results from structural similarity and fingerprint studies, molecular docking experiments were achieved for 6-demethoxy-4`-O-capillarsine (1) and tenuflorin C (2) against PLP (PDB ID: 3E9S) using TTT, the co-crystallized ligand, as a reference. We considered the binding free energies (∆G) besides the binding modes as the bases of evaluation. Firstly, a validation process has proceeded via the re-docking of TTT against PLP. The obtained RMSD value between the two poses was 2.1 ᵒ A confirming the validity of the docking method ( Figure 6 ). The binding mode of TTT showed a free energy of −20.32 kcal/mol. The p-toluidine part occupied the first pocket of PLP forming one H-bond with the amino acid, Leu163. The same moiety was engaged hydrophobically in two interactions with the amino acids Asp165 and Tyr269. In addition, the amide linker formed two H-bonds with the amino acids Asp165 and Gln270. The naphthalene moiety was buried in the second pocket making five hydrophobic interactions with Pro249, Pro248, and Tyr269 ( Figure 7A The binding mode of 6-demethoxy-4'-O-capillarsine (1) exhibited a binding free energy of -18.86 kcal/mol. The anisole moiety was directed into the first pocket of PLP forming one H-bond with the amino acid, Gln270, in addition to two hydrophobic interactions with the amino acids, Tyr269 and Asp165. The 5,7-dihydroxy-4H-chromen-4-one part occupied the second pocket of PLP forming a H-bond with the amino acid Pro249 besides two hydrophobic bonds with the amino acids Pro248 and Pro249 ( Figure 8A-C) . Tenuflorin C (2) bonded to PLP exhibiting a binding energy of −18.37 kcal/mol. The 2-methoxyphenol part was directed into the first pocket of PLP with two H-bonds with the amino acids Gln270 and Leu163. Additionally, the same moiety was engaged hydrophobically in two interactions with the amino acids Tyr269 and Asp165. Further, the 5,7-dihydroxy-4H-chromen-4-one part occupied the second pocket of PLP forming two H- The binding mode of 6-demethoxy-4`-O-capillarsine (1) exhibited a binding free energy of -18.86 kcal/mol. The anisole moiety was directed into the first pocket of PLP forming one H-bond with the amino acid, Gln270, in addition to two hydrophobic interactions with the amino acids, Tyr269 and Asp165. The 5,7-dihydroxy-4H-chromen-4-one part occupied the second pocket of PLP forming a H-bond with the amino acid Pro249 besides two hydrophobic bonds with the amino acids Pro248 and Pro249 ( Figure 8A Tenuflorin C (2) bonded to PLP exhibiting a binding energy of −18.37 kcal/mol. The 2-methoxyphenol part was directed into the first pocket of PLP with two H-bonds with the amino acids Gln270 and Leu163. Additionally, the same moiety was engaged hydrophobically in two interactions with the amino acids Tyr269 and Asp165. Further, the 5,7dihydroxy-4H-chromen-4-one part occupied the second pocket of PLP forming two Hbonds with Tyr274 and Ala247. Finally, it was incorporated hydrophobically in three bonds with the amino acids Pro248, Asp165 and Met209 ( Figure 9A-C) . There is always a need to examine the ADMET properties of new molecules especially in the stage of designing to decrease the possibility of late-stage attrition. Also, the identification of ADMET properties provides good information regarding the amount and frequency of dosing as well as the toxicity [68] . The in silico ADMET parameters were determined for 6-demethoxy-4`-O-capillarsine (1) and tenuflorin C (2) using Discovery Studio software and utilizing Remdesivir as a reference. The results were illustrated in Figure 10 and Table 3 . The results revealed that 6-demethoxy-4`-O-capillarsine (1) and tenuflorin C (2) have low and very low chance to penetrate BBB, respectively. Further, 6-demethoxy-4`-O-capillarsine (1) and tenuflorin C (2) exhibited good aqueous solubility and intestinal absorption. Finally, 6-demethoxy-4`-O-capillarsine (1) and tenuflorin C (2) were predicted to be a CYP2D6 non-inhibitor and to be able to bind plasma protein by a ratio of more than 90%. There is always a need to examine the ADMET properties of new molecules especially in the stage of designing to decrease the possibility of late-stage attrition. Also, the identification of ADMET properties provides good information regarding the amount and frequency of dosing as well as the toxicity [68] . The in silico ADMET parameters were determined for 6-demethoxy-4'-O-capillarsine (1) and tenuflorin C (2) using Discovery Studio software and utilizing Remdesivir as a reference. The results were illustrated in Figure 10 and Table 3 . The results revealed that 6-demethoxy-4'-O-capillarsine (1) and tenuflorin C (2) have low and very low chance to penetrate BBB, respectively. Further, 6-demethoxy-4'-Ocapillarsine (1) and tenuflorin C (2) exhibited good aqueous solubility and intestinal absorption. Finally, 6-demethoxy-4'-O-capillarsine (1) and tenuflorin C (2) were predicted to be a CYP2D6 non-inhibitor and to be able to bind plasma protein by a ratio of more than 90%. The predicted toxicity levels of 6-demethoxy-4'-O-capillarsine (1) and tenuflorin C (2) were investigated in silico using Discovery Studio software against seven toxicity models. Remdesivir was used as a reference. As shown in Table 4 , 6-demethoxy-4'-O-capillarsine (1) and tenuflorin C (2) were predicted to be non-carcinogenic and have high TD 50 values of 144.939 and 113.277 mg/kg/day, respectively. The values were more than 12-fold that of remdesivir (9.2458 mg/kg/day). Further, 6-demethoxy-4'-O-capillarsine (1) and tenuflorin C (2) showed high rat maximum tolerated dose (MTD) values of 0.289735 and 0.44318 g/kg, respectively, compared to remdesivir (0.234965 g/kg). The rat oral LD 50 values for 1 and 2 were 0.363122 and 0.549081 g/kg, respectively, and regarding the rat chronic LOAEL model, their values were 0.0221904 and 0.0362493, respectively. Finally, 6-demethoxy-4'-O-capillarsine (1) and tenuflorin C (2) were predicted to have mild ocular irritancy and to have no irritancy against the skin model. The main advantage of molecular dynamics (MD) simulation is its ability to recognize the flexibility of the protein-ligand complex. This advantage enables the MD studies to estimate accurately the thermodynamics and kinetics that take place during the drugenzyme binding [69] . MD simulations can explore the dynamic structural information of the protein-ligand complex and provide plentiful information on the energetic changes that resulted from the interactions between protein and ligand. Such information could be very useful to understand the structure-function relationship of the protein-ligand complex interactions [70] . In order to authenticate the binding and to investigate the thermodynamic characters of the isolated compounds, molecular dynamics (MD) simulation studies were conducted to 6-demethoxy-4'-O-capillarsine, 1, as it showed a bitter binding energy, against PLP. The atomical dynamic movements and conformational variations of backbone atoms of the 6-demethoxy-4'-O-capillarsine -PLP complex were calculated by RMSD to detect their stability upon ligand bonded and apo states. It was observed that the protein and ligand exhibit lower RMSD with no major fluctuations indicating their greater stability. The 6-demethoxy-4'-O-capillarsine -PLP complex was slightly fluctuating till 70 ns~and stabled later ( Figure 11A ). The flexibility of PLP was evaluated in terms of RMSF ( Figure 11B) to get more information about the protein regions that had been fluctuating during the simulation. It can be understood that the binding of 6-demethoxy-4'-O-capillarsine does not make PLP very flexible. The compactness of the 6-demethoxy-4'-O-capillarsine -PLP complex was evaluated through the examination of the radius of gyration (Rg) ( Figure 11C ). Decreasing fluctuation during the simulation study is an indication of the higher compactness of the system. The Rg of the 6-demethoxy-4'-O-capillarsine -PLP complex was found to be similar throughout the simulation period. Interaction between 6-demethoxy-4'-O-capillarsine -PLP complexes and solvents was evaluated through the calculation of solvent accessible surface area (SASA) over the simulation period. So, SASA of the 6-demethoxy-4'-O-capillarsine -PLP complex was calculated to analyze the extent of the conformational changes that occurred during the interaction. Interestingly, PLP featured a reduction of the surface area showing a relatively lower SASA value than the starting period ( Figure 11D ). Hydrogen bonding in the 6-demethoxy-4'-O-capillarsine -PLP complex is essential to stabilize the structure. It was noticed that the highest number of conformations of the PLP formed up to three H-bonds with the 6-demethoxy-4'-O-capillarsine ( Figure 11E ). A total of 0.94 kg of Artemisia commutata and 1.04 kg of Artemisia glauca areal parts were extracted with chloroform. Successive chromatographic methods led to the isolation of compounds (1 and 2) . Detailed description in the method part in supporting data. Molecular Similarity of 6-demethoxy-4`-O-capillarsine (1), tenuflorin C (2) against the co-crystallized ligands of SARS-CoV-2 was carried out calculated using Discovery Studio 4.0 (See the method part in supporting data). A total of 0.94 kg of Artemisia commutata and 1.04 kg of Artemisia glauca areal parts were extracted with chloroform. Successive chromatographic methods led to the isolation of compounds (1 and 2) . Detailed description in the method part in Supplementary Materials. Molecular Similarity of 6-demethoxy-4'-O-capillarsine (1), tenuflorin C (2) against the co-crystallized ligands of SARS-CoV-2 was carried out calculated using Discovery Studio 4.0 (See the method part in Supplementary Materials). A fingerprint study of 6-demethoxy-4'-O-capillarsine (1), tenuflorin C (2) against the eight co-crystallized ligands of SARS-CoV-2 was carried out calculated using Discovery The DFT parameters were calculated for 6-demethoxy-4'-O-capillarsine (1), tenuflorin C (2) using Discovery Studio software [71] . (See the method part in Supplementary Materials). The docking investigation was carried out for 6-demethoxy-4'-O-capillarsine (1), tenuflorin C (2) using MOE2014 software and Discovery Studio 4.0 [72] [73] [74] (See the method part in Supplementary Materials). ADMET descriptors of 6-demethoxy-4'-O-capillarsine (1), tenuflorin C (2) were determined using Discovery Studio 4.0. [75, 76] (See the method part in Supplementary Materials). Seven toxicity parameters of 6-demethoxy-4'-O-capillarsine (1), tenuflorin C (2) were calculated using Discovery Studio 4.0 [77] [78] [79] (See the method part in Supplementary Materials). The system was prepared using the web-based CHARMM-GUI [80] [81] [82] interface utilizing CHARMM36 force field [83] and NAMD 2.13 [84] package. The TIP3P explicit solvation model was used (See the method part in Supplementary Materials). 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 Quantitative knowledge presentation models of traditional Chinese medicine (TCM): A review In Vitro and In Silico Cytotoxic and Antibacterial Activities of a Diterpene from Cousinia alata Schrenk Biological and chemical evaluation of some African plants belonging to Kalanchoe species: Antitrypanosomal, cytotoxic, antitopoisomerase I activities and chemical profiling using ultra-performance liquid chromatography/quadrupoletime-of-flight mass spectrometer Antileishmanial Derivatives of Humulene from Asteriscus hierochunticus with in silico Tubulin Inhibition Potential 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 Biological evaluation and molecular docking study of metabolites from Salvadora Persica L. Growing in Egypt Chemical Constituents of Tagetes patula and Their Neuroprotecting Action Two new antileishmanial diketopiperazine alkaloids from the endophytic fungus Trichosporum sp 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) Black ginseng and its saponins: Preparation, phytochemistry and pharmacological effects Nigrosphaerin A a new isochromene derivative from the endophytic fungus Nigrospora sphaerica Chemical constituents of stems and leaves of Tagetespatula L. and its fingerprint Antileukemic α-pyrone derivatives from the endophytic fungus Alternaria phragmospora New α-Pyrone derivatives from the endophytic fungus Embellisia sp Anhydroaustricin from Artemisia albida 4-Epiashantin from Artemisia sieversiana Constituents of Artemisia austriaca and their biological activity Cirsilineol and cubreuva lactone from Artemisia umbrosa and their biological activity Constituents of Artemisia tschernieviana and their biological activity Sesquiterpene lactones and flavonoids from Artemisia albida Flavonoids from Artemisia santolinifolia 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 Composition and Bioactivity of Artemisia umbrosa Essential Oil 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 Chemical Diversity and Biological Activity of the Volatiles of Five Artemisia Species from Far East Russia WHO. WHO Coronavirus (COVID-19) Dashboard. Available online Computer-Aided Drug Design Current trends in computer aided drug design and a highlight of drugs discovered via computational techniques: A review 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 New 2-Phenoxychromones from the Leaves of Mimosa tenuiflora 2-Phenoxychromones and a structurally related flavone from leaves ofRosa rugosa A 2-phenoxychromone from Artemisia rupestris 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 Feature trees: A new molecular similarity measure based on tree matching 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 Silico Screening of Semi-Synthesized Compounds as Potential Inhibitors for SARS-CoV-2 Papain-Like Protease: Pharmacophoric Features, Molecular Docking, ADMET, Toxicity and DFT Studies Prediction of ADMET properties Molecular dynamics simulations Molecular dynamics simulations and novel drug discovery 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 Data Availability Statement: Data is contained within the article. The authors declare no conflict of interest.