key: cord-0807871-lrn8b78r
authors: Petit, Léna; Vernès, Léa; Cadoret, Jean-Paul
title: Docking and in silico toxicity assessment of Arthrospira compounds as potential antiviral agents against SARS-CoV-2
date: 2021-03-20
journal: J Appl Phycol
DOI: 10.1007/s10811-021-02372-9
sha: ea58c0a9de68cd514216616ea5cdc9388430a19a
doc_id: 807871
cord_uid: lrn8b78r

A race is currently being launched as a result of the international health situation. This race aims to find, by various means, weapons to counter the Covid-19 pandemic now widespread on all continents. The aquatic world and in particular that of photosynthetic organisms is regularly highlighted but paradoxically little exploited in view of the tremendous possibilities it offers. Computational tools allow not only to clear the existence and activity of many molecules but also to model their relationships with receptors identified in potential hosts. On a routine basis, our laboratory carries out a research activity on functionalities of molecules derived from algae using in silico tools. We have implemented our skills in algae biology and in modeling, as tests in order to identify molecules expressed by the genus Arthrospira showing an antiviral potential and more particularly anti-SARS-CoV-2. Using consensus docking and redocking with Autodock Vina and SwissDock, we were able to identify several promising molecules from Arthrospira: phycocyanobilin, phycoerythrobilin, phycourobilin, and folic acid. These four compounds showed reliable binding energies comprised between − 6.95 and − 7.45 kcal.mol(−1) in Autodock Vina and between − 9.285 and − 10.35 kcal.mol(−1) with SwissDock. Toxicity prediction as well as current regulations provided promising arguments for the inclusion of these compounds in further studies to assess their ability to compete with the SARS-CoV-2/ACE2 complex both in vitro and in vivo. SUPPLEMENTARY INFORMATION: The online version contains supplementary material available at 10.1007/s10811-021-02372-9.

Covid-19 is probably the most written, published acronym in this year 2020. This is for objective and obvious reasons. Covid-19 is the pandemic caused by SARS-CoV-2 (severe acute respiratory syndrome coronavirus 2) and it has overturned the entire planet as few pandemics have done, and this, in a modern era which has deployed modern and high-level medicine. Research on these pandemics takes multiple paths whether on technical, chemical, or biological levels. We then attack the biology of viruses or the defense mechanisms of their hosts. Different methods can be used for antiviral research to act against SARS-CoV-2. We can first target structural proteins in view to inhibit virus entry in human host cells or looking for functional internal viral protein inhibitors to block its replication for example. A third method would be to target human receptor for virus entry prevention.

In view to limit the impact on the human metabolism of such antivirals, we rather focused on the first critical step for the virus, its recognition by the host cell, mediated via specific structural proteins, spike proteins. As a food tech company, Algama produces innovative food ingredients made with microalgae. As part of our daily research, we developed a tool enabling to select the most relevant molecules from microalgae for different food applications according to various criteria (bioavailability, safety, drug interaction…). In the context of the Covid-19 health crisis and to contribute to the collective effort, we tested our methodology toward antiviral compound identification from microalgae. Based on a literature framework about antivirals from Arthrospira spp, the aim of the present study was to identify the most relevant SARS-CoV-2 antiviral molecules thanks to docking, and to conduct an in silico toxicity assessment of the most efficient antiviral compounds. The first part of this work focuses on the description of the virus and presents a state of the art of Arthrospira use as an antiviral agent. The second part of the article deals with the evaluation of the anti-SARS-CoV-2 activity by docking Arthrospira compounds. Finally, an in silico study of the selected molecules toxicity will be proposed.

SARS-CoV-2 belongs to the coronavirus family, composed of enveloped virus with positive strain of RNA. Its genome is contained in a capsid formed by nucleocapsid proteins, itself included in an envelope. Three structural proteins are characteristic of coronavirus: membrane protein, envelope protein, and spike protein (S), a glycoprotein responsible for virus host cell attachment (Li 2016) . Coronavirus took their name from "corona" (crown in Latin), formed by S protein protuberances on their surface.

The spike (S) glycoprotein ectodomain consists in a trimeric structure. Each monomer is composed of two subunits: S1, involved in the host cell receptor recognition; and S2, responsible for the membrane-fusion mechanism. S1 subunit contains two major domains, N-terminal domain (S1-NTD) and C-terminal domain (S1-CTD) (Li 2016) (Fig. 1) . The latter constitutes the receptor binding domain (RBD) and is composed of two subdomains: a core structure and a receptorbinding motif (RBM) (Li et al. 2005) . Zhou et al. (2020) recently determined the ability of SARS-CoV-2 to use human angiotensin-converting enzyme 2 (ACE2) as receptor to engage virus attachment to the host cell via its RBD, as for SARS-CoV. However, they excluded binding with other coronavirus receptors such as APN and DPP4 (Li 2016) . ACE2 identification as main receptor of SARS-CoV-2 in human has been confirmed shortly after by two other research teams (Hoffmann et al. 2020; Ou et al. 2020) .

ACE2 is a metallocarboxypeptidase involved in arterial pressure system via cleavage of angiotensin peptides (Donoghue et al. 2000) . It has been shown to be expressed in lung alveolar epithelial cells, in enterocytes of small intestine and in vascular endothelium (Hamming et al. 2004) , explaining the location of Covid-19 pathologic symptoms. ACE2 binding site to SARS-CoV-2 was found to be within its extracellular peptidase domain .

Once coronavirus RBD is linked to ACE2 peptidase domain, the prefusion S protein becomes unstable and is submitted to a first proteolytic activation via a host cell protease, between S1 and S2 subunits (Walls et al. 2017 ). In addition to S1/S2 cleavage, a second site, critical for viral fusion, has been determined within the S2 domain, on the S2′ site ). This second proteolysis results in the shedding of S1 subunit and exposes an internal fusion peptide. This one is essential for fusion of the virus and host membrane via S2 subunit, leading to the viral genome release in ACE2. This peptide is highly conserved among coronavirus (Madu et al. 2009 ) and composed of serine residues 798 to phenylalanine 815 in SARS-CoV.

TMPRSS2 was identified in lung cells as the main human protease used by SARS-CoV-2 for S protein modification (Hoffmann et al. 2020) . TMPRSS2 belongs to the type II transmembrane serine proteases (TTSPs) which have been involved in the spread of various respiratory viruses (Choi et al. 2009 ). Expressed in airway epithelial cells (Böttcher et al. 2006) , Shulla et al. (2011) further pointed the colocalization of TMPRSS2 and ACE2 during an immunoprecipitation assay. In this context, SARS-CoV-2 using the same receptor (ACE2) and host cell protease (TMPRSS2) as SARS-CoV, this colocalization appears to sorely facilitate its entry and spreading in human host cells. Gui et al. (2017) determined different conformational prefusion states of SARS-CoV S protein. Each Spike RBD on the S1 subunit of SARS-CoV can be found in two conformational states, one "up" and one "down." In down position, SARS-CoV RBD was reported as inaccessible for ACE2 recognition due to steric clashes. On the contrary, when a trimer is in "up" position, RBD is then exposed, allowing its binding with host cell ACE2 (Gui et al. 2017; Yuan et al. 2017) .

