Plants produce a variety of structurally diverse and often complex secondary metabolites (SM) with a range of biological functions [1]. Plant SM can be classified according to several criteria, such as chemical structure (presence of rings or sugars), composition (containing nitrogen or not), solubility in organic solvents or water, and the biosynthetic pathway [2]. They include alkaloids, carotenoids, organosulfur compounds, phytosterols, nitrogen compounds, and phenolics (Figure 1).

Each class is then divided into further classes. In particular, phenolics comprise a large group of different compounds, with the phenolic hydroxyl groups being the common structural feature. These are usually found conjugated with sugars and organic acids and can be divided into seven main classes such as hydrolyzable tannins, coumarins, lignans, phenolic acids, stilbenes, curcuminoids, and flavonoids [3] (Figure 1). Flavonoids exist broadly in nature, and to date, more than 9000 flavonoids have been reported. They can be divided into the following subclasses: flavonols, flavan-3-ols, isoflavones, anthocyanins, flavanones, and flavones. Flavonoids share a common structure of two benzene rings connected by three carbon atoms forming an oxygenated heterocycle and are responsible for the red, blue, and yellow coloration of plants and are found in foods and beverages of plant origin, such as fruits, vegetables, tea, cocoa, and wine ^[4][5][4,5]. Recent studies have reported that some of these natural compounds can reduce the risk of many chronic diseases and can significantly modulate and diversify the composition of the human gut microbiome [6]. Consumption of phenolic compounds can help to prevent metabolic, cardiovascular, respiratory, neurological, and cancerous diseases ^[7][8][7,8] due to their high antimicrobial, antioxidant, and anti-inflammatory immunomodulatory activities ^[9][10][11][9,10,11]. Eating plant-based foods is part of many diverse dietary patterns, including the well-studied Mediterranean diet [12], vegan, and vegetarian approaches. Despite the widely known health benefits of consuming fruits and vegetables, the intake is often inadequate and a large number of adults worldwide do not consume the daily-recommended servings. So, in the attempt to improve health, dietary guidelines and health promotion campaigns have advocated for individuals to “eat a rainbow” of fruit and vegetables, i.e., to take a qualitative color rather than a quantitative servings approach, based on the association of each color with a health benefit. For example, red foods include antioxidants that can contribute to decreased inflammation in the body; orange foods are abundant in carotenoids and have been linked to endocrine-regulating activities; yellow foods have been found to aid in digestion due to their fiber content and bioflavonoids that promote healthy gut bacteria; green foods, especially green leafy vegetables, contain an abundance of polyphenols that aid in reducing cardiovascular risk factors such as high blood pressure; blue or purple foods have been found to improve memory and mood because of flavonoids, flavonols, and phenolic acids that promote cognitive functions.

Recent technological advances in analytical methods such as metabolomics, metabolic engineering, and synthetic biology, as well as the ad hoc designed computational tools and databases, are providing powerful tools for drug discovery based on natural compounds [6]. It has been known that natural products and plant extracts exhibit potent anti-viral activities, often against multiple virus families, thus suggesting that they may be useful as broad-spectrum antiviral agents ^[13][14][13,14].

The COVID-19 pandemic caused by SARS-CoV-2 has prompted researchers to conduct many studies aimed at identifying useful compounds to counter the viral agent. As reported by the World Health Organization’s weekly epidemiological update dated 22 February 2023, “over 757 million confirmed cases and over 6.8 million deaths have been reported globally” [15], with nearly 5.3 million new cases and over 48,000 deaths in the 4 weeks preceding the report. For a state-of-the-art review of the pandemic two years after its appearance and the related lessons learned, an interesting article was published in September 2022 [16].

2. Antiviral Potential of Selected Natural Phytochemicals

Natural products serve as excellent sources for discovering antiviral agents due to their diversity and complexity and can offer remarkable efficacy and specificity to target viral infections ^[17][118]. The antiviral activity of a drug can affect viral entry, viral DNA and RNA synthesis, or viral reproduction, and is conditioned by viral structure and its replication cycle. For example, the presence or absence of a viral envelope, and, consequently, the different modes of entry of the virus into the host cell, play a significant role in the effectiveness of the virucidal activity, as the surface of the external membrane of the virus establishes the first contact with the drug. The antiviral activity of secondary metabolites can be evaluated by various biological assays to test the cytotoxicity, cytopathic effect inhibition, and the capability to block viral spread from cell to cell, thus limiting and/or fighting viral circulation ^[18][132]. Several natural products exhibit antiviral activities towards DNA and RNA viruses by acting on different cellular/viral targets and interfering with the infection and the viral replication cycle. Many of the secondary metabolites that are object of this review have shown antiviral activity against the Herpes simplex virus (HSV). HSV belongs to Herpesviridae, which is a broad family of enveloped-DNA viruses and is categorized into two types: HSV-1, commonly associated with orofacial ulcers, and HSV-2, which mainly causes genital ulcers ^[19][101]. For example, caffeic acid inhibits Herpes simplex type 1 (HSV-1) multiplication mainly before the completion of viral DNA replication but not thereafter ^[20][17]. Gallic acid, also belonging to the class of phenolic acids, has shown activity against HSV ^[21][29]. In the context of flavonoids, luteolin—belonging to the class of flavones—exhibited potent antiviral activity against wild-type and clinical isolates of HSV-2 ^[22][43], as well as kaempferol and its glycosides kaempferol-3-O-robinobioside and kaempferol-3-O-rutinoside, myricetin, quercetin, and its glycosides quercetin-3-O-rutinoside (rutin), and isorhamnetin, all belonging to the class of flavonols that have shown inhibitory activity against HSV-1 and/or HSV-2 ^{[21][23][24][25][26]}[29,38,39,60,71]. The flavanols epicatechin, (-)-epicatechin-3-O-gallate (ECG), (-)-epigallocatechin (EGC), and epigallocatechin-3-gallate (EGCG) and the flavanone naringenin possess antiherpetic activity ^[21][27][29,83]. In vitro studies have also been performed on molecules belonging to other classes of phenols such as curcumin ^[19][101], resveratrol ^[28][105], chebulagic acid ^[29][113], and other secondary metabolites such as allicin ^[30][126] and eugenol ^[31][128], demonstrating their antiviral activity against HSV-1 and HSV-2. Furthermore, studies have investigated the antiviral activity of some secondary metabolites against the Hepatitis C virus (HCV), a virus belonging to the Flaviviridae family that represents one of the major causes of chronic liver disease. It has been demonstrated that caffeic acid, as well as the flavone apigenin, could inhibit HCV replication ^[23][32][19,38], while gallic acid, quercetin-3-O-rutinoside (rutin), epigallocatechin-3-gallate (EGCG), and curcumin could inhibit HCV entry ^{[19][27][33][34]}[25,77,83,101]. In hepatitis virus infection, coumarin has been shown to target a wide range of proteins, like binding antigens present at the cell surface, proteins involved in viral replication, and interferon signaling pathways ^[35][44]. Furthermore, it has been demonstrated that sulforaphane, an isothiocyanate widely distributed in Brassicaceae plants, suppresses the replication of HCV by inducing the heme-oxygenase-1 (HO-1) expression, an enzyme that interferes with the replication of the virus through the activation of nuclear factor (erythroid-derived 2)-like (Nrf2) pathway, which regulates the expression of antioxidant proteins ^[36][124]. Adenoviruses (ADVs) are DNA viruses that typically cause mild infections involving the upper or lower respiratory tract, gastrointestinal tract, or conjunctiva. More than 50 serotypes have been identified showing different tissue tropisms that correlate with clinical manifestations of infection ^[37][133]. In a study conducted to examine the antiviral activity of aqueous extract and pure compounds of Plantago major L., a popular traditional Chinese medicine, it was found that chlorogenic acid was active against ADV-3, ADV-8, and ADV-11; instead, caffeic acid was active against ADV-3 ^[38][18]. Epigallocatechin-3-gallate (EGCG) has been shown to reduce the virus titers of adenovirus in two cell infection models, to inactivate purified adenovirions and to inhibit the attachment of adenovirus by interacting with virion surface proteins ^[27][83]. Caffeic acid also inhibits the multiplication of the Influenza A (IAV) virus in the early stage of viral multiplication ^[39][20]. Influenza virus (IV) is a single-stranded RNA virus, which belongs to the Orthomyxoviridae family, that causes acute respiratory infection. There are four types of influenza viruses: A, B, C, and D. Human influenza A and B viruses cause seasonal epidemics of disease; influenza C virus infections generally cause mild illness and are not thought to cause human epidemics; influenza D viruses primarily affect cattle and are not known to infect or cause human illness. Type A virus is divided into subtypes according to the antigenic properties of surface glycoproteins, hemagglutinin (HA), and neuraminidase (NA). In relation to genetic and antigenic variation, the current 18 subtypes of IAV HA (H1–H18) are divided into group 1 (H1, H2, H5, H6, H8, H9, H11, H12, H13, H16, H17, and H18) and group 2 (H3, H4, H7, H10, H14, and H15), while IAV of NA has 11 subtypes (N1–N11) ^[40][134]. Gallic acid has been shown to exhibit inhibitory effects against both influenza type A and B viruses, disrupting the viral particles ^[41][28], while hydroxytyrosol inactivates influenza A viruses, including H1N1, H3N2, H5N1, and H9N2 subtypes with an effect virucidal rather than antiviral ^[42][35]. Apigenin, luteolin, kaempferol, epicatechin, (-)-epicatechin-3-O-gallate (ECG), epigallocatechin-3-gallate (EGCG), cirsimaritin, and eugenol also show activity against IV inhibiting replication ^{[23][24][27][43][44]}[38,39,50,83,129]. An in vitro and in vivo study demonstrated that isorhamnetin exerts anti-influenza effects via direct HA and NA inhibition, direct or indirect inhibition of the expression of viral HA and NA genes as well as reduces virus-induced ROS (reactive oxygen species) generation ^[45][53]. Several in vitro studies demonstrated that curcumin inhibits the uptake, replication, and particle production of different IAVs ^[19][101]. In an in vitro and in silico study, Sang et al. ^[46][63] demonstrated that myricetin has not only a direct antiviral activity against IAV-inhibiting virus replication but also has immune modulatory effects. Wu et al. ^[47][73] demonstrated that quercetin inhibits a wide spectrum of influenza groups, including H1N1, H3N2, and H5N1, binding to influenza hemagglutinin protein and inhibiting viral-cell fusion. Instead, quercetin 3-rhamnoside (quercitrin) has a strong antiviral activity against the influenza A/WS/33 (H1N1) virus by inhibiting replication in the initial phase of the infection by indirect interaction with viral particles ^[48][78]. Sahoo et al. ^[49][87] used in silico docking approaches to explore the molecular interaction of 13 active compounds extracted from plants against NA protein of H1N1 and observed that theaflavin, found in green tea, could be a potential inhibitor of H1N1 NA proteins, as strongly suggested by lowest docking energy. Molecular docking studies have also shown the inhibitory role against IV infection of coumarin and its derivatives that can target various enzymes and pathways essential for viral entry, survival, and infection ^[35][44]. An in vitro study on the effect of sulforaphane on IAV replication in Madin-Darby canine kidney cells showed a decrease in replication ^[36][124]. The Human immunodeficiency virus (HIV), a human retrovirus belonging to the lentivirus family that attacks the body’s immune system, is responsible for causing acquired immunodeficiency syndrome (AIDS) once the immune system is compromised. There are two known types of HIV, which are deeply related, named HIV-1 and HIV-2. HIV-1 is the main cause of AIDS and responsible for the worldwide HIV pandemic; HIV-2, which differs in genomic structure and antigenicity, also causes the disease but with slower progression ^[50][135]. Different plants and phytochemicals exert activity against HIV. For example, gallic acid and ellagic acid possess anti-HIV activity and may act as post-entry inhibitors by affecting HIV protease activity ^[51][26]. Hydroxytyrosol also shows anti-HIV activity in vitro, inhibiting the viral integrase enzyme and fusion of the viral envelope with host cells ^[52][34]. Apigenin, as well as epicatechin, (-)-epicatechin-3-O-gallate (ECG), and epigallocatechin-3-gallate (EGCG), exhibit antiviral activities against HIV, presumably inhibiting HIV-1 protease enzyme and HIV reverse transcriptase ^[23][27][38,83]. Luteolin inhibits HIV-1 at the transcription level after the viral integration step ^[24][39], while myricetin inhibits the HIV-1 integrase (IN) and the HIV-1 reverse transcriptase ^[53][62]. The results from Ortega and co-workers ^[54][68] showed that the differences in the antiviral activity and the RT inhibition of the myricetin and myricetin-3-rhamnoside (myricitrin) together with the in silico analysis suggested that the glycosyl moiety could play a role in the entry of flavonoids into the cell and then, after enzymatic cleavage of the glycosyl moiety, the myricetin aglycone would ultimately be responsible for the anti-HIV activity. Quercetin has been investigated in vitro as an antiviral agent for HIV due to its ability to inhibit crucial enzymes such as reverse transcriptase, integrase, and protease ^[55][72]. Resveratrol also has effects on viral infection by inhibiting HIV-1 strain replication ^[24][39]. Several studies have reported that curcumin exhibits anti-HIV activity by directly targeting viral proteins. In particular, curcumin shows inhibition of the HIV-1 and HIV-2 proteases, the HIV integrase and the HIV trans-activator of transcription (tat), a viral transcription regulator ^[19][101]. Coumarin derivatives also exhibit their anti-HIV activity by targeting viral proteins inhibiting HIV protease, integrase, reverse transcriptase, and tat but also by inhibiting viral DNA replication ^[35][44]. Sulforaphane protects macrophages from HIV infection by mobilizing the transcription factor and antioxidant response regulator nuclear factor E2-related factor 2 (Nrf2) ^[36][124]. Enterovirus 71 (EV71), which belongs to the genus Enterovirus, is the major pathogen of hand, foot, and mouth disease (HFMD) that often occurs in children under the age of 10 ^[56][27]. Gallic acid was found to exert strong anti-EV71 activity ^[56][27], while isorhamnetin inhibits EV71 RNA replication and protein synthesis ^[23][38]. It has been demonstrated that resveratrol inhibits EV71 replication ^[57][106]. Antiviral activity of epigallocatechin-3-gallate (EGCG) and coumarin and its derivatives has been observed against EV71 ^[27][35][44,83]. Chebulagic acid and Punicalagin exhibit antiviral activity both in vitro and in vivo against EV71 ^[58][59][112,115].

3. Natural Phytocompounds with Potential to Inhibit the Coronavirus SARS-CoV-2 According to In Silico Approaches

In silico molecular modeling uses the crystal structure of viral targets from SARS-CoV-2, or targets of related viruses such as SARS-CoV, docked with molecules under study to establish the binding strength of competitive inhibitors of the target that is expressed in terms of Kcal/mol. The binding strength of a test compound is compared with that of a known inhibitor or control. The in silico docking is no indication that a compound will absolutely inhibit a target, but it offers a good start for engaging in drug design ^[60][146]. RIn this researchersview, we have performed an analysis of studies carried out with an in silico approach of possible interactions of the same molecules analyzed above for their antiviral activity with the major targets of SARS-CoV-2, like the main protease (M^Pro), RdRp, membrane (M) glycoprotein, envelope (E) glycoprotein, spike (S) glycoprotein, PL^pro, non-structural proteins (NSP) and ACE2 receptor. In particular, for those molecules that have shown activity against SARS-CoV, it has been postulated that they could potentially act against SARS-CoV-2, taking into account that SARS-CoV and SARS-CoV-2 have a high sequence similarity (79.5%) ^[61][147]. For example, Alrasheid and co-workers ^[62][30] have conducted an in silico study by using Molecular Operating Environment (MOE) drug discovery software platform in order to evaluate the antiviral activity targeting SARS-CoV-2 of twenty-one compounds from medicinal plants. They have found that gallic acid, quercetin, and capsaicin were among the best compounds interacting with SARS-CoV-2 M^pro, with rank scores ranging from −17.45 to −13.90 Kcal/mol. Moreover, a screening via molecular docking to test the binding affinity of various selected bioactive compounds of honey and propolis as inhibitors against the SARS-CoV-2 M^pro and RdRp revealed that ellagic acid and kaempferol have the strongest interaction with RdRp, instead the previously cited compounds and quercetin with the SARS-CoV-2 M^{pro [63]}[32]. In another study, in which one hundred secondary metabolites from Aframomum melegueta, a Zingiberaceae family plant spice widely spread in Africa, have been computationally evaluated for inhibition of SARS-CoV-2 targets, it is reported that quercetin could be a potential SARS-CoV-2 2′-O-methyltransferase (NSP16) enzyme inhibitor while apigenin could inhibit SARS-CoV-2 main protease (M^pro) ^[64][41]. In a 2020 study published as a preprint, kaempferol, quercetin, naringenin, oleuropein, curcumin, catechin, and epicatechin-gallate appeared to have the best potential to act as SARS-CoV-2 M^pro inhibitors with binding energies values obtained from the docking of main protease of −8.58, −8.47, −7.89, −7.31, −7.05, −7.24, and −6.67 Kcal/mol, respectively ^[65][57]. Although the preprint was not published elsewhere to date, other studies confirm the results. Bilginer et al. ^[66][58] also showed similar data for kaempferol, quercetin, and oleuropein; instead, Halder et al. ^[67][81] for curcumin and catechin. Moreover, Mukheriee et al. ^[68][59] confirmed a possible binding of epicatechin-gallate, quercetin, and kaempferol to the M^pro active site, with binding energies values of −8.5, −7.6, and −7.3 Kcal/mol respectively, and underline as epigallocatechin-gallate shows the best affinity in comparison to the re-docked inhibitor binding energy (−7.7 Kcal/mol). In another in silico study, published as a preprint, luteolin was also screened against SARS-CoV-2 main protease (M^pro), showing binding energy values of −8.35 Kcal/mol ^[35][44]; also, in this case, further studies corroborate this hypothesis ^[69][45]. Based on the inhibitory activity against SARS-CoV 3CL^pro exerted by pectolinarin, rhoifolin, and herbacetin ^[70][71][47,49] performed an in silico docking study to deduce their binding mode and binding affinity with SARS-CoV-2 3CL^pro and found that the affinity of rhoifolin became weakened while the efficiency of herbacetin and pectolinarin was still promising. Vicidomini et al. ^[72][54] found that for several flavonols and flavonol glucosides isolated from Opuntia ficus-indica, a dicotyledonous angiosperm cactaceous plant widespread worldwide in tropical and subtropical regions, including isorhamnetin, possible affinities for SARS-CoV-2 M^pro, with binding energies lower than −7.0 Kcal/mol. Myricetin-3-rhamnoside (myricitrin) was found to have a strong binding with the active site residues of the SARS-CoV-2 M^pro with a binding score of −8.9 Kcal/mol. In particular, its active site residues Tyr54, Phe140, Gly143, His163, and Glu166 are involved in the formation of 6 hydrogen bonds with hydroxyl and carboxyl oxygen of ligand while His41 and Cys145 form Pi-alkyl and Pi-sulfur interactions with the aromatic scaffold of ligand, respectively. Other additional interactions, such as van der Waals forces, stabilize the binding of the ligand to the enzyme ^[73][69]. In a study in which thirty-eight flavonoids have been tested by molecular docking against the active site of the SARS-CoV-2 M^pro, it has been demonstrated that natural aglycone flavonoids possess higher docking energies than flavonoids with sugar moieties. As a matter of fact, in these docking experiments, the aglycone myricetin displays a binding score of −7.4 Kcal/mol while myricetin-3-rhamnoside (myricitrin), also in this study, displays a binding score of −8.9 Kcal/mol. Similarly, quercetin-3-O-rhamnoside (quercitrin), quercetin-3-O-rutinoside (Rutin), quercetin-3-beta-galactoside, and quercetin 3-O-β-glucuronide exhibit lower binding energies than their aglycone quercetin (binding score of −9.7, −9.2, −8.4, −8.1, −7.5 Kcal/mol, respectively) ^[74][65]. Catechins and some of their derivatives showed a strong affinity for SARS-CoV-2 M^pro because they could form a stable ligand-receptor complex, and this leads to the hypothesis of a potential antiviral activity ^[75][79]. An in silico analysis has depicted that ellagic acid, epicatechin, and capsaicin exhibit significant binding and interaction with the most vital active site residue Cys145 of SARS-CoV-2 M^pro and can inhibit its activity in an extremely effective manner ^[76][82]. Ghosh et al. ^[77][80] selected eight polyphenols from green tea and elucidated the binding affinities and binding modes between these polyphenols and SARS-CoV-2 M^pro using molecular docking studies. All eight polyphenols exhibit good binding affinity, including catechin and (-)-epigallocatechin (EGC), but epigallocatechin-3-gallate (EGCG) and (-)-epicatechin-3-O-gallate (ECG) interact strongly with one or both the catalytic residues His41 and Cys145 of SARS-CoV-2 M^pro with a binding score of −7.6 and −8.2 Kcal/mol. Another study aimed to evaluate bioactive compounds found in plants to interact and hopefully inhibit SARS-CoV-2 M^pro and spike (S) glycoprotein using a molecular docking approach reported that pectolinarin, epigallocatechin-3-gallate (EGCG) and rhoifolin are among potential candidates to bind spike (S) glycoprotein surface, while rhoifolin and pectolinarin are among the candidates to become drugs targeting SARS-CoV-2 M^{pro [78]}[48]. The role of tea polyphenols in prophylaxis and treatment of COVID-19 has been investigated by Mhatre et al. ^[79][84] in a preliminary in silico study in which molecular docking interactions of two tea polyphenols with some of the possible binding sites of SARS-CoV-2 were performed. In particular, the receptors 3CL^pro, RdRp, PL^pro, spike (S) glycoprotein RBD, and ACE2 receptor with spike (S) glycoprotein RBD were docked against epigallocatechin-3-gallate (EGCG) from green tea. For all the receptors studied, EGCG exhibited good docking scores with various types of interactions, and except for RdRp, the binding energies were consistently better than −8.0 Kcal/mol. Mostafa et al. ^[80][88] have conducted an in silico molecular docking experiment on the effects of some natural antiviral phytoconstituents on the crystal structure of SARS-CoV-2 M^pro. Many of the docked compounds revealed good binding affinity, with theaflavin showing a good binding score (−8.5 Kcal/mol). The bioactive phenolic phytocompound 6-gingerol, found in the fresh ginger rhizome, was also tested against different viral proteins by using the molecular docking technique in order to evaluate its interaction with SARS-CoV-2 targets. This study suggests a high binding affinity and interaction with surface regions of multiple targets of COVID-19, including viral proteases, RdRp, and spike (S) glycoprotein ^[81][99]. Du et al. ^[82][114] investigated in vitro and in silico the potent antiviral activity of chebulagic acid and punicalagin against SARS-CoV-2 viral replication. In silico docking of chebulagic acid and punicalagin to SARS-CoV-2 M^pro was performed in order to search potential allosteric binding sites, and both chebulagic acid and punicalagin may interact with the cleft between domain II and domain III within SARS-CoV-2 M^pro with stable binding free energy. It is reported that SARS-CoV-2 can bind to Toll-like receptor 4 (TLR-4), which would eventually lead to α-synuclein aggregation in neurons and stimulation of neurodegeneration pathways. Oleuropein was investigated against the SARS-CoV-2 target (main protease 3CL^pro), TLR-4, and Prolyl Oligopeptidases (POP), to explore oleuropein potency against the neurological complications associated with COVID-19. POP inhibition could reduce aggregation of alpha-synuclein in neurons, hence, investigating anti-COVID agent(s) against TLR-4 and POP might help in alleviating neurological complications linked with COVID-19. Docking analyses showed a binding score of −7.8, −8.3, and −8.5 Kcal/mol for oleuropein-3CL^pro, oleuropein-TLR4, and oleuropein-POP interactions, respectively, suggesting oleuropein as a potential candidate that can target SARS-CoV-2 and alleviate neurological manifestations associated with it ^[83][119]. In a recent study, the antiviral potential of 50 natural coumarin phytochemicals isolated from plants was investigated in order to identify the binding interactions of these phytochemicals against the coronavirus 3CL^pro by molecular modeling approaches. It was found that glycycoumarin, Inophyllum P, mesuol, and oxypeucedanin hydrate displayed the highest binding affinity with the best negative energy scores and interacted with one or both of the catalytic residues (His41 and Cys145) of 3CL^pro through hydrophilic and hydrophobic bonding ^[84][121]. Shekh et al. ^[85][127], using virtual screening methods, have analyzed allicin to assess its ability to covalently modify cysteine residues of SARS-CoV-2 M^pro. The results suggest that allicin may induce dual S-thioallylation of Cys145 and Cys85/Cys156 residues of SARS-CoV-2 M^pro and may be useful in attenuating the coronavirus infection. It is known that essential oils have anti-inflammatory, immunomodulatory, bronchodilatory, and antiviral properties and, for this reason, are being proposed to have activity against SARS-CoV-2 virus. Moreover, owing to their lipophilic nature, essential oils are advocated to penetrate viral membranes easily, leading to membrane disruption. Molecular docking techniques were applied to screen the anti-SARS-CoV-2 efficacies of eugenol, menthol, and carvacrol, major components of essential oils, against various protein targets of SARS-CoV-2. Results revealed that these compounds have binding affinities towards SARS-CoV-2 spike (S) glycoprotein, SARS-CoV-2 M^pro, RdRp, and human ACE2 proteins, respectively ^[31][128]. Paolacci and co-workers ^[86][36] investigated the potential effects of hydroxytyrosol in combination with alpha-cyclodextrin, another naturally occurring compound, on SARS-CoV-2 entry into human cells. Bioinformatic docking studies showed that hydroxytyrosol is captured in the hydrophobic cavity of alpha-cyclodextrin, and the resulting complex could interact with the SARS-CoV-2 spike (S) glycoprotein and its host cell receptor ACE2. Even if these bindings do not occur at the spike-ACE2 interface, they may induce structural perturbations that could potentially influence the endocytosis process. Sekiou et al. ^[87][51] have screened in silico the interaction between the main protease SARS-CoV-2 M^pro active site with natural compounds, including cirsimaritin and quercetin, and found that both have a better binding affinity to this target. Molecular docking was used by Pandey et al. ^[24][39] to describe the binding ability of naturally occurring phytochemicals, including kaempferol and quercetin, with SARS-CoV-2 spike (S) glycoprotein and found that these natural compounds are capable of binding to either the S1 or S2 domains of the SARS-CoV-2 S, most probably preventing it from binding to the ACE2 receptor or internalization during fusion. High-end molecular docking analysis and MD simulation were performed to characterize the binding affinity of natural and synthetic anti-viral compounds with the SARS-CoV-2 structural proteins and to analyze the stability of drug-protein interactions. Results identified rutin as a potent inhibitor of SARS-CoV-2 envelope (E) glycoprotein, while caffeic acid and ferulic acid were found to inhibit SARS-CoV-2 membrane (M) glycoprotein ^[88][22]. Orfali et al. ^[89][33] screened all available SARS-CoV-2 molecular targets using a multistep in silico protocol to find out the most probable one that mediates the proved in vitro anti-SARS-CoV-2 activity of sinapic acid and found that the viral envelope protein (E-protein) was suggested as the most probable hit for sinapic acid. A further in-depth molecular dynamic simulation-based investigation was performed to precisely explore the binding path and mode of sinapic acid with E-protein. The interactions of the two crucial proteins NSP9 and NSP10 of COVID-19 have been investigated with potential antiviral compounds from Moringa oleifera using molecular docking and dynamic methods. The results revealed that all selected ligands form stable complexes with the targeted proteins and showed the highest binding affinity values of apigenin (−7.1 Kcal/mol) for NSP10 and ellagic acid (−7.1 Kcal/mol) for NSP9 ^[90][40]. Several viral proteins have been found to possess deubiquitinating activity, and they can be used to antagonize or modulate the antiviral immune signaling pathway. Along with protease activity, SARS-CoV-2 PL^pro possesses deubiquitinating activity. Naphthalene-based inhibitors, such as the well-investigated GRL-0617 compound, have been shown to possess dual effects, inhibiting both protease and deubiquitinating activity of the PL^pro. Pitsillou et al. ^[91][52] investigated the binding characteristics of the same dietary compounds, including rutin and cyanidin-3-O-glucoside, to the PL^pro and evaluated the deubiquitinating activity. Analyses of molecular docking highlighted the relatively high affinity of GRL-0617 and dietary compounds. A study was conducted on the most abundant pomegranate (Punica granatum) peel extract constituents, including punicalagin, with the aim of exploring their anti SARS-CoV-2 properties using in silico tools. The protein targets used in this study were: SARS-CoV-2 spike (S) glycoprotein, ACE2, furin, and TMPRSS2. Molecular docking results showed that all the analyzed ligands interacted through hydrogen bonds with amino acid residues at the S glycoprotein predicted druggable site, and punicalagin presented free binding energy of −7.312 Kcal/mol, with the ligand-protein complex stabilized through three hydrogen bonds (Asn343, Asn370, and Ser371). All tested pomegranate peel extract constituents showed significant binding affinity at the ACE2 predicted druggable site, with all ligand–protein complexes stabilized through hydrogen bonds, and punicalagin presented free binding energy of −7.144 Kcal/mol. Amino acid residue Lys441 was found to be important for the stabilization of punicalagin–ACE2 complex. Moreover, all pomegranate peel extract constituents formed stable complexes with furin, with punicalagin-free binding energy of −9.385 Kcal/mol. Punicalagin showed intensive interactions with TMPRSS2 amino acid residues at the predicted binding site with binding energy values of −7.358 Kcal/mol and with four hydrogen bonds. Polar interactions with Asn97 and Arg405 residues were essential for the stabilization of TMPRSS2 complexes with punicalagin ^[92][117]. Özdemir et al. ^[93][122] carried out molecular docking studies on five different proteins (Spike S1-subunit, NSP5, NSP12, NSP15, and NSP16) of the SARS-CoV-2 and two proteins (ACE2 and Vitamin K epoxide reductase complex subunit 1-VKORC1, a key enzyme in recycling reduced vitamin K, that plays an essential role in γ-carboxylation of vitamin K-dependent coagulation factors) of human and found that the best binding scores for 17 coumarins were determined for NSP12, with the highest score (−10.01 Kcal/mol) showed by 2-morpholinoethan-1-amine substituted coumarin. In conclusion, it is necessary to emphasize that the potential strength of molecular simulations offers the possibility of carrying out both large-scale screening and targeted docking simulations quickly and with limited costs, but on the other hand, it is also necessary to apply rigorous procedures and carefully evaluate the results obtained. The COVID-19 pandemic has prompted many research groups to carry out docking simulations to provide useful information to the scientific community as soon as possible and to identify possible solutions to the emergency as quickly as possible. However, reswearchers found several studies published only in pre-print form, and therefore without having received an evaluation by experts; others were published in peer-reviewed journals but very few days after submission to the journal. Although peer-reviewed, the very short evaluation time makes the quality of the review at least questionable. The result is that in several literature studies, as either full articles or pre-prints, there are poor protocols of docking simulations, weak or erroneous interpretations of the results, use of protein structures that are not experimental results but obtained by insufficiently described modeling procedures, docking with portions of proteins without critical evaluation of the real utility. Here, we chose to quote only in silico simulation studies that meet sufficient quality criteria in the simulations performed and in the presentation and discussion of the results obtained.