key: cord-0872218-6sezewft
authors: Zígolo, María Antonela; Goytia, Matías Rivero; Poma, Hugo Ramiro; Rajal, Verónica Beatriz; Irazusta, Verónica Patricia
title: Virtual screening of plant-derived compounds against SARS-CoV-2 viral proteins using computational tools
date: 2021-03-17
journal: Sci Total Environ
DOI: 10.1016/j.scitotenv.2021.146400
sha: b9a06ac95c980b1b1d71cdf6a98cd10e869b4518
doc_id: 872218
cord_uid: 6sezewft

The new SARS-CoV-2, responsible for the COVID-19 pandemic, has been threatening public health worldwide for more than a year. The aim of this work was to evaluate compounds of natural origin, mainly from medicinal plants, as potential SARS-CoV-2 inhibitors through docking studies. The viral spike (S) glycoprotein and the main protease Mpro, involved in the recognition of virus by host cells and in viral replication, respectively, were the main molecular targets in this study. Molecular docking was performed using AutoDock, which allowed us to select the plant actives with the highest affinity towards the viral targets and to identify the interaction molecular sites with the SARS-CoV2 proteins. The best energy binding values for S protein were, in kcal/mol: −19.22 for glycyrrhizin, −17.84 for gitoxin, −12.05 for dicumarol, −10.75 for diosgenin, and − 8.12 for delphinidin. For Mpro were, in kcal/mol: −9.36 for spirostan, −8.75 for N-(3-acetylglycyrrhetinoyl)-2-amino-propanol, −8.41 for α-amyrin, −8.35 for oleanane, −8.11 for taraxasterol, and − 8.03 for glycyrrhetinic acid. In addition, the synthetic drugs umifenovir, chloroquine, and hydroxychloroquine were used as controls for S protein, while atazanavir and nelfinavir were used for Mpro. Key hydrogen bonds and hydrophobic interactions between natural compounds and the respective viral proteins were identified, allowing us to explain the great affinity obtained in those compounds with the lowest binding energies. These results suggest that these natural compounds could potentially be useful as drugs to be experimentally evaluated against COVID-19.

cases reported; however, death rates continue to increase with over 88,000 new deaths reported last week (World Health Organization, 2020) . As of February 9, there have been almost 105 million cases and over 2.3 million deaths reported globally since the beginning of the pandemic.

These data demonstrate the urgent need to find quick solutions to the problem.

To date, very vaccines against SARS-CoV-2 have been developed showing good percentages of effectiveness in clinical trials; however, the world population has not been massively vaccinated and the period of immunity that these vaccines will generate is yet unknown (Malik et al., 2020; Detoc et al., 2020; Kaur et al., 2020) . On the other hand, it is important to mention that therapies with monoclonal antibodies proved to be effective as modulators of the immune response against SARS-CoV-2 (Jahanshahlu et al., 2020; Pecetta et al., 2020; Zhao, 2020) . Some monoclonal antibodies such as infliximab are already in the stage of clinical trials in patients with COVID-19 (Clinical Trials, 2020) . Also, many synthetic antiviral molecules, already existing in the market were evaluated at different levels, many of which reached the stage of clinical trials. However, none of them showed to be effective enough against SARS-CoV-2 (Clinical Trials, 2020; Frediansyah et al., 2021; Whang et al., 2020) . In this context, the search for new therapeutic agents is imperative as another possible alternative for the treatment of COVID-19; indeed, great effort and important resources are nowadays being dedicated worldwide to achieve that goal.

Antivirals are among the compounds used to prevent the development of viruses, useful in restricting viral infection. Of special interest are the secondary metabolites of plant origin since many of them not only act as powerful antivirals but also have the ability to stimulate the immune system, which is essential in a response to SARS-CoV-2 infection (Luengo-Fernandez, 2008; Cinatl et al., 2003; Rawat et al., 2020) .

Compounds known for their effectiveness against other viruses are being investigated in search J o u r n a l P r e -p r o o f of a solution to the rapid and lethal advance of the new SARS-CoV-2. Examples of that are some protease inhibitors used for treatment of patients with Human Immunodeficiency Virus infection and acquired immune deficiency syndrome (HIV/AIDS) or neuraminidase inhibitors or nucleotide analogues like Remdesivir (GS-5734™) used for other pathogenic CoVs (Khaerunnisa et al., 2020) . Another interesting molecule is chloroquine/hydroxychloroquine, an acidotropic amine form of quinine, extensively used against malaria and previously described as a potent inhibitor of most coronaviruses in vitro, including SARS-CoV-1 (Devaux et al., 2020; White et al., 2014) . It is interesting to note that quinone is a metabolite produced by cinchona trees (Cinchona officinalis), native from Peru. Plants are great producers of secondary metabolites and as such, especially those from poorly studied regions, they captured the interest of many researchers trying to develop new medicines (Yang et al., 2018) . Antiviral bioactivities have been described in numerous medicinal plants and associated with compounds like flavonoids, heterosides, terpenes, and triterpenes, organic acids, alkaloids, saponins, and quaternary ammonium salts, among others (Cecil et al., 2011; Khaerunnisa et al., 2020; Roy et al., 2006; Wang et al., 2017) . Molecular docking is a useful tool in structural molecular biology to evaluate the possibility of interactions and to predict binding modes of ligand-proteins from knowledge of the threedimensional structure. The search is conducted in a way that all the candidate dockings are ranked through the Gibbs free energy (negative values indicate spontaneous chemical reactions; the greater the magnitude, the higher is the affinity of the ligand-macromolecule interaction), allowing us to make hypotheses or decisions about how, or which, ligands inhibit the target. This tool can play a key role in computer-assisted drug design or in the selection of the best candidates from a library of compounds with a certain activity (Morris et al., 2008) . Analyzing J o u r n a l P r e -p r o o f these virtual protein-ligand interactions becomes critical when rushing in the search of metabolites that can cure or ease symptoms of diseases when time and resources for extensive experimental work are scarce. Considering that in the last months many SARS-CoV-2 macromolecules have been crystallized, more effective simulations may be performed cross-out if model proteins from other closely related viruses have to be used (Jin et al., 2020; Wrapp et al., 2020; Zhang et al., 2020) .

In the initial stage of the SARS-CoV-2 replication cycle, the attachment of the virion to the host cell is initiated by interactions between the S protein (spike) and human angiotensin-converting enzyme (ACE2 receptor) (Hoffmann et al., 2020) . The S protein-receptor interaction is the primary determinant for a coronavirus to infect a host species and governs the tissue tropism of the virus. Structurally, the coronavirus spicule is a trimer arranged asymmetrically. Each monomeric unit has three segments: an extensive ectodomain, a single-pass transmembrane segment, and a small intracellular tail (Walls et al., 2016) . The ectodomain contains two subunits; the S1 subunit allows binding to the ACE2 receptor while the S2 allows the fusion of viral and host membranes. Two major domains in S1, N-terminal domain (S1-NTD) and Cterminal domain (S1-CTD), have been identified. The S1-CTD, also known as the receptorbinding domain (RBD), is responsible for recognizing protein receptors ACE2 (Liu et al., 2015) . The next step in the infection cycle is the translation of the replicase gene from the viral genome (positive polarity RNA). Replicase gene encodes two overlapping polyproteins, pp1a and pp1ab, required for viral replication and transcription (Baranov et al., 2005) . Polyproteins pp1a and pp1ab contain the non-structural proteins (nsps) 1-11 and 1-16, respectively. These polyproteins are subsequently cleaved into the individual nsps (Ziebuhr et al., 2000) . SARS-CoV-2 encodes two proteases that cleave the polyproteins; the papain-like proteases, PLpro (encoded by nsp3) J o u r n a l P r e -p r o o f and chymotrypsin-like cysteine protease known as 3C-like protease (3CL pro ). The 3CL pro enzyme, also called main protease (M pro ), hydrolyses the polyprotein at no less than 11 conserved sites between nsp4 and nsp16, starting with the autolytic cleavage of this enzyme itself (Hegyi et al., 2020) .

The aims of this work were to evaluate by molecular docking selected compounds obtained from natural sources for antiviral activity against SARS-CoV-2. Synthetic drugs that are currently being evaluated for the treatment of COVID-19 were used as controls. Two viral proteins with different and essential functions were tested: the SARS-CoV-2 main protease (M pro ) and the S protein (spike) in its prefusion conformation .

SARS-CoV-2 M pro (PDB ID: 6LU7) and S protein (PDB ID: 6VSB) were used in this study.

Three-dimensional structure are shown in Table 1 .

Secondary metabolites with different structures were analyzed to select the most promising to assess a possible relationship between structure and activity. Many of them already had proven antiviral activity against other viruses (Alonso J. and Desmarchelier, 2015; Cecil et al., 2011; Enkhtaivan et al., 2018; Shin et al., 2015) . 

Docking calculations were carried out with Autodock 4.2 program, which uses Lamarckian genetic algorithm (LGA) for calculations (Morris et al., 2009) . Number of genetic algorithm (GA) runs was set to 200 for each case analyzed and binding affinities were determined and reported in kcal/mol. The Autodock 4.2 program was applied considering all rotatable bonds for ligands and the whole protein as a rigid structure. For the location and extent of the 3D area the search space was defined by specifying a center, the number of points at each dimension and the points between spaces to focus the search space in the enzyme active site or putative binding region. Prior to docking M pro , the inhibitor N3 and water molecules were removed from the protein structure. Polar hydrogen atoms were added, and Gasteiger atom charges were assigned to protein atoms. Other parameters were set to default values.

The 200 conformers found for each compound with the Autodock program were grouped in clusters that were ordered according to a ranking, which was determined by the stability of the enzyme-ligand complexes within each cluster. Clustering of the 200 conformers was done according to the similarity with the conformation adopted inside the enzyme. The criterion used to evaluate such similarity was the residue mean quadratic square root deviation (RMSD). The residues were obtained through the difference between the atom coordinates of a given conformer respect to the cluster to which the most stable conformer belongs. Each cluster grouped conformers with RMSD lower or equal to 2.0 Å.

For visualization of protein-ligand complexes the software Visual Molecular Dynamics 1.9.1 (VMD) was used (Theoretical and Computational Biophysics Group, University of Illinois).

The frequency distribution (number in cluster) of plant-derived compounds were first plotted vs.

Gibbs free energy of conformers for a specific target or protein. Then, the best compounds for each protein were selected considering the following criteria for each compound: i) clusters with highest frequencies, ii) clusters with lowest Gibbs free energies, iii) conformational dispersion associated with the number of clusters (the higher the number of clusters was, there was more conformational dispersion for the interaction studied), iv) intra-cluster dispersion given by the variation around the mean binding energy, and v) intra-cluster dispersion given by the percentile of energy values close to the minimum binding energy (more negative free energy). Statistical J o u r n a l P r e -p r o o f software R (https://cran.r-project.org/) with the Rstudio IDE (https://rstudio.com/) was used for processing.

Different clusters were ordered according to the energy of the most stable conformer within each group: the ranking was ordered according to the increasing stability of the lowest energy conformers in each cluster.

From the 200 docking runs for all the evaluated compounds, both from natural origin and the synthetic controls, the binding energies for each conformer were obtained which, in turn, were grouped into clusters. The binding energy for each conformer was plotted as a function of the frequencies (number of conformers in each cluster) for the M pro protein (Figure 1, A and B; S1A) and for S glycoprotein (Figure 2 ; S2).

The interaction of all the natural compounds evaluated with each viral protein showed, in general, a low number of clusters although they were highly frequent, suggesting that there was a high probability for them to happen. Conversely, the synthetic antiviral nelfinavir and atazanavir used as controls for the M pro protein and umifenovir, chloroquine, and hydroxychloroquine for the spike protein, showed broad conformational dispersion or high number of clusters (Figures 1 B and 2, S1A).

Each compound was represented by the most populated clusters, which in some cases were also those with the lowest binding energies. Then, statistical criteria were applied to select the best clusters when more than one was considered for one compound, and the binding energies associated with their interactions with the viral proteins. Three different statistics were calculated: a) the standard deviation of the Gibbs free energy of all the cluster´s conformers (SD J o u r n a l P r e -p r o o f binding energy), b) the variation coefficient of the Gibbs free energy of all the cluster´s conformers (CV), and c) the intra-cluster dispersion given by the 1% percentile of the Gibbs free energy close to the lowest binding energy (Lowest 1% perc). These metrics were applied in order to assign to each compound the Gibbs free energy corresponding to the interaction with the viral proteins based on the following criteria: i) the Mean binding energy was assigned when CV was lower than or equal to 10% and the value of the Lowest 1% perc was not higher than 60%; ii) the Lowest binding energy was assigned when the Lowest 1% perc reached at least 60%, despite the CV value was lower than or equal to 10% (Table 2 for M pro and Table 3 for S protein).

When two or more clusters were initially selected for the same compound (for having similar frequencies and binding energies), the following criteria were applied: i) if the clusters compared had Lowest 1% perc value equal to or greater than 60%, the one with the highest Lowest 1% perc value was chosen, and ii) when the clusters compared had a Lowest 1% perc value lower than 60%, the CVs were compared and that with the lowest CV value was chosen.

Following the statistical criteria for the docking results, the binding energies for the interactions of each compound with the corresponding protein were selected (Table 4) .

When M pro protein interacted with natural compounds, similar values than those observed with the synthetic antivirals atazanavir and nelfinavir, were found. However, for synthetics control compounds greater dispersion and lower frequencies were observed (Figure 1 A and B ; S1A).

The strongest affinities for spirostan and for N-(3-acetylglycyrrhetinoyl)-2-amino-propanol, with ΔG of -9.36 and -8.75 kcal/mol, were obtained, respectively. Other compounds showed interesting affinities for M pro , for example: α-amyrin, oleanane, taraxasterol, and glycyrrhetinic acid (Table 4) .

On the other hand, the plant-derived natural compounds evaluated for S protein, yielded very J o u r n a l P r e -p r o o f high affinity values compared to controls, except for aloin (-5.94 kcal/mol). The highest affinities were observed for glycyrrhizin, gitoxin, dicumarol, diosgenin, and delphinidin, which also showed low dispersion and high frequencies (Figure 2 ; S2). The most interesting results were observed for glycyrrhizin and gitoxin (-19.22 and -17.84 kcal/mol, respectively), which presented more than twice the affinity for the viral S protein compared to the synthetic antivirals umifenovir, chloroquine, and hydroxychloroquine (-7.47, -7.25, and -5.30 kcal/mol, respectively) ( Table 4) .

Key interactions of possible candidate drugs with M pro were analyzed: α-amyrin, spirostan, taraxasterol, oleanane, lapachol, glycyrrhetinic acid, and N-(3-acetylglycyrrhetinoyl)-2-aminopropanol (Table 5 ). All the natural and semisynthetic active compounds bound to the enzyme active site, as well as the controls nelfinavir and atazanavir. Only interactions of nelfinavir with M pro with lower binding energy were shown (Table 5 ). The semi-synthetic compound N-(3acetylglycyrrhetinoyl) -2-amino-propanol presented five hydrogen bond interactions, four with the residues of the active site and one with the tyrosine residue 54. It also showed hydrophobic interactions with six residues of the protease in the vicinity of the catalytic site. All these interactions increased the binding affinity with respect to the control compound, which presented only two hydrogen bond interactions with the active site residues and hydrophobic interactions with other four protease residues. Glycyrrhetinic acid, a precursor of the above-mentioned semisynthetic compound, had four hydrogen bond interactions with catalytic residues and hydrophobic interactions with leucine residue 4.

It is important to note that five of the six candidate compounds (Table 5) In addition to glycyrrhizin and gitoxin, the interaction of dicumarol with the spike protein showed one of the most favorable Gibbs free energies (-12.05 kcal/mol) (Table 6 ). Glycyrrhizin, with the highest affinity for the glycoprotein, interacted with a great number of residues from A and C chains. There was a total of four hydrogen bonds with residues of the A chain (PHE855, VAL976 and ASP979), and another two with the B chain (ASP571 and ASP574). Hydrophobic interactions occurred with the remaining of the 14 residues.

Plants have a long evolutionary history of developing resistance against viruses. This adaptation, together with the ability to produce secondary metabolites, has generated a wide range of possible sources for antiviral drugs. Herbal infusions and essential oils from plants have been used ancestrally in different civilizations. Their low toxicity and minimal adverse effects at the doses used have long been known, representing an important advantage as compared with new synthetic drugs, which need to be evaluated for safety through strict clinical trials (Abbas et al., 2016; Casais et al., 1997; Horden et al., 2016) . Other advantages of natural vs. synthetic (pure) compounds are that: i) they are mixtures of active compounds named phytocomplexes, formed J o u r n a l P r e -p r o o f by isomers and analogue compounds, with the capacity for synergistic activity, and ii) the mixture usually contains compounds that help mitigate the toxicity of the main active drugs.

Also, microorganisms are less likely to develop resistance to a phytocomplex than to a unique compound in a synthetic drug (Alonso, 2008; Bruneton, 2001; Camponovo et al., 1955) . To simplify the analyses in this study, only the main active compounds were used since the chemical structures are known.

Liu and coworkers (2018) successfully crystallized the SARS-CoV-2 main M pro protease, which is a potential drug target for the inhibition of SARS-CoV-2 as it plays an essential role in viral replication and transcription (Liu et al., 2015) . From the docking of M pro with natural compounds, the most favorable binding energies were obtained for compounds with steroid and pentacyclic triterpene structures. However, it is important to note that compounds with different structures extracted from plants were evaluated. Thereby alkaloids (berberine, theophylline, sauroxin), terpenoids (limonene, geranial), phenols (lapachol, lapachenol), alcohols (menthol), steroids (spirostan), flavonoids (rutin), glycosides (sambunigrin), pentacyclic triterpenes (amirin, taraxasterol, oleanaane, glycyrrhetinic acid, N-(3-acetylglycyrrhetinoyl) -2-amino-propanol) and macrolides (maytansine), were assayed.

Berberine is an isoquinoline alkaloid extracted from plants of the Berberis family (Coptis chinensis Berberis buxifolia, Berberis vulgaris, Argemone mexicana, etc.). It has antimalarial, antibacterial, and antiviral effects, among other biological activities (Wang et al., 2017) .

Berberine can also inhibit the replication of influenza A in vitro in several different cell types and with two influenza A (H1N1) strains (Cecil et al., 2011) . As in this study, other researchers found that berberine displayed binding energy of -7.91 kcal/mol and had the ability to interact with the protease through the amino acid ASP187 producing different conformational changes in J o u r n a l P r e -p r o o f the target enzyme (Agrawal et al., 2020) .

In particular, the structural family of pentacyclic triterpenes, with high affinity values for viral protease, could be suggested as possible candidates for the treatment of COVID-19. Taraxasterol is a pentacyclic triterpene only found in a few plants like arnica, chicory, and dandelion, but with remarkably interesting biological activity (Krajcovicova et al., 2004; Zígolo et al., 2018) . (Cecil et al., 2011; Enkhtaivan et al., 2018; Shin et al., 2015) .

In the molecular docking for M pro protein of SARS-Cov-2, the semi-synthetic compound N-(3-J o u r n a l P r e -p r o o f acetylglycyrrhetinoyl)-2-amino-propanol performed better than glycyrrhizic acid. The former, obtained by biocatalytic reaction of aminolysis from glycyrrhetinic acid (derived from the plant Glycyrrhiza glabra L), already showed powerful antiviral activity against the herpes virus types I and II, including strains resistant to acyclovir (De Clercq et al., 2000; Zígolo et al., 2018) . In addition, Zígolo and coworkers (2018) showed low toxicity of the N -(3-acetylglycyrrhetinoyl)-2-amino-propanol in Vero cell cultures, making this compound a promising candidate for COVID-19 treatment 27 . In a molecular docking study of natural products that are commonly used in Indian and other cuisines against SARS-CoV-2 main protease (PDB: 6LU7 and PDB: 6Y2E), it was shown that several of these compounds had high binding affinity against M pro , and they were comparable with a known anti-HIV drug, saquinavir. In our study, we report plant-derived active compounds with much higher binding affinities than other natural compounds reported as possible antivirals for SARS-CoV-2. Glycyrrhizin (binding energy of -8.40 kcal/mol) was one of the best potential candidates assessed. It is remarkable that a diammonium salt of glycyrrhizin was recommended for the cure of COVID-19 in China. The binding affinity of the positive control, saquinavir, was -7.90 kcal/mol (Sampangi-Ramaiah et al., 2020) .

In a recent study, several compounds, mainly flavonoids (kaempferol, quercetin, luteolin-7glucoside, demethoxycurcumin, naringenin, apigenin-7-glucoside, oleuropein, curcumin, catechin, epicatechingallate, zingerol, gingerol, and allicin) and the synthetic compounds nelfinavir and lopinavir were analyzed as potential inhibitors of the SARS-CoV-2 M pro (PDB: 6LU7) by docking studies (Khaerunnisa et al., 2020) . In our hands, the binding energy found (-8.25 kcal/mol) for the interaction of nelfinavir with M pro protein was similar to that reported by Khaerunnisa and coworkers (2020) (-8.37 kcal/mol). Interestingly, despite the similarity in the ΔG obtained in both studies, differences in the hydrogen bonds formed were observed J o u r n a l P r e -p r o o f (Khaerunnisa et al., 2020) . While residues Glu166, Gln189, and Gln192 were used according to Khaerunnisa and coworkers (2020) , we only found Glu166 and His172 in the interaction between nelfinavir and M pro protein by docking (Khaerunnisa et al., 2020) . Differences observed would be consequences of the parameters considered by different authors, for example, the rotatable links in the coupling. Furthermore, this difference may also be due to the selection of different interaction conformations made by the authors. In this work, rigorous statistical criteria to make the best decision in choosing the correct cluster (detailed in item 3.1), were applied. In this sense, nelfinavir presented a broad conformational dispersion (number of clusters) and similar frequencies in the clusters with better binding energies.

All of the compounds evaluated against S glycoprotein showed very favorable binding energies, higher than those for the synthetic antiviral compounds, umifenovir (-7.47 kcal/mol), chloroquine (-7.25 kcal/mol), and hydroxychloroquine (-5.30 kcal/mol), directed to that viral target. The most favorable binding energies were obtained for glycyrrhizin (-19 .22 kcal/mol), gitoxin (-17.84 kcal/mol), dicumarol (-12.05 kcal/mol), diosgenin (-10.75 kcal/mol), and delphinidin (-8.12 kcal/mol). These results are very promising compared to those obtained in recent works, where the coupling of natural and synthetic compounds with viral glycoprotein S, were analyzed. Thus, the best binding energies obtained for S protein were: -11.55 kcal/mol to coenzyme A, -11.089 kcal/mol to flavin adenine dinucleotide, and -9.36 kcal/mol to tiludronate (Hall et al., 2020) .

In agreement with the low Gibbs free energy, numerous molecular interactions were observed for the glycyrrhizin: six hydrogen bonds and many hydrophobic interactions, all of them with residues from chains A and C from the glycoprotein. This is despite no interactions with the RBD region (recognition of ACE human receptors), chain B, were observed. Therefore, J o u r n a l P r e -p r o o f glycyrrhizin would be a good candidate to be evaluated experimentally to learn how it would affect the glycoprotein activity. Interestingly, a dietary supplement of Glycyrrhiza glabra L.

(SoriaNatural®) is already in the market demonstrating the safety of consuming this compound since it has already been approved by the corresponding organizations.

Gitoxin is a secondary cardiac glycoside from Digitalis purpurea and Digitalis lanata. This compound and others from Digitalis sp. have demonstrated action on heart muscle (Haustein et al., 1970) . The aglycone that is the active part of the compound has a cyclopentaneperhydrophenanthrene structure, which is common in steroids. This hydrophobic portion of the molecule established many hydrophobic interactions with the spike protein residues. The glycone portion of the gitoxin interacted with two N-acetyl-D-glucosamine residues from spike protein. It also formed hydrogen bonds. Many of these interactions of different nature were established in the RBD region of S glycoprotein. For this reason, we consider that it is the most promising compound with potential to inhibit spike protein.

Dicumarol is a derivative from cumarin and was successfully used as an inhibitor against influenza A . Moreover, it could be a particularly good candidate for inhibiting the S glycoprotein by blocking the RBD region, which contains human ACE2 receptors, as the compound interacted with several residues in that region. One important aspect to consider regarding dicumarol and other cumarins, is that they have already been approved in the pharmaceutical market for their effectiveness as: antidiabetic (Dwivedi et al., 2008) , antiviral (Yu et al., 2003 ), anti-inflammatory (Kontogiorgis et al., 2004 , and anticancer and antileukemia activities (Stefanachi et al., 2011; Yang et al., 2011) ; thus, they are considered safe compounds for human oral consumption. Considering that dicumarol and coumarins have powerful anticoagulant effects, their use in certain patients must be under medical control. In fact, being J o u r n a l P r e -p r o o f derived from medicinal plants with ancestral use proven for oral consumption is an important advantage of some natural origin compounds. The medicinal use of phytocomplexes allows not only increasing the desired action of the main active compounds but also decreasing adverse effects than if any of the pure active compounds were used for the same treatment.

Phytocomplexes are plant extracts or essential oils that contain more than one main active compounds and other substances that can help decreasing the adverse effects of the main compounds, contributing to the better absorption and solubility of the main compounds. This important advantage would facilitate the production and reduce the cost because fewer downstream operations would be needed as there is no need to obtain the high purity product (Cortés et al., 2004; Duraffourd et al., 1978; Hostettmann et al., 1995; Pellecuer et al., 1993; Piñeros et al., 1988; Sharapin et al., 200) .

COVID-19 is a global health problem. Despite all the efforts dedicated to stop its pandemic dissemination there is no conclusive treatment to date. This study also contains valuable information to increase the knowledge of certain structural families of compounds with great affinities by M pro protease and SARS-CoV-2 spike proteins.

Furthermore, the present study provides molecular details that set the basis to further propose structural modifications of some compounds to make the interaction between them and proteins even more effective. J o u r n a l P r e -p r o o f Table 6 . Three-dimensional (3D) molecular interaction of candidate compounds with S protein.

The S protein residues that interacted by hydrogen bonds are represented with lines and with balls and sticks those that interacted hydrophobically with the compounds. 

Ilaj bil hijamah (cupping therapy) in the Unani system of medicine: anecdotal practice to evidence-based therapy

Molecular Docking Study to Identify Potential Inhibitor of COVID-19 Main Protease Enzyme: An In-Silico Approach

Tratado de Fitofármacos y Nutracéuticos. 1ª. Reimpr. Rosario: Corpus Ediciones

Ciudad Autónoma de Buenos J o u r n a l P r e -p r o o f Journal Pre-proof

Aires: Corpus Libros Médicos y Científicos

Screening and identification of inhibitors against influenza A virus from a US drug collection of 1280 drugs

An antiviral compound from Melia azedarach L. inhibits Interferon production

Programmed ribosomal frameshifting in decoding the SARS-CoV genome

Farmacognosia, fitoquímica, plantas medicinales

Materia Médica y Terapéutica. Buenos Aires: López Etchegoyen Libreros Editores

Panorama Histórico de la Medicina Argentina. Número especial. Todo es Historia. 1ª Edición

Inhibition of H1N1 influenza A virus growth and induction of inflammatory mediators by the isoquinoline alkaloid berberine and extracts of goldenseal

Glycyrrhizin, an active component of liquorice roots, and replication of SARS-associated coronavirus

Clinical Trials, 2020. N.S. National Library of Medicine USA. COVID-19 and SARS-CoV-2

Farmacognosia: breve historia de sus orígenes y su relación con las ciencias médicas

Current lead natural products for the chemotherapy of human immunodeficiency virus (HIV) infection

Intention to participate in a COVID-19 vaccine clinical trial and to get vaccinated against COVID-19 in France during the pandemic

New insights on the antiviral effects of chloroquine against coronavirus: what to expect for COVID-19?

Synthesis and antihyperglycemic activity of novel N-acyl-2-arylethylamines and N-acyl-3-coumarylamines

Mistry BM Berberine-piperazine conjugates as potent influenza neuraminidase blocker

The species severe acute respiratory syndrome-related coronavirus: classifying 2019-nCoV and naming it SARS-CoV-2

A search for medications to treat COVID-19 via in silico molecular docking models of the SARS-CoV-2 spike glycoprotein and 3CL protease

Bioactive natural products with anti-herpes simplex virus properties

Different relationships between therapeutic and toxic actions of 16-epi-gitoxin, gitoxin and ouabain on isolated cardiac preparations

Conservation of substrate specificities among coronavirus main proteases

The novel coronavirus 2019 (2019-nCoV) uses the SARS-coronavirus receptor ACE2 and the cellular protease TMPRSS2 for entry into target cells

Religion and Gender in Medieval Culture

Software News and Updates. AutoDock4 and AutoDockTools4: Automated Docking with Selective Receptor Flexibility

Molecular Modeling of Proteins

Quantum leap of monoclonal antibody (mAb) discovery and development in the COVID-19 era

Investigación Científica en Fitoterapia: de la planta al medicamento

New insights into the structure-cytotoxicity relationship of spirostan saponins and related glycosides

Extractos Naturales de Plantas Medicinales

COVID-19 vaccine: A recent update in pipeline vaccines, their design and development strategies

Limonoids: overview of significant bioactive triterpenes distributed in plants kingdom

Molecular docking analysis of selected natural products from plants for inhibition of SARS-CoV-2 main protease

Pharmacologic Treatments for Coronavirus Disease 2019 (COVID-19): A Review

Fundamentos de tecnología de productos fitoterapéuticos

Inhibition of respiratory syncytial virus replication and virus-induced p38 kinase activity by berberine

Design, synthesis, and biological evaluation of imidazolyl derivatives of 4, 7-disubstituted coumarins as aromatase inhibitors selective over 17-α-hydroxylase/C17-20 lyase

Cryo-electron microscopy structure of a coronavirus spike glycoprotein trimer

The synthesis and antistaphylococcal activity of 9, 13-disubtituted berberine derivatives

Weekly epidemiological update -9

Cryo EM structure of the 2019-nCoV spike in the prefusion conformation

Synthesis and biological evaluation of hydroxylated 3-phenylcoumarins as antioxidants and antiproliferative agents

Response of plant secondary metabolites to environmental factors

Recent progress in the development of coumarin derivatives as potent anti-HIV agents

Crystal structure of SARS-CoV-2 main protease provides a basis for design of improved αketoamide inhibitors

Cytokine storm and immunomodulatory therapy in COVID-19: Role of chloroquine and anti-IL-6 monoclonal antibodies

Virus-encoded proteinases and proteolytic processing in the Nidovirales

Chemoenzymatic synthesis of new derivatives of glycyrrhetinic acid with antiviral activity. Molecular docking study

An overview of the safety, clinical application and antiviral research of the COVID-19 therapeutics

This research was funded by Project #2361 from Consejo de Investigaciones de la Universidad Nacional de Salta (CIUNSa, Argentina). María Antonela Zígolo is a recipient of a postdoctoral fellowship from CONICET. The authors would like to thank Dr. Jerold Last for his kind help with the language.

Declarations of interest: none.