At this day, two conformational states of the S protein have been found in SARS-CoV-2 (Wrapp et al. 2020; Walls et al. 2020) . As for SARS-CoV, one presents all RBDs in down state (not completely closed according the Cryo-EM assays) and the second shows one of the three RBDs exposed in an open state. Like for its predecessor, SARS-CoV-2 RBD binding with ACE2 seems to be impossible when the Spike trimer is entirely in down conformation. On the contrary, with one RBD in up position, no steric clash with S protein was S1 90°a S2 S1-NTD S1-CTD b c Fig. 1 a Cryo-electron microscopy structure of prefusion trimeric SARS-State of the art on antivirals derived from algae, in particular Arthrospira (Tables 1 and 2) The term "algae" actually includes a wide variety of organisms with a complex evolutionary history. These organisms have features in common such as the ability to do photosynthesis and being predominantly aquatic. We can distinguish two main categories of algae: macroalgae, multicellular eukaryotes; and microalgae, in the largest sense of the word, including unicellular eukaryotes and prokaryotes. Most macroalgae are found in seawater and are called "algae," but some can also be found in freshwater (e.g., Oedogonium, Cladophora and Spirogyra). Microalgae are found in both environments, but the most consumed species (Arthrospira and Chlorella) generally grow in freshwater. These aquatic plants are known to be rich in proteins and lipids, but also offer a wide range of immunostimulatory, antioxidant, and antiviral properties. Numerous publications have reviewed the generic properties of algae in these areas (Mimouni et al. 2012) . For example, the immune response in the cytokines which is tested by Talukdar et al. (2020a Talukdar et al. ( , 2020b with astaxanthin. Griffithsin, a lectin isolated the from red alga Griffithsia, has shown activity on protein S of MERS-CoV (Middle East Respiratory Syndrome Coronavirus) (Millet et al. 2016) . As for microalgae, the very numerous studies have been compiled by De la Jara et al. (2018) and Perumal and Sundararaj (2020) . In addition, few studies have been identified in the literature as dealing with the antiviral potential of microalgae compounds against SARS-CoV-2. For instance, astaxanthin as well as linolenic acid (C18:00) have been mentioned as potential candidates for the treatment or the prevention of Covid-19 (Talukdar et al. 2020b ). However, these are pre-published studies and the confirmation of the robustness will be given upon the final peer-reviewed publications.

The genus Arthrospira is a genus of prokaryotic filamentous cyanobacteria. The best-known species, Arthrospira platensis and Arthrospira maxima, are mostly used as food supplement due to their nutritional properties: high protein content, vitamins, minerals, etc. Moreover, Arthrospira also contains various metabolites with beneficial health properties such as polyphenols, carotenoids, and sterols (de la Jara et al. 2018 ). Due to drug-resistant virus strains emergence, and because of adverse effects caused by usual synthetic antiviral treatments, increasingly researchers are looking into more natural approaches. This explains a growing interest in antivirals from natural sources, in algae in general and in Arthrospira for example in particular. Hence, in vitro and in vivo studies have been conducted to evaluate Arthrospira compounds as antiviral against various viruses such as influenza, herpes simplex virus (HSV), and hepatitis C (HCV) as presented in Tables 1  and 2. Data from scientific literature were critically reviewed in order to highlight results of the most rigorous studies: presence of controls, replicates, measurement of toxicity, and inhibition efficacy. The following information were collected and analyzed: CC 50 , cytotoxic concentration of the compounds decreasing in vitro cell viability to 50%; EC 50 , effective concentration needed to inhibit 50% of virus; and therapeutic index (TI or SI), referring to the ratio CC 50 /EC 50 . A high therapeutic index (TI), or selectivity index (SI), is preferable for a drug to have a favorable safety and efficacy profile. Radonic et al. (2010) evaluate the effect of sulfurcontaining exopolysaccharides from A. platensis (TKV3). According to their study, anionic polysaccharides from Arthrospira showed antiviral activity against enveloped virus such as VACV (EC 50 (TK V3) = 0.78 μg.mL −1 ) via interactions with viral membrane glycoproteins and thus inhibition of virus binding to host cells. In the work of Chen et al. (2016) , Arthrospira extract was capable of inhibiting in vitro influenza viral replication and plaque formation (Table 1) , by targeting hemagglutinin, an influenza virus surface glycoprotein. Abdo et al. (2012) , El-Baz et al. (2013) , Deyab et al. (2020) , and Hetta et al. (2014) focused on A. platensis antiviral effect on adenovirus and Coxsackievirus (CV). Most of them evaluate Arthrospira methanol extract activity on Hep-2 cell in vitro model. In all cases, Arthrospira appears as a good antiviral agent: reducing virus titer from 50% and preventing virus from attaching to host cell receptor when used as pre-treatment with TI(CVB3) = 30 and TI(RV) = 45 (Deyab et al. 2020) . Efficiency of methanolic extract is attributed to the presence of polar compounds binding to viral capsid. Moreover, El-Baz et al. (2013) found that Arthrospira (Jang and Park 2016) ethanolic extract was active against non-enveloped RNA and DNA viruses. However, authors did not compare Arthrospira effect with a control drug since there is no available drug against enteric virus according to them. Elsewhere, Rechter et al. (2006) compared spirulan-like substance to ganciclovir (GCV), a reference drug used as anti-herpesvirus. According to their in vitro experiments, preincubation of substances enhances antiviral activity via virus entry and replication inhibition: EC 50 (TK-V3a) = 1.4 ± 0.3 μg.mL −1 , close to drug reference EC 50 = 0.7 ± 0.1 μg.mL −1 in post-incubation (Table 1) .

Due to the emergence of drug-resistant Herpes simplex strains, several studies have focused on the search for alternatives to synthetical acyclovir in the prevention and treatment of Herpes simplex virus (HSV) 1 and 2. For instance, hot water extract of A. maxima provided EC 50 of 0.069 mg.mL −1 against HSV-2. The authors then attributed this antiviral activity to high polar compounds (Hernández-Corona et al. 2002) . In a second study, phosphate buffer extract of A. fusiformis showed a viral infection inhibition almost as high as acyclovir (85% against 99%) during the intracellular replication period (Sharaf et al. 2013 ). In the same way, Chirasuwan et al. (2009) evaluated in vitro A. platensis extract as potential HSV antiviral agent to propose an alternative to the synthetic anti-herpes drug acyclovir. Sulfoquinovosyl diacylglycerol (SQDG) was responsible for antiviral activity of A. platensis lipid fraction (EC 50 = 6.8 μg.mL −1 ) with an effect comparable to the reference drug (EC 50 = 1.5 μg.mL −1 ), without toxicity for cells. This activity is due to DNA polymerase inhibition thanks to SQDG interaction with different regions of the enzyme. Furthermore, Lee et al. (2007) highlighted antiviral potency of spirulan molecules from A. platensis.

Best in vitro results were obtained when compounds were added to the medium during HSV-1 and 2 infections. Na-Sp then showed the highest selectivity indexes of 13,000 and 17,000, with very low ECs 50 (0.63 and 0.41 μg.mL −1 ). Otherwise, Shalaby and Shanab (2010) attributed the observed hepatitis A (RNA virus) and HSV-1 antiviral effect of aqueous and phosphate buffer extracts to the sulfated polysaccharides and tannins. According to them, these compounds may interfere at different viral stages such as attachment and penetration of the virus.

In the case of Kok et al.'s (2011) study, different extracts from microalgae were compared to conventional drugs (acyclovir and foscarnet) for their antiviral capacity against Epstein-Barr virus (EBV). Methanol extract from A. platensis reduced the cell-free EBV DNA load in B95-8 cells with an EC 50 of 0.021 μg.mL −1 with an CC 50 of 166 μg.mL −1 (TI = 7905) (Kok et al. 2011) .

Peptides were also evaluated as antiviral agent. Indeed, Jang and Park (2016) found that a peptide isolated from A. maxima displayed in vitro antiviral activity since it inhibited reverse transcriptase activity in HIV-1-infected cell by 90% compared to the no-peptide assay. Moreover, 0.75 mg.mL −1 of this peptide inhibited HIV-1 p24 antibody production by more than 95%, without toxicity for the cells. Similarly, Ayehunie et al. (1998) evaluate HIV antiviral activity of Arthrospira. According to their results, both polysaccharide fraction and the fraction depleted in polysaccharide and tannin of A. platensis inhibited in vitro HIV-1 replication. Preincubation of polysaccharide fraction extract with cells showed better results due to polysaccharides binding to CD4 receptors, thus preventing virus attachment to the host cell through its gp120 envelope glycoprotein.

According to in vitro studies presented above (Table 1) , Arthrospira has antiviral properties against different types of viruses: Herpes simplex, hepatitis, vaccinia, coxsakievirus…. It has no, or very little, cytotoxicity to cells, and low doses (EC 50 ) are sufficient to induce an effect (1-10 μg.mL −1 ) with high therapeutic indices. It seems that Arthrospira compounds act during the early stages of viral infection. According to the main hypotheses evoked, compounds of polar nature and mainly sulfated polysaccharides and peptides are responsible for interactions with viral envelope glycoproteins, thus preventing binding to the host cell. Calcium spirulan, a sulfated polysaccharide from Arthrospira, is the most cited compound with antiviral activity against various enveloped viruses. The sulfate groups as well as the carboxyl groups of polysaccharides, such as Ca-Sp isolated from the extract, have negative charges which can react with the basic amino acids of viral proteins and block the interaction with cellular receptors.

In order to confirm in vitro results, a second state of the art was carried out focusing on human clinical studies involving the use of Arthrospira as an antiviral agent ( Table 2) .

Most of the in vivo studies focused on Arthrospira supplementation effect on HIV-infected patients. Today, there is no treatment capable of eradicating AIDS. HIV-infected patients still carry the virus and they can spread it to others. So far, there are five major classes of anti-HIV drugs that target distinct steps in the virus life cycle: reverse transcriptase inhibitors, protease inhibitors, fusion inhibitors, integrase inhibitors, and multidrug combinations (Jang and Park 2016) . To improve therapeutic potential of existing medicines, research on Arthrospira-based novel treatment have been conducted.

In this context, Winter et al. (2014) conducted an in vivo study on 73 HIV-infected women to evaluate Arthrospira supplementation effect. The study was a 3-month pilot, randomized, double-blind, and placebo-controlled intervention. No significant clinical effect of Arthrospira supplementation was observed regarding viral load and CD4 cells; nevertheless, Arthrospira supplementation showed a positive effect on weight stabilization and protection against opportunistic infections. Similarly, Teas and Irhimeh (2012) found that Arthrospira strengthen the immune system of affected patients. Unfortunately, sample size of selected patients for the present study was too small (11) to conclude that there were significant effects of Arthrospira supplementation. Ngo-Matip et al. (2015) conducted a single-blind, randomized, multicenter study on 320 HIV-1 ARV-naïve participants for 12 months to analyze Arthrospira supplementation effect compared to local diet, both with standard HIV therapy. This work confirmed previous results, showing that Arthrospira supplementation intake positively and significantly stimulated the immune system and inhibited virus replication of HIV-infected subjects. As in the previous cases, Arthrospira contributed to the proper functioning of immune system and thus helped to limit the appearance of opportunistic diseases. Similar results were also obtained by Azabji-Kenfack et al. (2011) , indicating that Arthrospira supplementation helps HIV-infected and malnourished adults to improve their immune defenses (Azabji-Kenfack et al. 2011) . For this experiment, authors compared effects of Arthrospira and soya bean supplementation. Like other studies, Arthrospira supplement improved general state of health by promoting weight gain of malnourished adults infected with HIV. Regarding the baseline and the results obtained with soya beans, Arthrospira consumption significantly reduced viral load and increased CD4 count (Azabji-Kenfack et al. 2011) .

Two in vivo studies listed in Table 2 focused on Arthrospira effect in the case of hepatitis C infection (HCV). In the work of Gomaa et al. (2017) , Arthrospira was added to usual thalassemic drug for 6 months to evaluate the effect on thalassemic HCV-infected children. As patients are suffering from two different pathologies, it was difficult to affirm that Arthrospira was responsible for health improvement. However, the authors conclude that Arthrospira stimulated the immune system of thalassemic children infected with HCV since CD4 and CD8 number was increased after 6 months of Arthrospira intake. They also state that phycocyanin from Arthrospira stimulated hematopoiesis, by inducing erythropoietin hormone (EPO) release, thus helping thalassemia children. Moreover, Yakoot and Salem (2012) conducted an in vivo study comparing the effect of Arthrospira and silymarin against hepatitis C. For the present study, 66 patients with chronic hepatitis C virus infection had been randomized, double-blinded, and treated with Arthrospira or silymarin for a period of 6 months of treatment. According to their results, Arthrospira helped to improve the general well-being of patients and has led to a loss or reduction of detectable hepatitis C virus RNA for six patients. Compared to the second group treated with silymarin, these results were no significant (0.12) ( Table 2) . Nevertheless, the present study did not include a placebo control since silymarin possess itself some benefits in the treatment of viral hepatitis. Thus, it does not allow to conclude on these results.

According to these in vivo studies, it appears that Arthrospira improves the general health of AIDS and HCVinfected patients by contributing to restore their body's defense mechanism against infectious immune system disorder, thus limiting opportunistic diseases. According to these results, Arthrospira may significantly affect virus progression according to some studies. However, most studies conducted in vivo on human patients have evaluated the effect of daily Arthrospira supplementation in addition to the usual antiviral treatment. These studies were mostly randomized and doubleblind, but the study duration was often short with a small number of patients and did not always present a placebo control. It could be interesting to evaluate the effect of Arthrospira alone, without other antiviral treatment, in a rigorous long-term study to confirm its antiviral potential.

According to this literature review, Arthrospira acts as antiviral against various type of viruses. Indeed, significant antiviral effects were observed in vitro, and in vivo studies (Table 1 and 2) confirmed that Arthrospira is well tolerated by patients and even improved their general state of health by increasing immune response. Arthrospira can thus be administered in addition to current gold standard therapy in order to improve its effects, or it can also be considered as an interesting alternative to certain drugs with serious adverse effects. Many compounds such as sulfated polysaccharides (calcium-spirulan), fatty acids (Hetta et al. 2014) , and proteins (C-phycocyanin, cyanovirin-N, microvirin) have been identified as responsible for this antiviral activity. Indeed, these compounds can bind to glycoprotein virus envelope, and thus hinder virus attachment to its host cell. Nevertheless, according to Ayehunie et al. (1998) , it appears that compounds other than polysaccharides and tannins, present in the aqueous extract of Arthrospira, may also have an antiviral activity. To go further, it might therefore be interesting to evaluate the antiviral potential of other Arthrospira metabolites.

A pre-published work publicly available highlighted Arthrospira potential on two viral proteins of SARS-CoV. However, it was based on SARS-CoV-1 database, and it is p e n d i n g f o r p u b l i c a t i o n ( D O I : 1 0 . 2 6 4 3 4 / chemrxiv.12051927). Very few authors have taken an interest in the specific case of Arthrospira against SARS-CoV-2, and most of the articles present a literature survey without in silico, in vitro, or in vivo evaluation. The use of docking in order to test Arthrospira as a candidate against SARS-CoV-2 is thus original.

In order to propose new antiviral leads from Arthrospira, the present study evaluates in silico the antiviral potential of Arthrospira metabolites using an updated database specific to SARS-CoV-2 thanks to docking tools. Molecular docking is a computational method which simulates different orientations of a receptor and a ligand to predict the most stable complex, with a minimal free energy. Widely used to new drug discovery, it is composed of two main steps. First, programs perform a conformational sampling of ligands in the receptor active site. Then, an algorithm (scoring function) scores each position. The latter predicts binding energy of compounds and ranks all the docking positions by increasing energy. The lower the binding energy, the greater the affinity between the ligand and its receptor.

Autodock Vina, on one hand, a popular program for molecular docking, involves a hybrid scoring function based on both empirical scoring functions and knowledge-based potentials. Indeed, it is largely inspired of X-score which considers van der Waals (VDW) interactions, hydrogen bonds, hydrophobic effects, and effective rotatable bonds to determine a binding score (Huang et al. 2010; Trott and Olson 2010) . On the other hand, SwissDock is based on the EADock DSS (Dihedral Space Sampling), a forcefield scoring function. Indeed, it involves CHARMM22 protein forcefield using VDW and electrostatic grids to determine the interaction energy between the protein and ligands (Grosdidier et al. 2011a) .

In order to contribute to the collective effort to find a solution to the current global Covid-19 crisis and to further explore the antiviral potential of Arthrospira, this article investigates antiviral activity of Arthrospira molecules from Algama internal database. This database is composed of 48 high-value molecules from Arthrospira genus. As control molecules, we choose an antiviral (remdesivir) and two molecules with high docking score according to the literature (linoleic acid and pavetannin C1) (Prasanth et al. 2020; Toelzer et al. 2020 ).

Two crystal structures of the receptor-binding domain (RBD) of SARS-CoV-2 Spike protein bound to the cell receptor ACE2 have been determined (Lan et al. 2020; Wang et al. 2020 ) and deposited on the Protein Data Bank under the respective codes 6M0J and 6LZG. All the docking assays were executed with the firstly deposited crystal structure of the complex, corresponding to the 6LZG identifier . Complete RBD part (with all residues) was extracted from this crystal complex using Discovery Studio Visualizer, thus expecting having the most relevant conformation of the RBD in the Spike/ACE2 complex (Forli et al. 2016) . Receptor structure was then prepared, heteroatoms (water included) were deleted, and hydrogens were added. Protonation states were assigned thanks to H++ 3.2v at physiological pH (7.4) with AMBER ff 14SB parameters (Gordon et al. 2005; Myers et al. 2006; Anandakrishnan et al. 2012 ). Finally, protein structure was energy minimized with Chiron tool (Ramachandran et al. 2011 ) and the file was saved as PDB file.

Ligands were divided into three categories: small molecules (mainly secondary metabolites), Arthrospira peptides, and the antiviral drug remdesivir. Small molecule and drug structures were obtained from PubChem and downloaded as SDF format. Peptide sequences were taken from literature and converted in SMILES format via NIH Cactus translator, and finally in SDF format in Open Babel (O'Boyle et al. 2011 ). Ligands were then prepared for docking with both Autodock Vina and SwissDock programs according to the following protocol. Each molecule was redocked with both software in view to discriminate false-positive results.

All the 52 ligands from the Algama database were energy minimized before being finally converted in PDBQT format on PyRx (Dallakyan and Olson 2015) . Docking was executed with Autodock Vina (Trott and Olson 2010) implemented in the PyRx platform (Dallakyan and Olson 2015) . The grid box was manually adjusted around active residues identified above, with the coordinates X: − 37. Exhaustiveness was set on 8, equivalent to the short mode, since no statistically significant difference in accuracy was detected using the short, medium, or long mode of Autodock Vina (Nguyen et al. 2020) .

Binding affinity result files were exported as CSV for analysis and docked structures were visualized in Discovery Studio Visualizer in view to obtain a complete 3D structure of each complex. Ligand interactions were determined with this same tool and analyzed to target better candidates for SARS-CoV-2 RBD binding. In literature, amino acids of Spike RBD are classically numbered from 333 to 527, as mentioned in Table 3 . Nevertheless, during docking preparation of the protein, amino acids numbers were modified by the H++ program, thus going from 1 to 195. Interaction prediction files were also generated with these new numbers (Figs. 2, 3, 4, 5, 6, 7, 8, 9, and 10) . However, for a better understanding and being aligned with other papers, residue numbers were post-docking reassigned from 333 to 527 in the results and discussion parts, thanks to a translation table (Table 9 Supplemental Data).

SDF files were converted into .mol2 format with 3D coordinates and protonation at physiological pH with Open Babel software (O'Boyle et al. 2011 ). All ligands were submitted to a control with SwissParam (Zoete et al. 2011 ) before docking on SwissDock (Grosdidier et al. 2011a (Grosdidier et al. , 2011b . Previous PDB file of Spike RBD was used as receptor file again. Docking was set in accurate mode with the same grid box coordinates and dimensions used for the Autodock Vina assays. Pose with the best binding energy was selected for each ligand and respective result file was exported as CSV format.

In silico toxicity assessment Toxicity of selected molecules was assessed via the AdmetSAR 2 prediction tool (Yang et al. 2019 ), using canonical smiles from PubChem for compound identification. SwissADME and Osiris platforms were also used to complete toxicity data (Daina et al. 2017) .

Acute oral toxicity was determined through toxicity class. Mutagenicity, carcinogenicity, hepatotoxicity, and reprotoxicity were also predicted. In order to avoid interaction with other potential treatment, study of the main cytochrome as well as P-glycoprotein interactions was needed. Finally, druglikeness parameters were assessed such as gastrointestinal absorption, blood-brain barrier permeability, or Lipinski's rule.

After docking twice the 52 molecules from Algama's internal database with Spike RBD in both Autodock Vina and SwissDock, the first mode was selected for each compound, corresponding to the best binding affinity (Table 8 Supplemental data). Both their RMSD/ub and RMSD/lb were lower than 2 Å (Trott and Olson 2010) .

In order to identify best compounds and false-positive results, averages, SDs as well as coefficients of variation were calculated for the two assays of each program (Table 8 Supplemental data). Thus, for one molecule, binding energies from Autodock Vina trials were averaged, like those of SwissDock. Two lists of averaged binding energies were then obtained and ranked from the smallest to the largest. In general, redocking confirmed the results since coefficients of variation were almost all behind 5% with both Autodock Vina and SwissDock. Redocking can be therefore considered as a reliable method regarding the two conducted assays for each program.

The first 15 compounds were selected on each list and then compared to each other. Autodock Vina provided the following 15 compounds: dieckol, phycoerythrobilin, pavetannin C1, phycocanobilin, rutin, catechin, β-carotene, phycourobilin, rosmarinic acid, epicatechin, astaxanthin, apigenin, folic acid, remdesivir, and poriferasterol (Table 4) . Averaged binding energies were comprised between − 6.7 and − 8.1 kcal.mol −1 . Regarding substances with lowest binding energies obtained with SwissDock, the first 15-molecule list w a s s l i g h t l y d i f f e r e n t , c o n s i s t i n g o f p e p t i d e GGTCVIRGCVPKKLM, folic acid, phycoerythrobilin, phycocyanobilin, phycourobilin, nervonic acid, pavetannin C1, (Table 5) . Binding scores went from − 7.9 to − 11.5 kcal.mol −1 . Consensus docking method aimed to identify false-positive results. To do this, common molecules of both tops were selected and grouped in a unique list. Table 6 therefore presents the nine best molecules regarding their binding energy selected for the rest of the study, according to docking results. Among them can be found eight Arthrospira molecules: phycoerythrobilin, phycocyanobilin, dieckol, phycourobilin, folic acid, rutin, β-carotene, and astaxanthin. Pavetannin C1, which was a positive control found to bind with Spike RBD in a previous paper (Prasanth et al. 2020) , also showed competitive values with low coefficients of variation (< 5%) in the present study, thus supporting experimental results. The 12 remaining compounds were considered as probably false-positive values. Nine common compounds stand out among the top 15 molecules with the best binding affinity to the RBD domain for Autodock Vina and SwissDock. The remaining 12 compounds from the two software programs are therefore considered false positives. These are catechin, rosmarinic acid, apigenin, remdesivir, and poriferasterol, which are ranked high with Autodock Vina, but have high binding energies in the case of the SwissDock models. Similarly, pGGTCVIRGCVPKKLM, nervonic acid, eicosadienoic acid and pLDAVNR, arachidonic acid, and dihomo-gamma-linolenic acid present in the top 15 molecules with Swissdock are considered poor candidates in the case of Autodock Vina because of their low binding energies (higher than − 5 kcal.mol −1 ). Analysis of specific interactions between Spike RBD and selected substances was then performed. Binding residues are recorded in Table 6 .

Phycoerythrobilin binds to Spike RBD via six van der Waals interactions, with the following residues: LYS417, TYR453, GLN493, TYR495, PHE497, and GLN498. In addition, 2 alkyl and π-alkyl interactions engage TYR449 and LEU455, and a π-cation interaction binds phycoerythrobilin with ARG403 of Spike RBD. TYR505 is involved in two different kinds of bonds: one π-π T-shaped bond and one H-bond. Two other Hbonds link respectively GLY496 and ASN501 of Spike RBD. The last two H-bonds are located on the SER494 residue (Fig. 2) . Thus, with three key binding amino acids involved in different interactions including two H-bonds on them, phycoerythrobilin shows interesting binding affinities from − 7.4 to − 9.63 kcal.mol −1 with coefficients of variation < 1% for AutodockVina and Swissdock, respectively. These results provide us very good arguments for deeper research on a phycoerythrobilin-based antiviral.

Phycocyanobilin/Spike RBD complex is composed of eight van der Waals interactions, with residues ARG403, TYR449, TYR453, GLN493, TYR495, PHE497, ASN501, and TYR505 (Fig. 3) . Four H-bonds are involved between phycocyanobilin and Spike RBD, Non-bold residues reflect unfavorable interactions between the molecule studied and the spike RBD, those in bold are directly involved in the link between ACE2 and spike RBD. Stars were added to amino acids involved in one (*) or two (**) H-bond ( 

Dieckol has seven van der Waals interactions spread on residues ARG403, TYR449, PHE456, GLN493, TYR495, PHE497, and GLN498. Two π-sigma, one π-alkyl as well as one π-π T-shaped bonds engage respectively LEU455, LYS417, and TYR505 (Fig. 4) .

The four present H-bonds involve key binding residues of Spike RBD namely TYR453, GLY496, THR500, and ASN501. Despite an unfavorable acceptor-acceptor interaction on GLU406, dieckol has six different bonds with Spike RBD and possesses the lowest binding energy in Autodock Vina program with 8.1 kcal.mol −1 . SwissDock values are also low with −8.375 kcal.mol −1 , thus positioning this molecule as one of the best candidates resulting from the present study.

Six van der Waals as well as five alkyl, π-alkyl, πcation, and π-π T-shaped interactions link Spike RBD to phycourobilin, involving residues ARG403, GLU406, LYS417, TYR449, TYR453, LEU455, GLN493, TYR495, PHE497, GLN498, and TYR505. In addition, five H-bonds interact with the compound and SER494, GLY496, and ASN501 amino acids (Fig. 5) . Three of them are directly in the binding pocket of Spike RBD, leading to great binding energies with − 7.1 kcal.mol −1 for Autodock Vina and − 9.3 kcal.mol − 1 with SwissDock.

Three kinds of interactions were predicted between folic acid and Spike RBD. It accounts the greater number of H-bonds with a total of seven, engaged respectively with ARG403, LEU492, GLY496, ASN501, and TYR505. Among them, three link directly with key binding residues. Moreover, eight VDW and one carbon hydrogen bonds were also detected, in connection with residues TYR449, LEU452, TYR453, GLU484, PHE490, GLN493, SER494, TYR495, PHE497, and GLN498 (Fig. 6) . Binding energies went from − 6.9 to − 10.35 kcal.mol −1 with both coefficients of variation lesser than 1%, making folic acid a reliable candidate for Covid-19 solutions.

Pavetannin C1 presents seven VDW bonds as well as one carbon-hydrogen bond and one π-π T-shaped interaction with respectively ARG403, TYR449, PHE490, TYR495, PHE497, GLN498, ASN501, GLN493, and TYR505. The latter residue can also bind with pavetannin C1 thanks to a H-bond, with three others on LEU492, SER494, and GLY496 (Fig. 7) . One of them is engaged with key binding residues, namely GLN496, and GLY493 is connected thanks to a carbon hydrogen bond. Binding scores go from − 7.3 to − 8.52 a b Fig. 3 Phycocyanobilin docked with spike RBD of SARS-CoV-2. a 3D structure of the complex, phycocyanobilin (ligand) is represented in sticks, Spike RBD (receptor) is in blue. b Interaction map of the ligand/receptor complex kcal.mol −1 with coefficients of variation larger than for other compounds but still lesser than 5%.

Rutin binds to Spike RBD via seven van der Waals interactions, with the following residues: TYR449, TYR453, SER494, TYR495, PHE497, GLN498, and GLY502. In addition, six H-bonds link respectively ARG403, GLN493, GLY496, ASN501, and TYR505 of Spike RBD (Fig. 8) .

Those located on the GLN493, GLY496, and ASN501 residues belong to the binding pocket of Spike RBD (Table 3) . Thus, with three key binding amino acids involved in four different H-bonds on them, rutin shows interesting binding affinities from − 7.25 to − 8.08 kcal.mol −1 with coefficients of variation < 1%. These results provide us good arguments for deeper research on this compound for the elaboration of a new antiviral. 

In the same way as β-carotene, the algorithm does not predict any H-bond for astaxanthin. Astaxanthin/Spike RBD complex is then composed of eight VDW interactions, with residues TYR449, LEU452, LEU455, GLU484, ASN487, LEU492, GLN493, and SER494. Two π-sigma interactions are also involved between the protein and the ligand, on TYR489 and PHE490 residues.

Finally, four alkyl and π-alkyl interactions are located on PHE456, ALA475, TYR489, and PHE490. Only ALA475 and TYR489 share bonds with astaxanthin among key binding residues. This compound therefore appears to have less strong bonds with Spike RBD, but it still has good binding energies. Indeed, as mentioned previously, binding energies calculated for astaxanthin are similar to those for β-carotene and rutin with values included between −7.05 and − 8.19 kcal.mol −1 . 

After docking assays and binding interaction analysis, an in silico toxicity study was conducted on the top 10 molecules with AdmetSAR 2, SwissADME, and Osiris tools. Substances were analyzed regarding specific toxicity criteria including both general toxicity criteria (mutagenicity, acute oral toxicity, carcinogenicity, hepatotoxicity, and reprotoxicity), druglikeness parameters (GI absorption, BBB permeability, or Lipinski's rule), and ability of each compound to inhibit human cytochromes, thus identifying possible risk of drug interactions. Each toxicity tool provided different types of results, which are compiled in Table 7 . AdmetSAR 2 focused on toxicity risks, cytochrome a b inhibitions, and blood-brain barrier permeability. A positive or negative result was provided with a p value for each compound. When a negative result was obtained, corresponding to a low-risk probability, "low" was recorded in the table. In the same way, when a criterion responded positively but with a p value < 0.7, it was considered as not enough statistical significance so the risk was determined as "mid," and finally, in case of positive response combined with a p value > 0.7, a high risk is registered in Table 7 . For SwissADME, no p value was provided; results were treated on a binary way. Therefore, a positive result was associated with a "high" risk whereas a negative result gave a "low" risk. Finally, Osiris results consisted in a range of colors: red orange or green. These colors were gradually assigned to risk levels: red color, "high" risk/probability; orange color, "mid" risk/probability; and green color, "low" risk/probability. Some criteria were covered by two toxicity tools. In that case, the higher risk was kept and recorded in Table 7 . Nevertheless, detailed data are available in Table 10 (Supplemental data). All of the compounds presented acute oral toxicities included in the third class of toxicity except pavetannin C1 which belongs in the fourth class. The third class corresponds to LD 50 that comprised between 500 and 5000 mg.kg −1 , whereas the fourth class included LD 50 greater than 5000 mg.kg −1 . These ones are respectively equivalent to the fourth and fifth classes of toxicity in the Hodge and Sterner scale. Selected molecules can be therefore considered as safe for this criterion. Regarding other toxicity risks, none of the tested compounds presented carcinogenicity, but rutin responds positively to the mutagenicity and hepatotoxicity criteria with high probability levels. Hepatotoxicity risk is also detected for folic acid while astaxanthin and pavetannin C1 were predicted to be potential reprotoxic compounds. It is thus necessary to remain vigilant regarding these parameters.

Except folic acid and rutin, most of the compounds may cause P-glycoprotein inhibition, a major actor in drug metabolism (Rautio et al. 2006) . Such inhibition may enhance bioavailability of drugs by accumulation in cells which, if it is not the wanted effect, could lead to toxic effects. Thus, particular attention should be paid to drug interactions with the seven concerned molecules, and more particularly with dieckol and pavetannin C1. Indeed, both substances show significant risk of interaction with human cytochromes CYP1A2, CYP2C19, and CYP2C9 for dieckol and CYP2C9 for pavetannin C1. Also involved in drug metabolism, interactions could modulate activation of other drugs, their absorption, or their elimination. Nevertheless, all the tested compounds showed no relevant probability of inhibition of cytochromes CYP2D6 and CYP3A4. The latter is the most represented of CYP450s in human, distributed in both hepatocytes and enterocytes. Involved in the metabolism of almost half of drugs, it is therefore essential (Martínez-Jiménez et al. 2007) . Moreover, none of folic acid and rutin present significant risk of cytochrome inhibition which make these compounds safe for a potential administration in patients already consuming other treatments.

In another way, all the candidates were predicted to be absorbed by the gastrointestinal tract but at a low level. Six molecules showed favorable parameters regarding druglikeness, namely astaxanthin, pavetannin C1, phycocyanobilin, phycoerythrobilin, phycourobilin, and rutin. Among them, three compounds also responded to Lipinski's rule: phycocyanobilin, phycoerythrobilin, and phycourobilin. Combined with their docking results, these candidates seem to be among the most relevant for the fight against SARS-CoV-2. Nevertheless, in vitro and in vivo studies remain essential to complete these results.

In the global context of the Covid-19 pandemic, joint advances in research are essential. Spike protein, a major structural protein of the SARS-CoV-2, constitutes a privileged target for antiviral research, as the first interactions between virus and human host cell take place via the Spike protein RBD, interacting with human ACE2 receptor.

Taking advantage of our skills with in silico studies, we choose to apply this knowledge as a test in the Covid-19 frame. Using consensus docking and redocking with largely widespread software (Autodock Vina and SwissDock), 8 Arthrospira molecules among 48 were identified as potential candidates with antiviral activity against SARS-CoV-2: phycoerythrobilin, phycocyanobilin, dieckol, phycourobilin, folic acid, rutin, β-carotene, and astaxanthin. Thanks to these methods, false-positive results have been discarded to keep the most reliable results.

Phycocyanobilin, phycourobilin, and phycoerythrobilin are phycobilins. The tetrapyrrolic structure of these chromophores is involved in light harvesting of cyanobacteria such as Arthrospira (Akimoto et al. 2012; Rodrigues et al. 2018) . Great binding energies were obtained for these compounds from Autodock Vina and SwissDock with respectively − 7.45 and − 10.35 kcal.mol −1 , − 7.25 and − 9.355 kcal.mol −1 , and − 7.1 and − 9.285 kcal.mol −1 for phycoerythrobilin, phycocyanobilin, and phycourobilin. Low toxicity risks were also predicted with in silico methods for these three compounds. Moreover, phycocyanin and phycoerythrin are authorized by both FDA and European Commission and are usually used as food coloring (Bratinova 2015) , thus substantiating previous results. No clinical trial on the website ClinicalTrials.com or docking study directly reported the use of those compounds with the aim to cure or prevent viral infections, thus opening a new way of research.

Phycocyanins are mainly known for their great antioxidant properties (Zhou et al. 2005; Wu et al. 2016 ). An in vitro study showed that phycocyanobilin represents the main Arthrospira compound responsible for its antioxidant effect (Hirata et al. 2000) . Phycocyanobilin was shown to have a similar structure to biliverdin, an NADPH inhibitor (Zheng et al. 2013) . Phycocyanobilin is converted into phycocyanorubin by biliverdin reductase (Terry et al. 1993 ) with a structure close to bilirubin (McCarty 2007) . In other in vivo studies, both phycocyanin and phycocyanobilin were able to normalize oxidative stress markers and expression of components from NADPH pathway (Zheng et al. 2013) . McCarty et al. (2010) suggested NADPH inhibition as a way of treatment in influenza viral infection. Indeed, it has been identified as the main oxidative stress source in lung epithelial cells, participating in viral symptoms. Thus, phycobilins and especially phycocyanobilin present additional arguments to support their antiviral potential.

One clinical study focuses on dieckol. This trial in progress relies on safety issues of several phlorotannins including dieckol (NCT04335045). It should provide further information to confirm or infirm our in silico toxicity assessment. Drug interaction study provided high risk of inhibition of several cytochromes. No scientific paper studied specific inhibition of CYP1A2, 2C9, and 2C19 by dieckol. However, Sadeeshkumar et al. (2017) found that dieckol could reduce activities of liver cytochromes P450 after being previously increased with N-nitrosodiethylamine. Antiviral potential of dieckol has been studied in several papers. Kwon et al. (2013) showed that this phlorotannin was able to block viral entry of porcine epidemic diarrhea virus via hemagglutination and replication inhibitions, with an IC 50 of 14.6 μM. Kim (2013) reported an in vitro activity of phlorotannins, especially dieckol, against influenza virus neuraminidase with about 72% inhibition at 30 μg.L −1 . Other studies highlighted its in vitro antiviral activity against both human papilloma virus and murine norovirus (Eom et al. 2015; Kim and Kwak 2015) . Regarding coronaviruses, dieckol was submitted to an in vitro study wherein it showed an IC 50 of 2.7 μM against 3CLPro of SARS-CoV-1 . Moreover, in a recent docking paper using Autodock Vina, dieckol bonds the main protease (Mpro) of SARS-CoV-2 with a binding energy of − 12.0 kcal.mol −1 , a greater result than those obtained with Autodock Vina and SwissDock on Spike RBD (− 8.1 and − 8.375 kcal.mol −1 ) . Dieckol therefore arouses the interest of other researchers in the global context of Covid-19 and would be probably a more potent inhibitor against Mpro than Spike RBD but could be limited by its low bioavailability according to absorption prediction data.

Folic acid, also known as vitamin B 9 , is a largely used molecule, approved and under health claims by both FDA and EFSA (Electronic Code of Federal Regulation 2000; EFSA 2010). Unlike our absorption prediction results, EFSA reported a great bioavailability of different forms of folic acid (EFSA 2014). Indeed, folic acid is usually used during pregnancy to prevent the risk of neural tube defects (Greenberg et al. 2011) . It is also widely used as concomitant treatment with methotrexate to offset its side effects including hepatotoxicity. Indeed, methotrexate (an antirheumatic drug) tends to exert a competitive activity against folic acid, necessitating a substantial supplementary administration of this folic acid (Shea et al. 2013) . In this context, folic acid has been included in two clinical studies: one in HIV-infected patients taking both methotrexate and folic acid (Clinical trial no. NCT01949116) and a second one on hepatitis C-infected patients (Clinical trial no. NCT02150291). The latter studied how folic acid could act as prophylactic treatment in a viral infection to counteract side effects of the main treatment. It should be noted that predicted hepatotoxicity risk was not verified in literature. Coupled with the fact that folic acid is used to counteract the hepatotoxicity effects of methotrexate, this suggests hepatotoxicity risk is actually not relevant. Besides, folic acid provided strong results during docking assays, with binding energies of − 6.95 and − 10.35 kcal.mol −1 with Autodock Vina and SwissDock. These arguments are in favor of a safety and well-known consumption of folic acid which thus can be seriously considered as a candidate in the fight against Covid-19.

Rutin has already been studied for its potential as antiviral for SARS-CoV-2. Indeed, several docking papers mentioned its promising effect as inhibitor of some virus targets, especially against the main protease (Mpro). Das et al. (2020a) explored the efficiency of Calendula officinalis compounds as ligand for Mpro, using Autodock Vina. With a binding energy of − 8.8 kcal.mol −1 , rutin was the best molecule of the assay. Similar results have been obtained in the scientific literature, still with Autodock Vina, with − 9.0 kcal.mol −1 (Huynh et al. 2020) . Two other teams confirmed these results, within rutin remained the best tested compound among tens, including antivirals (Das et al. 2020b; Hu et al. 2020) . Only one study (pending for publication) focused on the interaction between rutin and ACE2/Spike complex. Among an FDAapproved molecule list, rutin was ranked in the best results. Moreover, molecular dynamic simulation (MD) showed its ability to destabilize the ACE2/Spike complex, thus consolidating our results (binding energies of − 7.2 and − 8.085 kcal.mol −1 with Autodock Vina and SwissDock) (https:// europepmc.org/article/ppr/ppr188118). Elsewhere, several studies showed in vitro effects of rutin against other viruses belonging to enterovirus, norovirus, or hepatitis virus families (Lin et al. 2012; Chéron et al. 2015; Parvez et al. 2019 . These data highlight the rise of rutin in the field of antivirals. Regarding hepatotoxicity prediction, literature infirmed this result since several studies showed a protective effect of rutin on hepatotoxicity such as paracetamol-induced hepatotoxicity (Janbaz et al. 2002; Reddy et al. 2017) . However, controversial studies underlined mutagenic activity of rutin and its main metabolite quercetin, corroborating in silico prediction of our study (Crebelli et al. 1987; Okamoto 2005) . If the conclusion of the debate is still unclear, it is yet necessary to use this compound with precautions and further safety assessment would be required to clarify this point.

At this day, no clinical study assesses antiviral abilities of β-carotene and few papers mention it. Nevertheless, this molecule has been included in a placebo-controlled study where it was given in supplementation in HIV patients to analyze enzymatic antioxidant system response. Like selenium, also tested in the experiments, positive conclusions have been made with an improvement of antioxidant enzyme activities (Delmas-Beauvieux et al. 1996) . More recently, a docking study has been conducted on interactions between this compound and hepatitis C virus (HCV). However, no interaction was highlighted with the active site of the HCV helicase (Kaneez et al. 2014) . In SARS-CoV-2 context, Lianhua Qingwen Capsule (LQC) showed some therapeutic effects on Covid-19 patients . In view to understand the molecular mechanism of such activity, docking analysis has been executed with Autodock Vina ). Among the six active compounds of LQC with great binding energies, a − 6.5 kcal.mol −1 binding score was predicted for β-carotene, showing its ability to bind with the active pocket of Akt1. This human kinase has been suggested as a potential therapeutic target for Covid-19 treatment. Indeed, its overexpression could promote viral protein synthesis and thus help in virus spreading. Regarding Spike RBD interaction, βcarotene presented greater binding energies than with Akt1, with − 7.15 and − 8275 kcal.mol −1 , suggesting a complementary activity of this molecule on the two targets. As folic acid previously mentioned, β-carotene is well known as a nutritive molecule and precursor of vitamin A. Used in plenty of clinical trials and widely commercialized as a dietary supplement, β-carotene consumption is generally assumed as safe. Nevertheless, EFSA highlighted some studies with preoccupying side effects on the lung, especially in heavy smokers (EFSA 2006) . Moreover, according to our experiments and literature, this compound shows a low bioavailability due to its long carbon chain (EFSA 2006) . Thus, an important dose would be required for this compound to act as antiviral compound and inhibit SARS-CoV-2 Spike protein.

Astaxanthin is not a well-documented substance regarding antiviral activities since no clinical study deals with this topic. Indeed, astaxanthin is mainly known for its antioxidant and anti-inflammatory properties. In this way, it has been included in two clinical trials (NCT03945526 and NCT03906825). The first one is related to effects of astaxanthin on antioxidant parameters in stroke patients, whereas the second one relied on effects of a dietary supplement containing astaxanthin on inflammation markers (Birudaraju et al. 2020) . Talukdar et al. (2020a) brought arguments to consider this compound as a serious natural molecule to alleviate the risk of cytokine storm in some Covid-19 patients thanks to its anti-inflammatory properties. Elsewhere, in a study conducted in vitro on human papilloma virus (HPV), astaxanthin prevented HPV L1 protein (capsid protein) by binding on human sperm membrane. Such antiviral activity could protect men from idiopathic infertility caused by HPV (Donà et al. 2018) . In silico toxicity study highlighted a high risk of reprotoxicity, but this has not been verified in literature. Indeed, acute and subchronic toxicity studies showed no adverse effects of astaxanthin (Stewart et al. 2008; Katsumata et al. 2014) and no reprotoxicity study has been conducted yet to our knowledge. Despite interesting binding energies (− 7.05 and − 8.195 kcal.mol −1 ), astaxanthin is, such as β-carotene, a low bioavailable molecule (Mercke Odeberg et al. 2003) . For the same reasons previously mentioned, these data therefore suggest that astaxanthin would be more effective as an anti-inflammatory agent than a direct antiviral compound in Covid-19 treatment.

Many articles referred to docking techniques in view to assess the potential of a plethora of molecules as inhibitors against different SARS-CoV-2 targets (Elfiky 2020; Quimque et al. 2020; Ton et al. 2020 (475), none of the two others interacted with key binding residues involved in this critical bind. Another publication presented hesperidin as the only potent inhibitor of Spike/ACE2 interface despite a quite large list of compounds with high binding energies . This result was consolidated by another article pending publication where hesperidin has been tested on a closed state of Spike protein (https://www.preprints.org/manuscript/202004.0102/v3). The ligand obtained a binding energy of − 10.4 kcal.mol −1 with Autodock Vina software. Nevertheless, this compound, as mentioned by the authors, is characterized by a low oral bioavailability and then does not constitute the ideal candidate for a new antiviral development.

SARS-CoV-2 can recognize several cellular receptors such as ACE2 or HSPA5 (Heat Shock Protein A5) also known as GRP78. Souza et al. (2020) found several compounds among estrogenic compounds capable of competing with SARS-CoV-2 on the substrate binding domain of HSPA5. Autodock Vina provided great binding energies of about − 7.8 kcal.mol −1 for some compounds. Nevertheless, phytoestrogens are controversial molecules, suspected of being involved in some estrogen-dependent cancers (Kwack et al. 2009; Touillaud et al. 2019) . Moreover, consequences of human receptor inhibition are unknown, which is why we preferred working on viral protein instead of human targets.

Glycyrrhizic acid has been elsewhere identified after a docking assay with Vina on Spike protein, with a − 9.2 kcal.mol −1 binding energy (Sinha et al. 2020b) . Despite this promising result, analysis of binding residues showed that engaged amino acids were not located on key binding residues. Indeed, they were positioned on HIS625, GLU298, ARG319, GLN321, VAL320, TYR343, THR323, ILE624, and ASN317 and thus do not belong to the Spike RBD (i.e., residues 333-527). This suggests that such interactions between glycyrrhizic acid and Spike may not impede the link between Spike and ACE2.

Spike RBD has been isolated from the Spike/ACE2 complex for docking studies in two papers, thus proposing a similar method to this publication. The first one analyzed essential oil components and no compounds showed binding energy greater than − 5.5 kcal.mol −1 . Moreover, few H-bonds were predicted and no key binding residues were engaged (Kulkarni et al. 2020) . Schrödinger software was used in the second study, focused on Indian ginger compounds (Chikhale et al. 2020) . The authors highlighted the potential of QGRG (quercetin-3-o-galactosyl-rhamnosyl-glucoside) for Spike RBD inhibition, with a − 9.25 kcal.mol −1 binding score. This compound presented three H-bonds with Spike RBD, located on GLU340, GLN245, and VAL339 which are not among the key binding residues, suggesting better candidates could be found for Spike RBD inhibition. During cryoelectron microscopy analysis, Toelzer et al. (2020) found that linoleic acid was bound to the Spike RBD in the closed conformation. The authors hypothesized that the fatty acid could stabilize the closed state of RBD and thus limiting ACE2 binding for which open state is required. In the same study, the antiviral remdesivir was able to inhibit RNA polymerase of SARS-CoV-2 in vitro. Because of possible adverse effects, WHO does not recommend conducting clinical trials on hospitalized patients who are weakened. However, clinical trials are still encouraged in populations with non-severe symptoms to identify effective antiviral therapy. As such, remdesivir has been included in several clinical trials (NCT04365725, NCT04582266, NCT04431453). To analyze the potential of linoleic acid as inhibitor of Spike RBD in an open state as well as the impact of an important antiviral molecule, both compounds were added to our docking study. Results revealed that linoleic acid weakly bonded with Spike RBD, with binding energies of − 4.85 and − 7.025 kcal.mol −1 via Autodock Vina and SwissDock (Table 8 Supplemental Data). This suggests the exposition of key binding residues is modified between the two conformational states and is therefore no more favorable for linoleic acid binding. This corroborates observations of Toelzer since no fatty acid has been recovered in the open state during their study. Regarding remdesivir docking, moderate binding affinities were obtained with respectively − 6.9 and − 8.02 kcal.mol −1 with Autodock Vina and SwissDock (Table 8 Supplemental Data). Chikhale et al. (2020) also included remdesivir in their docking study with Schrodinger software. Similar to our results, it presented a low/moderate binding energy of − 4.65 kcal.mol −1 , thus suggesting that its antiviral activity is not based on Spike RBD inhibition. P a v et an n i n C 1 ha s b e en c o n f i r m e d t w i c e b y docking assays as potential SARS-CoV-2 antiviral and has been thus chosen as positive control for the present study. This compound has been highlighted with MPro bonding abilities (− 9.5 kcal.mol −1 with Autodock Vina) in a study on phytochemicals from Indian cuisine. Rajendran et al. (2020) . Another study presented pavetannin C1 as a great Spike inhibitor (Prasanth et al. 2020) . Docking assay with Autodock 4.0 provided binding energy of − 11.1 kcal.mol −1 and different interactions were predicted with Spike RBD region. Indeed, four hydrophobic interactions engaged THR347, TRP349, GLU398, and TYR510, as well as nine H-bonds located on residues SER44, HIS345, ALA348, ASP350, GLU375, GLY395, ASN397, ARG514, and TYR515. However, none of the detected interactions engaged key binding residues. Results obtained with our experiments provided binding energies of − 7.3 and − 8.52 kcal.mol −1 with Autodock Vina and SwissDock and allowed to select this compound among the best candidates for antiviral development, thus corroborating literature results. Moreover, two interactions including one H-bond was predicted on key binding residues (GLN493 and GLY496). The slight variation of binding energies with our study could be explained in part by the use of different docking programs. Moreover, in Prasanth's work, docking has been executed with the same initial PDB file: 6LZG as in our work, but Spike RBD has not been isolated and thus remained complexed with human ACE2. This could explain the difference between predicted interactions. Pavetannin C1 could not bind directly key binding residues as they were still bonded with human ACE2. However, despite promising binding scores, little is known about this molecule. Indeed, no toxicity data would allow us to infirm or confirm toxicity prediction, especially regarding reprotoxicity and cytochrome interactions. It would be therefore wise to promote molecules on which with have more hindsight.

On the basis of these in silico results, our work provides a list of four reliable compounds with potential antiviral activity via Spike RBD inhibition (phycoerythrobilin, phycocyanobilin, phycourobilin, and folic acid). Interactions such as H-bonds with RBD key binding residues constitute the strength of these results since such bonds could destabilize or impede the link between SARS-CoV-2 and one of its human receptors, ACE2. It will be more interesting to follow the future work that will deal with the interaction of our four molecules in viral episodes and further in silico studies will improve and guide bioactive molecule selection for antiviral compound research.

Arthrospira is one of the most well-known genera in the world of microalgae. In several articles reviewed here, Arthrospira or Arthrospira extracts presented significant antiviral activities against many virus families. In the context of the Covid-19 pandemic, we were interested in the potential of this microalgae to fight against SARS-CoV-2. The present study considered docking techniques on a selection of 48 Arthrospira components which were submitted to two different docking software (Autodock Vina and SwissDock). The main purpose of the present paper was to test their ability to inhibit Spike protein, a structural protein of the virus. In a second time, in silico toxicity analysis has been conducted to assess the safety of selected compounds.

After docking assays, eight Arthrospira molecules have been identified for their great aptitude to link with Spike RBD. Among them, four provided good arguments to be considered for an oral administration according to absorption prediction and existing literature. These compounds were namely phycoerythrobilin, phycocyanobilin, phycourobilin, and folic acid. The other compounds have been discarded due to their low oral bioavailability or potential toxicity risks. Consensus docking and redocking methods allowed us to obtain reliable results by eliminating false-positive results. Our results were therefore consistent with those currently found in the literature.

Our study provides promising arguments for the inclusion of these four compounds in further studies to assess their ability to compete with the SARS-CoV-2/ACE2 complex both in vitro and in vivo.

Antiviral activity of freshwater algae

Adaptation of light-harvesting systems of Arthrospira platensis to light conditions, probed by time-resolved fluorescence spectroscopy

H++ 3.0: automating pK prediction and the preparation of biomolecular structures for atomistic molecular modeling and simulations

Inhibition of HIV-1 replication by an aqueous extract of Spirulina platensis (Arthrospira platensis)

Potential of Spirulina platensis as a nutritional supplement in malnourished HIV-infected adults in Sub-Saharan Africa: a randomised, single-blind study

Activation of the SARS coronavirus spike protein via sequential proteolytic cleavage at two distinct sites

A combined effect of cavacurcumin, eicosapentaenoic acid (omega-3s), astaxanthin and gamma-linoleic acid (omega-6) (CEAG) in healthy volunteers -a randomized, double-blind, placebo-controlled study

Proteolytic activation of influenza viruses by serine proteases TMPRSS2 and HAT from human airway epithelium

Provision of scientific and technical support with respect to the classification of extracts/concentrates with colouring properties either as food colours (food additives falling under Regulation (EC) No 1333/2008) or colouring foods

Repurposing of rutin for the inhibition of norovirus replication

Sars-cov-2 host entry and replication inhibitors from Indian ginseng: an in-silico approach

Anti HSV-1 activity of sulphoquinovosyl diacylglycerol isolated from Spirulina platensis

Type II transmembrane serine proteases in cancer and viral infections

Urinary and faecal mutagenicity in Sprague-Dawley rats dosed with the food mutagens quercetin and rutin

SwissADME: a free web tool to evaluate pharmacokinetics, drug-likeness and medicinal chemistry friendliness of small molecules

Small-molecule library screening by docking with PyRx

In-Silico approach for identification of effective and stable inhibitors for COVID-19 main protease (Mpro) from flavonoid based phytochemical constituents of Calendula officinalis

An investigation into the identification of potential inhibitors of SARS-CoV-2 main protease using molecular docking study

Impact of dietary Arthrospira (Spirulina) biomass consumption on human health: main health targets and systematic review

The enzymatic antioxidant system in blood and glutathione status in human immunodeficiency virus (HIV)-infected patients: effects of supplementation with selenium or beta-carotene

Antiviral activity of five filamentous cyanobacteria against coxsackievirus B3 and rotavirus

Astaxanthin prevents human papillomavirus L1 protein binding in human sperm membranes

A novel angiotensin-converting enzyme-related carboxypeptidase (ACE2) converts angiotensin I

Tolerable upper intake levels for vitamins and minerals

Scientific Opinion on the substantiation of health claims related to folate and contribution to normal psychological functions

Scientific Opinion on Dietary Reference Values for folate

In vitro antiviral and antimicrobial activities of Spirulina platensis extract

Ribavirin, remdesivir, sofosbuvir, galidesivir, and tenofovir against SARS-CoV-2 RNA dependent RNA polymerase (RdRp): a molecular docking study

In vitro antiviral activity of dieckol and phlorofucofuroeckol-A isolated from edible brown alga Eisenia bicyclis against murine norovirus

Computational protein-ligand docking and virtual drug screening with the AutoDock suite

Putative inhibitors of SARS-CoV-2 main protease from a library of marine natural products: a virtual screening and molecular modeling study

Immunological studies in Egyptian thalassemic children infected with hepatitis C virus

H++: a server for estimating pK a s and adding missing hydrogens to macromolecules

Folic acid supplementation and pregnancy: more than just neural tube defect prevention

Fast docking using the CHARMM force field with EADock DSS

SwissDock, a protein-small molecule docking web service based on EADock DSS

Cryoelectron microscopy structures of the SARS-CoV spike glycoprotein reveal a prerequisite conformational state for receptor binding

Tissue distribution of ACE2 protein, the functional receptor for SARS coronavirus. A first step in understanding SARS pathogenesis

Antiviral activity of Spirulina maxima against herpes simplex virus type 2

Antiviral and antimicrobial activities of Spirulina platensis

Antioxidant activities of phycocyanobilin prepared from Spirulina platensis

SARS-CoV-2 cell entry depends on ACE2 and TMPRSS2 and is blocked by a clinically proven protease inhibitor

Possible SARS-coronavirus 2 inhibitor revealed by simulated molecular docking to viral main protease and host toll-like receptor

Scoring functions and their evaluation methods for protein-ligand docking: recent advances and future directions

Structure-based lead optimization of herbal medicine rutin for inhibiting SARS-CoV-2's main protease

Protective effect of rutin on paracetamol-and CCl 4 -induced hepatotoxicity in rodents

A Spirulina maxima-derived peptide inhibits HIV-1 infection in a human T cell line MT4

Docking studies of Pakistani HCV NS3 helicase: a possible antiviral drug target

A sub-chronic toxicity evaluation of a natural astaxanthin-rich carotenoid extract of Paracoccus carotinifaciens in rats

Antiviral phlorotannin from Eisenia bicyclis against human papilloma virus in vitro

Inhibitory activities of microalgal extracts against Epstein-Barr virus DNA release from lymphoblastoid cells

Computational evaluation of major components from plant essential oils as potent inhibitors of SARS-CoV-2 spike protein

Risk assessment of soybean-based phytoestrogens

In vitro antiviral activity of phlorotannins isolated from Ecklonia cava against porcine epidemic diarrhea coronavirus infection and hemagglutination

Structure of the SARS-CoV-2 spike receptor-binding domain bound to the ACE2 receptor

Effect of partial desulfation and oversulfation of sodium spirulan on the potency of anti-herpetic activities

Structure, function, and evolution of coronavirus spike proteins

Structure of SARS coronavirus spike receptor-binding domain complexed with receptor

Lianhua Qingwen prescription for Coronavirus disease 2019 (COVID-19) treatment: advances and prospects

Fisetin and rutin as 3C protease inhibitors of enterovirus A71

Characterization of a highly conserved domain within the severe acute respiratory syndrome coronavirus spike protein S2 domain with characteristics of a viral fusion peptide

Transcriptional regulation and expression of CYP3A4 in hepatocytes

Clinical potential of Spirulina as a source of phycocyanobilin

Practical strategies for targeting NF-kappaB and NADPH oxidase may improve survival during lethal influenza epidemics

Oral bioavailability of the antioxidant astaxanthin in humans is enhanced by incorporation of lipid based formulations

Middle East respiratory syndrome coronavirus infection is inhibited by griffithsin

The potential of microalgae for the production of bioactive molecules of pharmaceutical interest

A simple clustering algorithm can be accurate enough for use in calculations of pKs in macromolecules

Impact of daily supplementation of Spirulina platensis on the immune system of naïve HIV-1 patients in Cameroon: a 12-months single blind, randomized, multicenter trial

Autodock Vina adopts more accurate binding poses but Autodock4 forms better binding affinity

Open Babel: an open chemical toolbox

Safety of quercetin for clinical application (review)

Characterization of spike glycoprotein of SARS-CoV-2 on virus entry and its immune cross-reactivity with SARS-CoV

Dieckol, a SARS-CoV 3CLpro inhibitor, isolated from the edible brown algae Ecklonia cava

Plant-derived antiviral drugs as novel hepatitis B virus inhibitors: cell culture and molecular docking study

Algae: a potential source to prevent and cure the novel coronavirus -a review

In silico identification of potential inhibitors from Cinnamon against main protease and spike glycoprotein of SARS CoV-2

Virtual screening-driven drug discovery of SARS-CoV2 enzyme inhibitors targeting viral attachment, replication, post-translational modification and host immunity evasion infection mechanisms

Anionic polysaccharides from phototrophic microorganisms exhibit antiviral activities to Vaccinia virus

Arul Selvaraj RC, Provaznik I (2020) In silico screening and molecular dynamics of phytochemicals from Indian cuisine against SARS-CoV-2 M Pro

Automated minimization of steric clashes in protein structures

In vitro P-glycoprotein inhibition assays for assessment of clinical drug interaction potential of new drug candidates: a recommendation for probe substrates

Antiviral activity of Arthrospiraderived spirulan-like substances

Protective effect of rutin in comparison to silymarin against induced hepatotoxicity in rats

Ultrasound-assisted extraction of phycobiliproteins from Spirulina (Arthrospira) platensis using protic ionic liquids as solvent

Modulatory efficacy of dieckol on xenobiotic-metabolizing enzymes, cell proliferation, apoptosis, invasion and angiogenesis during NDEA-induced rat hepatocarcinogenesis

Salt stress enhancement of antioxidant and antiviral efficiency of Spirulina platensis

Structural basis of receptor recognition by SARS-CoV-2

Antiherpetic efficacy of aqueous extracts of the cyanobacterium Arthrospira fusiformis from Chad

Folic acid and folinic acid for reducing side effects in patients receiving methotrexate for rheumatoid arthritis

A transmembrane serine protease is linked to the severe acute respiratory syndrome coronavirus receptor and activates virus entry

An in-silico evaluation of different saikosaponins for their potency against SARS-CoV-2 using NSP15 and fusion spike glycoprotein as targets

Identification of bioactive compounds from Glycyrrhiza glabra as possible inhibitor of SARS-CoV-2 spike glycoprotein and non-structural protein-15: a pharmacoinformatics study

A molecular docking study revealed that synthetic peptides induced conformational changes in the structure of SARS-CoV-2 spike glycoprotein, disrupting the interaction with human ACE2 receptor

Safety assessment of astaxanthin-rich microalgae biomass: Acute and subchronic toxicity studies in rats

Potential of natural astaxanthin in alleviating the risk of cytokine storm in COVID-19

COVID-19: Potential of microalgae derived natural astaxanthin as adjuntive supplement in alleviating cytokine storm

Dietary algae and HIV/AIDS: proof of concept clinical data

Inactivation of phytochromeand phycobiliprotein-chromophore precursors by rat liver biliverdin reductase

Free fatty acid binding pocket in the locked structure of SARS-CoV-2 spike protein

Rapid identification of potential inhibitors of SARS-CoV-2 main protease by deep docking of 1.3 billion compounds

Use of dietary supplements containing soy isoflavones and breast cancer risk among women aged >50 y: a prospective study

AutoDock Vina: improving the speed and accuracy of docking with a new scoring function, efficient optimization, and multithreading

Tectonic conformational changes of a coronavirus spike glycoprotein promote membrane fusion

Structure, function, and antigenicity of the SARS-CoV-2 spike glycoprotein

Structural and functional basis of SARS-CoV-2 entry by using human ACE2

The effect of Arthrospira platensis capsules on CD4 T-cells and antioxidative capacity in a randomized pilot study of adult women infected with human immunodeficiency virus not under HAART in Yaoundé

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

Stability and antioxidant activity of food-grade phycocyanin isolated from Spirulina platensis

Analysis of therapeutic targets for SARS-CoV-2 and discovery of potential drugs by computational methods

Network pharmacology and molecular docking analyses on Lianhua Qingwen capsule indicate Akt1 is a potential target to treat and prevent COVID-19

Spirulina platensis versus silymarin in the treatment of chronic hepatitis C virus infection. A pilot randomized, comparative clinical trial

Structural basis for the recognition of SARS-CoV-2 by full-length human ACE2

admetSAR 2.0: web-service for prediction and optimization of chemical ADMET properties

Cryo-EM structures of MERS-CoV and SARS-CoV spike glycoproteins reveal the dynamic receptor binding domains

Phycocyanin and phycocyanobilin from Spirulina platensis protect against diabetic nephropathy by inhibiting oxidative stress

) factors that affect antioxidant activity of C-phycocyanins from Spirulina platensis

SwissParam: a fast force field generation tool for small organic molecules

Publisher's note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations