key: cord-1055143-z59idt2b authors: Wong, Fai-Chu; Ong, Joe-Hui; Chai, Tsun-Thai title: SARS-CoV-2 spike protein-, main protease- and papain-like-protease-targeting peptides from seed proteins following gastrointestinal digestion: An in silico study date: 2021-02-28 journal: Phytomedicine Plus DOI: 10.1016/j.phyplu.2020.100016 sha: 0ff61fc838b562405c89b65df34912f86f21ec31 doc_id: 1055143 cord_uid: z59idt2b Background The anti-COVID-19 potential of phytochemicals was investigated in numerous studies, but efficacy of peptides released by seed proteins upon gastrointestinal (GI) digestion is underexplored. Purpose This study investigated whether multi-target anti-COVID-19 peptides could be released from edible seeds following GI digestion, by using in silico and molecular docking approaches. Methods Nineteen seed storage proteins from Chenopodium quinoa (quinoa), Sesamum indicum (sesame), Brassica napus (rape), Helianthus annuus (sunflower) and Cucurbita maxima (pumpkin) seeds were subjected to in silico GI digestion, in order to detect the released peptides with high GI absorption that concurrently target the spike protein, main protease and papain-like protease of SARS-CoV-2. Results Molecular docking study revealed that 36 peptides with high GI absorption, out of the 1593 peptides released from seed proteins, could bind to the binding or catalytic sites of the spike protein, main protease and papain-like protease of SARS-CoV-2, after GI digestion. Among the five seeds, quinoa was predicted to release the largest number (27) of multi-target peptides. When compared with PIY (Pro-Ile-Tyr), a high-GI-absorption fragment released from a potential anti-COVID-19 peptide, pumpkin-derived peptide PW (Pro-Trp) could bind more strongly to SARS-CoV-2 spike protein. PW was superior to some previously reported anti-SARS-CoV-2 phytochemicals when binding affinities towards the three viral targets were compared. Conclusion Edible seeds are a potential source of anti-COVID-19 peptides upon GI digestion, hence they should be considered as an alternative to assist in the treatment and management of COVID-19. The pandemic of the coronavirus disease 2019 , caused by a novel pathogenic virus named as Severe Acute Respiratory Syndrome Coronavirus type 2 (SARS-CoV-2), has escalated globally. The great pandemic potential of SARS-CoV-2 could be attributed to its high human-to-human transmissible efficiency ( Zehra et al., 2020 ) . Patients affected by COVID-19 most commonly experience cough, runny nose, fever and headache, although some suffer from serious complications such as severe pneumonia, bronchitis, multi-organ failure, respiratory illness, and even deaths ( Zehra et al., 2020 ) . SARS-CoV-2 is a betacoronavirus whose genome codes for four major structural proteins, namely spike, envelope, membrane and nucleocapsid proteins, as well as non-structural proteins, for instance the main protease (M pro ) and papain-like protease (PL pro ) ( Zehra et al., 2020 ) . The interaction between the receptor-binding domain (RBD) of SARS-CoV-2 spike glyco-could serve as specific targets in the formulation of inhibition strategies and pharmacological treatments against SARS-CoV-2 ( Ortega et al., 2020 ) . In response to the COVID-19 pandemic, there is a need to develop countermeasures urgently. This has fueled current interest among researchers to adopt molecular docking and other in silico approaches to evaluate the potential of existing antiviral drugs and bioactive natural products in targeting SARS-CoV-2 ( Fakih, 2020 ) . Molecular docking has been widely exploited as structure-based drug discovery tool owing to its quick, low-cost and precise prediction of the biological activity and interactions of ligands with target molecules. With this approach, information on the bond types, binding affinity and interactions between two molecules can be predicted from the best-matching conformation ( Tao et al., 2020 ) . To date, molecular docking studies have identified some natural and synthetic inhibitors that could potentially block the binding between SARS-CoV-2 and hACE2 and potentially inhibit the catalytic activities of SARS-CoV-2 M pro and PL pro , as presented in Tables S1-S6. Notwithstanding the foregoing, studies focusing on food-derived and other natural peptides, specifically those that could concurrently target the spike glycoproteins, M pro and PL pro of SARS-CoV-2 are scarce. Food peptides exhibit diverse bioactivities, including antiviral activity ( Agarwal and Gabrani, 2020 ) . The potential of food-derived peptides in the development of nutraceuticals and functional foods is wellrecognized ( Chai et al., 2017 ; Wong et al., 2020b ) . Peptides are also preferred as therapeutic agents due to their higher specificity, efficiency and biodegradability, as well as lower molecular weight and toxicity ( Haggag et al., 2018 ) . Considering the aforementioned, in this study, we investigated whether food-derived peptides could serve as multitarget inhibitors against host cell entry and viral replication of SARS-CoV-2. Specifically, we explored the potential release of such multifunctional peptides from the seed storage proteins of Cucurbita maxima L. (Cucurbitaceae) (pumpkin), Chenopodium quinoa Willd. (Amaranthaceae) (quinoa), Brassica napus L. (Brassicaceae) (rape), Helianthus annuus L. (Asteraceae) (sunflower) and Sesamum indicum L. (Pedaliaceae) (sesame) seeds by using in silico gastrointestinal (GI) digestion, followed by molecular docking analysis. Seed storage proteins are the major proteins in the five seeds and are likely to be the main sources of bioactive peptides upon oral consumption and GI digestion. The nutritional values and the bioactivity of peptides derived from the five seeds have been reported ( Bao et al., 2020 ; Guo et al., 2020 ; Han et al., 2019 ; Innocent-Ukachi, 2019 ; Vilcacundo et al., 2018 ; Wang et al., 2020b ) . Currently, there is still no effective approved treatment or vaccine that can be used to tackle or prevent COVID-19. The outcome of this research could relevant information that could be used by medical professionals for formulating strategy to assist in the treatment and/or management of COVID-19. Importantly, this research could serve as the basis for future studies to elucidate the health-promoting effects of seed peptides upon oral ingestion, especially in relation to coronavirus diseases. Nineteen seed storage proteins found in pumpkin, quinoa, rape, sunflower and sesame seeds ( Guo et al., 2020 ; Han et al., 2019 ) , were targeted in this study. Quinoa seed protein sequences were retrieved from NCBI ( https://www.ncbi.nlm.nih.gov ) ( Clark et al., 2015 ) on 5 August 2020. The protein sequences of the other four seeds were downloaded from the UniProt web server ( http://www.uniprot.org ) ( The UniProt Consortium, 2019 ) on 5 August 2020. In silico GI digestion of peptides and the screening of GI absorption were performed by using BIOPEP-UWM ( http://www.uwm.edu.pl/biochemia/index.php/en/biopep ) ( Minkiewicz et al., 2019 ) and SwissADME ( http://www.swissadme.ch ) ( Daina et al., 2017 ) web servers, respectively, as previously described ( Wong et al., 2020a ) . Besides, phytochemicals used as references in molecular docking studies were also screened for GI absorption, entering their canonical Simplified Molecular-Input Line-Entry System (SMILES) from PubChem ( https://pubchem.ncbi.nlm.nih.gov ) ( Kim et al., 2016 ) as input to be analyzed by SwissADME. Following screening, peptides and reference phytochemicals predicted to have high GI absorption (Tables S4-S6) were taken for molecular docking analysis. The crystal structures of SARS-CoV-2 proteins were obtained from the Protein Data Bank ( http://www.rcsb.org/pdb ) ( Berman et al., 2000 ; Burley et al., 2018 ) . Three SARS-CoV-2 viral proteins were chosen for this study, namely spike glycoprotein RBD (PDB ID: 6LZG) ( Wang et al., 2020a ) , M pro (PDB ID: 6LU7) ( Jin et al., 2020 ) and PL pro (PDB ID: 6W9C) ( Osipiuk et al., 2020 ) . Proteins and ligands were separated and prepared for molecular docking by using BIOVIA Discovery Studio Visualizer (BIOVIA, Dassault Systèmes, BIOVIA Discovery Studio Visualizer, Version 20.1.0.192, San Diego: Dassault Systèmes, 2020), as previously reported ( Wong et al., 2020a ) . Crystallographic water molecules were removed when the target protein structures were prepared for molecular docking analysis. There is no evidence in the literature that such water molecules are part of the viral proteins structure and tightly bound to viral proteins. Hence, maintaining the water molecules in the target protein structures for molecular docking simulation is not justified. Docking energy scores of high-GI-absorption peptides, types of seeds they were released from, and capability to bind to key binding residues on RBD as well as to key catalytic residues of Mpro and PLpro. Docking energy scores Seeds Binding to key residues a M pro PL pro Spike protein RBD Pumpkin Quinoa Rape Sunflower Sesame a Peptides predicted to bind to at least one key binding residue on RBD as well as to at least one key catalytic residue of M pro and PL pro . To perform molecular docking on the HPEPDOCK server ( http://huanglab.phys.hust.edu.cn/hpepdock ) ( Zhou et al., 2018 ) , the three prepared viral proteins were uploaded as receptor inputs. The prepared hACE2 and N3 were uploaded as binding site references for RBD and M pro respectively. For PL pro , three receptor binding residues Cys111, His272 and Asp286 were entered in the docking option for HPEPDOCK analysis ( Wlodawer et al., 2020 ) . High-GI-absorption seed peptides released from in silico GI digestion and selected reference peptides that were previously reported to bind to SARS-CoV-2 viral proteins, where 3D structures are not available, were entered in the FASTA format as peptide input. Tables S1-S3 list the reference peptides used and sources of their crystal structures, where available. The peptide model giving the lowest (most negative) docking score as indicated by the HPEPDOCK analysis was chosen to combine with the respective SARS-CoV-2 viral protein to form a docked model. The protein-peptide interactions of the docked model were analyzed by using LigPlot + v .2.2 ( Laskowski and Swindells, 2011 ; Wallace et al., 1995 ) . The binding poses of selected peptides on the surface of the viral proteins were visualized by using BIOVIA Discovery Studio Visualizer. Following HPEPDOCK analysis, the high-GI-absorption seed peptide with the lowest docking energy score when docked on three SARS-CoV-2 target proteins was further compared with selected phytochemicals (Tables S4-S6 ). The reference phytochemicals were chosen from those reported to bind to the three viral proteins and also predicted in this study to have high GI absorption. Molecular docking analysis of the selected seed peptide and reference phytochemicals was performed by using the AutoDock Vina option in PyRx 0.8 ( Dallakyan and Olson, 2015 ; Trott and Olson, 2010 ) as previously described ( Wong et al., 2020a ) . Docking was carried out for the spike protein RBD, with the grid box set to center x-, y-, z-axis values of − 40.34, 30.85, 6.61 with dimension of 35 Å x 50 Å x 30 Å, encompassing 21 binding site residues (Lys417, Gly446, Tyr449, Tyr453, Leu455, Phe456, Tyr473, Ala475, Gly476, Glu484, Phe486, Asn487, Tyr489, Phe490, Gln493, Gly496, Gln498, Thr500, Asn501, Gly502 and Tyr505) ( Wang et al., 2020a ) and three other RBD residues (Ser459, Gln474 and Pro499) reported to be important for RBD-hACE2 binding ( Yi et al., 2020 ) . For M pro , grid box was set with center x-, y-, z-axis values of − 10.71, 12.41, 68.83 with dimension of 15 Å x 25 Å x 15 Å, surrounding the catalytic-dyad His41 and Cys145 ( Ortega et al., 2020 ) . For the docking of PL pro , the grid box was set with center x-, y-, z-axis values of − 37.19, 24.28, 37.99 with each dimension of 20 Å, encircling the catalytic-triad Cys111, His272 and Asp286 ( Wlodawer et al., 2020 ) . The resultant binding affinities of the peptide and reference phytochemicals were tabulated. For binding affinities, data are presented as mean ± standard errors based on analysis in triplicates. Statistical analysis was performed by using the SAS University Edition Software (Version 9.4) ( https://www.sas.com/en_us/software/university-edition.html ). Pvalue of less than 0.05 was considered as statistically significant according to Tukey's test. In silico GI digestion generated a total of 1593 peptides fragments from 19 seed storage proteins taken from the five seed samples. The largest number of peptide fragments (801 fragments) were released from quinoa seed storage proteins, whereas the smallest number of peptide fragments (194 fragments) were from rape seed ( Table 1 ) . Regardless of seed type, more peptide fragments were released from the globulintype proteins than from the albumin-type. For example, the 13S and 11S globulins of the quinoa seed collectively released about 93% of all GIreleased peptide fragments. Our finding agrees with the observation that most peptide fragments found in the gastroduodenal digest of quinoa proteins were derived from 11S globulin of the seed ( Vilcacundo et al., 2018 ) . The distribution of peptide length in the in silico GI digestionreleased fragments varied among the five seeds. For each seed storage protein listed in Table 1 , the total pool of 2-and 3-residue peptides accounts for about 27-50% of the total number of peptides released. Such short peptides are of great interest as they may be absorbed intact through the epithelium of GI tract, entering the blood stream and producing biological effects at the tissue level ( Wang et al., 2019 ) . Thus, our findings suggest the five seeds analyzed in this study are likely to release a significant portion of potentially bioavailable peptides upon GI digestion. Besides resistance to GI digestion, GI absorption is another determinant of the bioavailability of short peptides. As predicted by the SwissADME tool, the greatest number of high-GI-absorption peptides was released from quinoa (108 fragments), in contrast with the smallest number from sunflower seeds. For all five seeds, the globulin-type seed storage proteins are the main contributors of high-GI-absorption peptides ( Table 1 ) . Together, our finding suggests that the quinoa seed, in particular its globulin-type storage proteins, is a promising source of potentially bioavailable peptides, which deserves more research attention in future. The same peptide sequences might be released from different proteins of same seed, or from proteins of different seeds. Hence, the 199 # Residues in bold are the catalytic site residues of M pro and PL pro , or key binding sites on the spike protein RBD. The underlined residues are other RBD residues that also bind to hACE2. high-GI-absorption peptides in Table 1 can be consolidated into a set of 56 unique sequences. Molecular docking analysis by HPEPDOCK revealed that the docking energy scores of the 56 peptides range from − 49.066 to − 144.652, − 36.734 to − 132.325, and − 38.642 to − 128.975, corresponding to viral proteins M pro , PL pro and spike glycoprotein RBD, respectively ( Table 2 ). In general, lower (more negative) scores suggest more stable binding between a peptide and its target protein (receptor). Based on docking energy scores, the five most promising peptides, in descending order, are PW, PSF, IVF, PGF and PIL, when docked with M pro ; PW, PSF, PGF, IVF and PY when docked with PL pro ; and PW, PSF, PGF, APY and IVF when docked with spike glycoprotein RBD ( Table 2 ) . Overall, our analysis implies that all five seeds analyzed in this study could potentially release not only high-GI-absorption peptides, but some such peptides may concurrently target the M pro , PL pro and spike glycoprotein RBD of SARS-CoV-2. We also analyzed potential anti-SARS-CoV-2 peptides reported in the literature for comparison and found that unlike PW, the bivalve mollusk Mizugopecten yessoensis myosin-derived QRPR (Gln-Arg-Pro-Arg) and LPIY (Leu-Pro-Ile-Tyr) were predicted to be susceptible to GI digestion by BIOPEP-UWM (data not shown). Fragments released from GI-digested QRPR were predicted to have low GI absorption by Swis-sADME. Only PIY (Pro-Ile-Tyr) released from GI-digested LPIY was predicted to exhibit high GI absorption. The docking energy score computed for PIY was lower than that of PW, suggesting stronger M pro -binding potential of PIY relative to PW (data not shown). Our analysis found that for binding to PL pro , bacterial azurin-derived p18 and p28 peptides had similar or superior docking scores relative to the seed peptides in our study. However, p18 and p28 were predicted to be susceptible to GI degradation and absorbed by the GI tract poorly (data not shown). While M. yessoensis myosin-derived CSNAIPEL (Cys-Ser-Asn-Ala-Ile-Pro-Glu-Leu) and LPIY could potentially bind to spike protein RBD, they were predicted to be susceptible to GI digestion. The fragments released from GI-digested CSNAIPEL were still predicted to be poorly absorbed. Although PIY released from GI-digested LPIY potentially has high GI absorption, its docking energy score when bound to RBD was weaker (less negative) than that of PW. Thus, comparison with other natural peptides in the literature highlights the promising potential of seed peptides in this study, particularly when considering their GI resistance and hence potential bioavailability upon oral consumption. In particular, pumpkin seed-derived PW is the most promising peptide in our high-GI-absorption peptide dataset, being a putative multi-target antagonist against SARS-CoV-2 M pro , PL pro , and spike glycoproteins. Visualization of protein-peptide docking results by using LigPlot + revealed that molecular interactions with SARS-CoV-2 viral proteins varied among the 56 high-GI-absorption seed peptides (data not shown). Interestingly, 36 of 56 high-GI-absorption peptides analyzed were found to be capable of binding to at least one key binding residue on spike glycoprotein RBD, as well as binding to at least one residue in the catalytic dyad of M pro and in the catalytic triad of PL pro . The 36 peptides include 19 dipeptides and 17 tripeptides ( Table 2 ). The 36 peptides may have the versatility of targeting different phases that are critical to the development of COVID-19 infection. Remarkably, 27 of the 36 multitarget, high-GI-absorption peptides were released from quinoa seed storage proteins. Thus the quinoa seed could be a very promising source of anti-COVID-19 peptides upon oral consumption and GI digestion. The binding poses of the five top-score seed peptides on the surfaces of the M pro , PL pro and spike glycoprotein RBD are shown in Fig. 1 . Detailed interactions between the three viral proteins and the five topscore seed peptides, as revealed by LigPlot + analysis, are presented in Table 3 and Fig. S1 -S3. Seed peptides PW, PSF, IVF and PIL could interact with both catalytic-dyad residues His41 and Cys145 on M pro either via hydrophobic interactions, hydrogen bond, or both. The same ap-plies to the PIY fragment released from reference peptide LPIY following GI digestion (data not shown). The seed peptides were predicted to also interact with the M pro active site residues Met49, Leu141, Asn142, His163, His164, Met165, Glu166, Asp187 and Gln189 ( Narkhede et al., 2020 ) ; such interactions were also observed between SARS-CoV-2 M pro and its inhibitor N3 ( Jin et al., 2020 ) . LigPlot + analysis on PL pro revealed that PW, PSF, PGF, IVF and PY could interact with catalytic-triad residues Cys111, His272 and Asp286. Other than that, the main interactions of peptides with PL pro involved residues Trp106, Asn109, Leu162 and Gly271 ( Table 3 ). The four aforementioned residues are part of the substrate-binding pocket of SARS-CoV-2 PL pro ( Ortega et al., 2020 ) . The same interactions with these four residues were also be detected between PL pro and reference peptide p28 (data not shown). Seed peptides PSF, PGF and IVF could each form at least one interaction with Gln493, a key binding residue on RBD, besides binding to Tyr453, Gly496 and Tyr505 that were previously reported to participate in binding to hACE2 ( Wang et al., 2020a ) . PW and APY were also found to interact with Tyr453, Gly496, Asn501 and Tyr505 which participates in binding to hACE2 ( Wang et al., 2020a ) , despite no interactions with any key binding residues on RBD. In short, all the five top-score peptides analyzed in this study were predicted to bind to the catalytic site residues on M pro and PL pro , as well as to the hACE2-binding residues on the spike glycoprotein RBD. The top-score peptides corresponding to M pro , PL pro and spike glycoprotein RBD, namely PW, PSF, IVF, PGF, PIL, PY and APY, comprise mainly hydrophobic residues (Ala, Ile, Leu, Phe, Pro, Trp and Val). At least 50% of the residues in the aforementioned peptides are hydrophobic. Hydrophobic residues also make up of 25-75% of the reference peptides LPIY, QRPR, CSNAIPEL, p18 and p28. Notably, Pro is the most common residue among the top-score seed peptides and the reference peptides investigated in this study. Hydrophobic interactions contributed about 80% of the total number of interaction between the seed peptides and the viral proteins ( Table 3 ) , suggesting their significance in forming stable protein-peptide complexes. We also observed that hydrophobic interactions account for more than 60% of the total number of protein-peptide interaction between all the reference peptides and the three viral proteins (data not shown). PW, the multi-target, high-GI-absorption seed peptide with the lowest (most negative) docking energy score as computed by HPEPDOCK, was also compared with phytochemicals from edible and/or medicinal plants. For fair comparison, we focused on only phytochemicals that were predicted to show high GI absorption. In a preliminary step, we screened 28, 25 and 26 phytochemicals previously reported as M pro , PL pro and spike glycoprotein RBD inhibitors, respectively, on the SwissADME server for high GI absorption. Consequently, 15 high-GIabsorption phytochemicals were chosen for M pro docking analysis, 17 N = 3) . Different lowercase letters on the left of the bars indicate statistically significant difference ( P < 0.05) between mean values. Black bars represent PW. The numbers on the vertical axis in (A) represent: 1, Dihomo--linolenic acid; 2, p-Coumaroyltyramine; 3, Lignan; 4, Sugiol; 5, PW; 6, Ursodeoxycholic acid; 7, N-cis-feruloyltyramine; 8, Deserpidine; 9, Quercetin; 10, Neoandrographolide; 11, Rhein; 12, Tanshinone IIA; 13, Berberine; 14, Kaempferol; 15, Cryptotanshinone; 16, Absinthin. The numbers on the vertical axis in (B) represent: 1, N-cis-feruloyltyramine; 2, Magnolol; 3, 14-deoxyandrographolide; 4, Piceatannol; 5, Kaempferol; 6, Chrysin; 7, PW; 8, 14-deoxy-11,12-didehydroandrographolide; 9, Quercetin; 10, p-Coumaroyltyramine; 11, Moupinamide 12, Andrographolide; 13, Neoandrographolide; 14, Phaitanthrin D; 15, GRL0617; 16, Cryptotanshinone; 17, Tanshinone IIA; 18, Withanolide A. The numbers on the vertical axis in (C) represent: 1, Geraniol; 2, Anethole; 3, Cinnamaldehyde; 4, Pulegone; 5, l -4-terpineol; 6, Cinnamyl acetate; 7, Thymol; 8, Carvacrol; 9, Curcumin; 10, Chrysophanol; 11, Neoandrographolide; 12, Quercetin; 13, Apigenin; 14, Kaempferol; 15, PW; 16, Dihydrotanshinone I; 17, Fisetin. for PL pro and 16 for spike glycoprotein RBD (Tables S4-S6 ). The binding affinity of PW ( − 6.9 kcal/mol) when docked on M pro , was 28% stronger than dihomo--linolenic acid, but 26% weaker than absinthin ( Fig. 2 A) . When docked on PL pro , the binding affinity of PW ( − 5.8 kcal/mol) was up to 16% stronger than magnolol and N-cis-feruloyltyramine, but 29% weaker than 11 phytochemicals ( Fig. 2 B) . Thus, when compared with other phytochemicals, PW has moderate binding affinities to M pro and PL pro . Evaluation of binding to the spike protein RBD ( Fig. 2 C) highlights the superior binding affinity of PW over most of the reference phytochemicals analyzed. Specifically, the binding affinity of PW towards spike protein RBD was comparable to those of kaempferol and apigenin, surpassing 12 other reference phytochemicals. Taken together, seed peptide PW is comparable or superior to some potentially high-GI-absorption phytochemicals previously reported to bind to the M pro , PL pro and spike protein RBD of SARS-CoV-2. Our findings thus suggest that like the natural products derived from edible and medicinal plants, seed storage peptides could also be a natural resource of SARS-CoV-2 inhibitors for combating COVID-19. Notwithstanding, bioinformatic studies have their limitations. For example, such a theoretical or computational approach cannot reveal the relative quantities of the different types of peptides that could be produced upon in vivo GI digestion. Thus, future in vitro and in vivo evaluations are still necessary to validate our findings, focusing on the seven most promising peptides highlighted in this study. By combining the methods of in silico GI digestion, GI absorption screening, and molecular docking, this study discovered 36 high-GIabsorption peptides with putative anti-COVID-19 potential from the seed storage proteins of five edible seeds. These multi-target peptides could bind to the key binding or catalytic site residues of SARS-CoV-2 spike glycoproteins, M pro and PL pro . Such peptides could potentially suppress viral attachment and replication concurrently. Quinoa was the most promising among the five seeds, releasing the largest number of multi-target peptides. Thus, its potential as a medicinal food deserves greater attention when designing nutrition approach for curbing COVID-19. Overall, the outcome of this study lays a foundation for future search of multi-target anti-COVID-19 peptides from seed storage proteins. Antiviral peptides: identification and validation Sensory and structural characterization of umami peptides derived from sunflower seed The protein data bank Enzyme-assisted discovery of antioxidant peptides from edible marine invertebrates: a review 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 Dermaseptin-based antiviral peptides to prevent COVID-19 through in silico molecular docking studies against SARS-Cov-2 spike protein Novel dipeptidyl peptidase-IV and angiotensin-I-converting enzyme inhibitory peptides released from quinoa protein by in silico proteolysis Peptides as drug candidates: limitations and recent development perspectives Identification of angiotensin converting enzyme and dipeptidyl peptidase-IV inhibitory peptides derived from oilseed proteins using two integrated bioinformatic approaches Chemical and nutritional evaluation of pumpkin ( Cucurbita pepo ) seed proteins Structure of Mpro from SARS-CoV-2 and discovery of its inhibitors PubChem substance and compound databases LigPlot + : multiple ligand-protein interaction diagrams for drug discovery BIOPEP-UWM database of bioactive peptides: current opportunities Insights into SARS-CoV-2 genome, structure, evolution, pathogenesis and therapies: structural genomics approach Recognition of natural products as potential inhibitors of COVID-19 main protease (Mpro): in-silico evidences Class A G protein-coupled receptor antagonist famotidine as a therapeutic alternative against SARS-CoV2: an in silico analysis The crystal structure of papain-like protease of SARS CoV-2 Recent developments in molecular docking technology applied in food science: a review UniProt: a worldwide hub of protein knowledge AutoDock Vina: improving the speed and accuracy of docking with a new scoring function, efficient optimization, and multithreading In vitro chemopreventive properties of peptides released from quinoa ( Chenopodium quinoa Willd.) protein under simulated gastrointestinal digestion LIGPLOT: a program to generate schematic diagrams of protein-ligand interactions Influence of peptide characteristics on their stability, intestinal transport, and in vitro bioavailability: a review Structural and functional basis of SARS-CoV-2 entry by using human ACE2 Novel ACE inhibitory peptides derived from simulated gastrointestinal digestion in vitro of sesame ( Sesamum indicum L.) protein and molecular docking study Ligand-centered assessment of SARS-CoV-2 drug target models in the protein data bank Identification of putative cell-entry-inhibitory peptides against SARS-CoV-2 from edible insects: an in silico study Advances on the antioxidant peptides from edible plant sources Key residues of the receptor binding motif in the spike protein of SARS-CoV-2 that interact with ACE2 and neutralizing antibodies Corona virus versus existence of human on the earth: a computational and biophysical approach HPEPDOCK: a web server for blind peptide-protein docking based on a hierarchical algorithm The authors declare no financial or commercial conflict of interest Supplementary material associated with this article can be found, in the online version, at doi:10.1016/j.phyplu.2020.100016 .