Next Article in Journal
Consumption of Antibiotics in Primary Care Setting before and during COVID-19 Pandemic in Republic of Srpska, Bosnia and Herzegovina
Next Article in Special Issue
Recent Progress in the Discovery and Design of Antimicrobial Peptides Using Traditional Machine Learning and Deep Learning
Previous Article in Journal
In Silico Study for Algerian Essential Oils as Antimicrobial Agents against Multidrug-Resistant Bacteria Isolated from Pus Samples
Previous Article in Special Issue
Overcoming Methicillin-Resistance Staphylococcus aureus (MRSA) Using Antimicrobial Peptides-Silver Nanoparticles
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Antibiotics
Volume 11
Issue 10
10.3390/antibiotics11101318
antibiotics-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Virtual Screening for SARS-CoV-2 Main Protease Inhibitory Peptides from the Putative Hydrolyzed Peptidome of Rice Bran

by
Nathaphat Harnkit
1,†,
Thanakamol Khongsonthi
2,†,
Noprada Masuwan
2,†,
Pornpinit Prasartkul
2,†,
Tipanart Noikaew
3 and
Pramote Chumnanpuen
4,5,*
1
Medicinal Plant Research Institute, Department of Medical Sciences, Ministry of Public Health, Nonthaburi 11000, Thailand
2
Mahidol Wittayanusorn School, 364 Salaya, Phuttamonthon District, Nakhon Prathom 73170, Thailand
3
Department of Biology and Health Science, Mahidol Wittayanusorn School, 364 Salaya, Phuttamonthon District, Nakhon Prathom 73170, Thailand
4
Omics Center for Agriculture, Bioresources, Food and Health, Kasetsart University (OmiKU), Bangkok 10900, Thailand
5
Department of Zoology, Faculty of Science, Kasetsart University, Bangkok 10900, Thailand
*
Author to whom correspondence should be addressed.
†
These authors contributed equally to this work.
Antibiotics 2022, 11(10), 1318; https://doi.org/10.3390/antibiotics11101318
Submission received: 31 August 2022 / Revised: 21 September 2022 / Accepted: 26 September 2022 / Published: 27 September 2022
(This article belongs to the Special Issue Antimicrobial Peptides - Discovery, Structure, Function, and Application)
Browse Figures
Versions Notes

Abstract

:
The Coronavirus Disease 2019 (COVID-19) caused by the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) has led to the loss of life and has affected the life quality, economy, and lifestyle. The SARS-CoV-2 main protease (Mpro), which hydrolyzes the polyprotein, is an interesting antiviral target to inhibit the spreading mechanism of COVID-19. Through predictive digestion, the peptidomes of the four major proteins in rice bran, albumin, glutelin, globulin, and prolamin, with three protease enzymes (pepsin, trypsin, and chymotrypsin), the putative hydrolyzed peptidome was established and used as the input dataset. Then, the prediction of the antiviral peptides (AVPs) was performed by online bioinformatics tools, i.e., AVPpred, Meta-iAVP, AMPfun, and ENNAVIA programs. The amino acid composition and cytotoxicity of candidate AVPs were analyzed by COPid and ToxinPred, respectively. The ten top-ranked antiviral peptides were selected and docked to the SARS-CoV-2 main protease using GalaxyPepDock. Only the top docking scored candidate (AVP4) was further analyzed by molecular dynamics simulation for one nanosecond. According to the bioinformatic analysis results, the candidate SARS-CoV-2 main protease inhibitory peptides were 7–33 amino acid residues and formed hydrogen bonds at Thr22–24, Glu154, and Thr178 in domain 2 with short bonding distances. In addition, these top-ten candidate bioactive peptides contain hydrophilic amino acid residues and have a positive net charge. We hope that this study will provide a potential starting point for peptide-based therapeutic agents against COVID-19.

1. Introduction

The COVID-19 pandemic leads to a high mortality rate, affects the work system, food system, and economy, and causes changes in human lifestyles [1]. Moreover, infection with Severe Acute Respiratory Syndrome Coronavirus-2 (SARS-CoV-2) provokes a mutation of COVID-19 or Coronavirus Disease 2019 [1,2]. The virus passes through the body into the respiratory tract, where angiotensin-converting enzyme 2 (ACE2), which is present on the cell surface, binds to the viral spike protein and enables the virus to enter the cell. Since ACE2 can capture the virus on it spike protein, the virus can enter the host cell. The main protease (Mpro) of SARS-CoV-2 plays a pivotal role in mediating the replication and transcription of viral genes. Mpro hydrolyzes the polyprotein at at least eleven conserved sites and begins by cleaving the pp1a and pp1b of Mpro. For this reason, finding Mpro inhibitors can help reduce the spread of COVID-19. As a result, proteases are attractive antiviral targets [1]. One of the consequences of the COVID-19 epidemic is the shortage of resources and materials for drug research and development. Natural extracts represent an attractive option due to the ready availability of materials.
Since plant seeds naturally accumulate proteins for long-term storage, they are suitable for the production of stable and large amounts of antimicrobial peptides. In particular, rice seed and rice bran are considered good protein sources [3]. In rice (Oryza sativa), 60–80% of all storage proteins are glutelins, 20–30% are prolamines, and about 5% are globulins [4,5,6]. Interestingly, there are several reports on the broad biological activities of proteins and peptides derived from rice bran, i.e., anti-cancer, anti-inflammatory, anti-diabetic, therapy for chronic diseases, and anti-COVID-19 effects [7,8,9,10,11,12].
Bioinformatics is the application of tools of computation and analysis to capture, interpret, and evaluate data in molecular biology. A great benefit of using bioinformatics is that it can save time and resources. Moreover, it provides a base for creating three-dimensional models of complex molecules [13]. Furthermore, it is safer in terms of cost and time to use bioinformatics in research than in the laboratory. Recently, docking studies of some rice bran peptides with the integrin αIIbβ3 and ACE 2 receptors demonstrated that the rice peptide has potential against SARS-CoV-2 infection [12]. However, the antiviral peptide screening of the putative hydrolyzed rice bran peptidome and molecular docking of the SARS-CoV-2 main protease have not been studied.
Therefore, this study aimed to develop a bioinformatic workflow for the virtual screening of antiviral peptides and use in silico protein–peptide docking studies to select the inhibitory peptide for the SARS-CoV-2 main protease.

2. Results and Discussion

2.1. Rice Bran Putative Antiviral Peptide Screening using a Computational Method

All of the peptide sequences from this putative peptidome of rice bran (Oryza sativa) are shown in the supplementary file (Table S1) as the predicted cut results of specific protease enzyme and original protein sources. The peptides were predicted to be antiviral peptides (AVPs) if the cut-off criteria at the support vector machine (SVM) probability exceeded 50 in AVPpred (either in Model 3: Composition or Model 4: Physicochemical properties) and the probability was greater than 0.5 in Meta-iAVP and AMPfun (only antiviral activity). As a result, AVPpred predicted 28 AVPs; Meta-iAVP predicted 73 AVPs; AMPfun predicted 244 AVPs; and there are 71 AVPs in common from at least two prediction servers from 292 input peptide sequences (Figure 1). In terms of the peptides’ properties (Figure 1), the majority of putative AVPs were 5–14 amino acids long (76%), hydrophilic (63%), and cationic (56%) (Figure 2).
Focusing on the length distribution of candidate AVPs and upon further analysis of the peptide length determining antiviral activity, it was observed that 58 (76%) and 14 (19%) were 5–14 and 15–24 amino acids in length, respectively. In general, AVPs are usually short (8–40 amino acids in length) and consist of cationic amino acid residues [14,15,16,17,18,19]. The most optimal peptide length for AMPs, especially AVPs, is 10–40 amino acid residues and hence, it is of great value to optimize the AVP length [20,21].
A comparison of the amino acid composition of AVPs is shown in Figure 3. Compared to the non-AVPs, the preferential amino acid residues with higher percentages were Lys, Leu, Ile, Val, Trp, and Gln. This is consistent with what have been reported before for AVPs, which work as antimicrobial peptides and possess cationic and amphipathic characteristics with positive net charges [14,22,23]. In particular, this specific basic residue (Lys) is the preferential amino acid for antiviral activity and is abundantly found as a component of therapeutic peptides with enhanced electrostatic properties [14]. This also facilitates the interaction and insertion of peptides into the anionic cell walls and phospholipid membranes of microorganisms [22,24]. Therefore, Lys (not Arg or His) plays an important role in providing cationic properties to enhance the antiviral activities of AVPs [14,25]. However, the reason for AVPs bias toward lysine is not clear: either their potential to inhibit the viral enzymes or the inhibition of the viral-entering process [14,26]. On the other hand, the medium-sized hydrophobic residues (Leu and Ile) were also found to be the most abundant residues in the AVPs and play important role in the amphipathic characteristics of antimicrobial peptides [14,27].

2.2. Prediction Scores and Physicochemical Properties of the Ten Top-ranked AVPs

The selection of the ten top-ranked AVPs was based on their predictive probability scores or SVM scores as mentioned in the previous section. The predictive scores and amino acid sequence details of the ten top-ranked antiviral peptides predicted from the putative rice bran peptidome by using four bioinformatics programs are shown in Table 1. Since the secondary structure is one of the most important peptide sequence features for predicting AVPs, each peptide secondary structure was predicted by PEP-FOLD3. These peptides were structurally classified as two main groups; random coiled (AVP1, 2, 4, 9, and 10) and helix-containing loops (AVP3, 5, 6, 7, and 8). Consistent with previous studies, random coils and α-helices were reported to be the top two dominant secondary structures of AVPs compared to the β-sheet structure [14,28].
All of the predictive scores that were greater than the cut-off criteria (the support vector machine (SVM) probability exceeded 50 in AVPpred and the probability was higher than 0.5 in Meta-iAVP and AMPfun) were highlighted in red. Even though there are four SVM-based models available in the AVPpred program based on the peptide sequence features, only the amino acid compositions (Model 3) and physiochemical property (Model 4) models can provide predictive scores. Since the antiviral peptide motifs (Model 1) and sequence alignment (Model 2) could give only two classifications (AVPs and Non-AVPs), we only took the predictive SVM scores from Model 3 (M3) and Model 4 (M4) for this particular AVP candidate ranking. In the ENNAVIA program, there are two main sets of the neural network prediction models. The first two models were used for antiviral property classification (ENNAVIA-A and B, based on the antiviral vs. non-antiviral datasets and antiviral vs. random datasets, respectively) while the other two models were used for anti-coronavirus property prediction (ENNAVIA-C and D, based on the anti-coronavirus vs. non-antiviral datasets and anti-coronavirus vs. random datasets, respectively) [29].
The calculated scores of the physicochemical properties of all of the selected AVP candidates (AVP1 to AVP10) obtained from the ToxinPred web server are shown in Table 2. The largest peptide was AVP4 with a molecular weight of 3916.01 g/mole while AVP10 was the smallest peptide with a molecular weight of 792.01 g/mole. All top-ranked putative AVPs were classified as amphipathic characteristics with steric hindrance and sidebulk scores higher than 0.5.

2.3. Protein–peptide Docking Simulations

Bioinformatics prediction of the SARS-CoV-2 main protease inhibitory peptides is challenging due to lack of the computational predictors available. The molecular docking of these ten top-ranked AVPs to the crystal structure of the SARS-CoV-2 main protease demonstrated that AVP1 to AVP10 binds near the active site of the SARS-CoV-2 main protease structure (Figure 4). Based on the visualization of the protein–peptide docking simulation results, all hydrogen bonds observed from the molecular docking of ten peptides (AVP1 to AVP10) to the crystal structure of the SARS-CoV-2 main protease are listed in Table 3. In general, the proper hydrogen bond acceptor–donor pair is within the correct distance (2.7 to 3.3 Å) [30], and most hydrogen bonds between peptide molecules and protein macromolecules are relatively strong, starting with a 3 Å bond length [31]. According to our molecular docking results, the distance of the hydrogen bonds between selected AVPs and the SARS-CoV-2 main protease structure ranged from 1.8 to 3.8 Å. In particular, hydrogen bonds with donor–acceptor distances of 2.2–2.5 Å are considered to be “strong, mostly covalent” while those of 2.5–3.2 Å and 3.2–4.0 Å are considered to be “moderate, mostly electrostatic” and “weak, electrostatic”, respectively [32]. A comparison of the interaction of the SARS-CoV-2 main protease structure indicated that all selected AVPs showed similar binding positions to the SARS-CoV-2 main protease structure (Figure 4). All ten selected AVPs could interact with the amino acid residues near the active site of the SARS-CoV-2 main protease by hydrogen bonding (Figure 5 and Table 3). The top-ranked AVP candidate (AVP1) with the highest predictive anti-coronavirus activity score formed five hydrogen bonds (with Gln180, Thr178, Glu154, Cys42, and Gly131). Notably, the largest AVP candidate (AVP4) interacted with 19 residues by forming 24 hydrogen bonds near the active site area, mainly the Gln22 and Arg33 residues (Figure 5 and Table 3).
The docking positions of all selected AVPs closer to the active site of the SARS-CoV-2 main protease structure were very similar to the binding pocket of the known main protease inhibitory marine polyketides [33] and antiviral drugs [34]. The inhibitory peptides tended to form hydrogen bonds with the Thr22–24, Glu154, and Thr178 residues of the SARS-CoV-2 main protease. These hydrogen bonds were found in very short distances (1.8–2.5 Å) and are considered to be “strong, mostly covalent” interactions, indicating a high binding affinity.
According to the binding affinity and binding energy analysis by the PROGIDY and PIMA web servers (Table 4), AVP4 showed the strongest binding to the binding groove of the SARS-CoV-2 main protease enzyme with a molecular docking score of −363.04 kJ/mol and a binding affinity (ΔG) of −14.5 (kcal/mol). This high docking score was better than the docking scores of SARS-CoV-2 Mpro for several drugs, i.e., noscapine, chloroquine, ribavirin, and favipiravir (−292.42, −269.71, −214.17, and −153.91 kJ/mol, respectively) [34]. The Mpro–AVP4 interface interactions involved in affinity and binding energy distribution were obtained by high van der Waals energy, −325.48 kcal/mol (hydrophobic interactions) and hydrogen bond energy −37.56 kcal/mol, which are considered to be the most significant calculations to assess the binding stability. Based on the molecular docking results of this study, it was proven that the selected AVPs have potential as candidates for the SARS-CoV-2 main protease inhibitor in controlling the COVID-19 disease.

2.4. Simulation of the Molecular Dynamics of the Mpro–AVP4 Complex

Based on the highest ranked binding affinity and docking scores, the AVP4 binding to the Mpro enzyme was assessed by molecular dynamics simulation for one nanosecond to examine the conformational stability and fluctuation analysis of the complex. To determine the structural activity of the macromolecule, the radious of gyration (Rg) of resulting trajectories was calculated. The Rg level varied according to the folding state of the protein–peptide complex and fluctuated between 17.7 and 18.0 Å, and minimal fluctuations showed the stability of the Mpro macromolecule while binding to AVP4 (Figure 6A). After that, the fraction of the native contacts Q of trajectories was analyzed to define the transition states for protein folding (conformational changes) with free energy. The result depicted that the native contacts were favored by the coarse-grained theoretical models. Q values were determined to be above 96% and depicted the conformational dynamics of the Mpro enzyme along with the energetics of bound ligands (AVP4) (Figure 6B). Moreover, the hydrogen bonds were analyzed to reveal the dynamic equilibration of the complex system with a high number of hydrogen bonds and demonstrated the stable binding of AVP4 with the target Mpro enzyme (Figure 6C). These analyses demonstrated the prolonged and robust binding of AVP4 with target Mpro of coronavirus and the involvement of the potential binding energies with the correlation of the MD calculations and the stability of the Mpro enzyme−AVP4 complex.

2.5. In Silico Toxicity Analysis of the Selected AVPs

Toxic side effects are usually considered to ensure high specificity and low cross-reactivity when designing effective, safe, and theoretically infallible therapeutic molecules. In particular, the computational screening of non-toxic peptide approaches is required to improve the selectivity of therapeutic peptides with less cost and time [35]. The online bioinformatic tool, ToxinPred, was used to predict and estimate the toxicity of the putative AVPs to the host cell. There are eight predictive models available for peptide toxicity analysis; SVM (Swiss-Prot) based, SVM (Swiss-Prot) + Motif based, SVM (TrEMBL) based, SVM (TrEMBL) + Motif based, Monopeptide (Swiss-Prot), Monopeptide (TrEMBL), Di-peptide (Swiss-Prot), and Dipeptide (TrEMBL). In order to comprehensively profile the potential toxic side effects to the host cells, all eight models were used for this particular analysis. All of the predictive scores can be observed in Table 5. The negative scores indicate non-toxic classified results while the positive scores indicate the possible toxicity of the analyzed peptides (highlighted in red). Three selected AVPs (AVP1, 7, and 8) were predicted to be potentially toxic in some predictive models (only from the Quantitative Matrix (QM) method). The rest of the candidates (AVP2 – 6, 9, and 10) had ToxinPred scores lower than the others, suggesting less possible toxicity to the host cells.
Based on the fact that the ToxinPred server can predict the toxicity or non-toxicity of peptides with higher accuracy in the SVM method (93.92%) compared to the QM method (88.00%) [36], the SVM model seems to be more reliable. Since there are some other factors involved in the real biological systems (i.e., secondary structure, in vivo instability, the bioactive activity at a specific pH, temperature, and tonicity, etc.), laboratory experiments would still be required to confirm the real toxicity of all candidate peptides.
Even though the therapeutic peptides have a broad spectrum of targets and low toxicity in general, there are also some limitations and challenges of therapeutic peptide development to be considered for further applications [37]. For example, most of them are limited in oral bioavailability with a short half-life and rapid clearance, and some of them contain immunogenic sequences with some potential to cause an allergenic effect in some patients [37,38]. To overcome these challenging problems, several optimized solutions have been proposed by researchers. A previous study suggested that the side chain of non-polar aromatic amino acids (Trp and Phe) can promote peptide structural stability by restricting their conformation through hydrogen bond formation [39]. In addition, it has been reported that non-polar aliphatic amino acids (Ala, Ile, Leu, and Val) are responsible for the thermal stability of proteins and peptides [40]. Taking this constraint as a guideline, we could selectively screen and/or redesign peptides to improve not only the stability of the therapeutic peptides but also other physicochemical properties involving specific cellular and molecular functions. Beside ToxinPred, there are also some web servers for allergenic peptide prediction, i.e., AllerTOP [36] and AllerFP [41], to avoid possible side effects of therapeutic peptides on the host cells [42,43].

3. Materials and Methods

According to the pipeline illustrated in Figure 1, we proposed the bioinformatic virtual screening workflow with in silico validation by protein–peptide molecular docking. The workflow begins with the predicted digestion peptidomes of the four major proteins in rice bran; albumin, glutelin, globulin, prolamin, with three protease enzymes (pepsin, trypsin, and chymotrypsin). The reason for this putative peptidome of rice bran hydrolyzed by these three proteases was that there are several relevant research reports on the high potential bioactivities of enzymes prepared rice bran protein hydrolysates using digestive enzymes [44]. The functional peptides from rice bran prepared by pepsin, trypsin, and chymotrypsin showed high efficiency and antioxidant bioactivity [45], as well as ACE-inhibitory [46], antimicrobial [47], and tyrosinase-inhibition activity [48]. Having specific cut sites on polypeptide sequences, pepsin, trypsin, and chymotrypsin have been beneficially used for amphipathic and/or cationic therapeutic peptide screening from food-derived peptides [44]. In this study, the putative hydrolyzed peptidome was established and used as the input datasets of 292 peptides for the selection of the peptides with the ten top-ranked predicted scores of unique AVPs (also predicted as non-toxic to the host cells).

3.1. Preparation of the Rice Bran Putative Hydrolyzed Peptidome

Since several predictive programs (PEP-FOLD 3.5, ENNAVIA, and GalaxyPepDock) require a proper length for their machine learning-based analyses, only those peptides that consisted of at least 6 amino acid residues were further screened. After the removal of short peptides (fewer than 5 amino acids), 292 rice bran peptide sequences were obtained from the predicted cut sites of the four major proteins of rice bran (Oryza sativa: Uniprot taxonomic ID = 4530); albumin, glutelin, globulin, prolamin (from the National Center for Biotechnology Information: NCBI with specific accession numbers; Q01882, Q40689, O65042, and Q40714, respectively). The in silico pepsin, trypsin, and chymotrypsin digestions were performed by the cleaver R package (version: 1.34.1; [49]). There were 3 groups of peptide length distributions, i.e., 5–14 amino acid residues (233 sequences, 80%), 15–24 amino acid residues (42 sequences, 14%), and 25–34 amino acid residues (17 sequences, 6%) (Figure 1).

3.2. The Bioinformatic Prediction and Screening of Antiviral Peptides

The peptide sequences were arranged in FASTA format and were used as input to predict the antiviral properties using 3 online machine learning-based prediction programs, i.e., AVPpred (http://crdd.osdd.net/servers/avppred/submit.php, accessed on 10 January 2022) [50], Meta-iAVP (http://codes.bio/meta-iavp/, accessed on 10 January 2022) [22], and AMPfun (http://fdblab.csie.ncu.edu.tw/AMPfun/run.htm, accessed on 10 January 2022l) [51]. Since there are several bioactivity models to choose from in AMPfun, only the antiviral classification model was considered in this particular study. Venny 2.1.0. (https://bioinfogp.cnb.csic.es/tools/venny, accessed on 12 January 2022) [52] was used to generate Venn diagrams to visualize the significant positive AVP candidates that were classified as AVPs at least by 2 programs from 3 online bioinformatic prediction tools. To obtain the probability score of anti-coronavirus activity, all significant AVP candidates were analyzed with the neural network peptide antiviral and anti-coronavirus activity predictor ENNAVIA (https://research.timmons.eu/ennavia, accessed on 13 January 2022) [29]. The average scores from 4 predictive models of ENNAVIA (2 antiviral and 2 anti-coronavirus predictive models) together with 3 other AVP predictors were considered to rank all of the AVP candidates. Only the ten top-ranked AVPs were selected for further molecular docking performances regarding the SARS-CoV-2 main protease inhibition. To predict the secondary structure of the selected peptides, PEP-FOLD3.0 (https://bioserv.rpbs.univ-paris-diderot.fr/services/PEP-FOLD3, accessed on 15 January 2022) was used to simulate the feasible molecular structure of the ten top-ranked antiviral peptides. Finally, all 8 predictive models based on SVM methods and the Quantitative Matrix (QM) method in ToxinPred (https://webs.iiitd.edu.in/raghava/toxinpred/protein.php, accessed on 15 January 2022) were used to predict whether the peptides were cytotoxic to the host cells. All calculated scores of the physicochemical properties of the selected AVP candidates were also obtained by the ToxinPred web server in the batch submission option.

3.3. The Protein–peptide Molecular Docking Simulation

The crystal structure of the SARS-CoV-2 main protease in the apo state with PDB ID code 7C2Q [53] was accessed from the Protein Data Bank (PDB) (http://www.rcsb.org, accessed on 10 January 2022).
The amino acid sequences of each selected AVP candidate (AVP1 to AVP10) were docked to the SARS-CoV-2 main protease enzyme using the GalaxyPepDock (http://galaxy.seoklab.org/pepdock, accessed on 21 January 2022)) [53]. The docking results of the best model and hydrogen bond finding were visualized by the UCSF Chimera program [54]. The FindHBond tool (in the menu under Tools and Structure Analysis) was used to identify and analyze the H-bonding patterns between each AVP candidate and SARS-CoV-2 main protease. To investigate the protein–peptide interface interactions involved in affinity and binding energy, PIMA [55], available at http://caps.ncbs.res.in/pima (accessed on 24 January 2022), and PRODIGY [56] (https://wenmr.science.uu.nl/prodigy, accessed on 24 January 2022) web servers were employed.

3.4. Molecular Dynamics Simulation Analysis of AVP with the Mpro Enzyme

The molecular dynamics of the Mpro–AVP docked complexes were determined using the Ligand and Receptor Molecular Dynamics (LARMD) web server (http://chemyang.ccnu.edu.cn/ccb/server/LARMD/, accessed on 24 January 2022) [57]. The steered molecular dynamics simulation (Str_mod) was chosen for one-nanosecond molecular dynamics simulation in an explicit water model.

4. Conclusions

In conclusion, all 10 selected AVPs named AVP1 to AVP10 from our proposed bioinformatic virtual screening workflow were quite diverse in length (7–33 amino acid residues). These top-ten AVP candidates contain hydrophilic amino acid residues and have a positive net charge. The molecular docking performances infer that all AVP candidates had significant hydrogen bonding with the SARS-CoV-2 main protease (PDB ID = 7C2Q) active site at Thr22–24, Glu154, and Thr178 in domain 2 with short bonding distances. AVP4 was the best candidate for SARS-CoV-2 Mpro inhibition with the highest affinity and binding energy among all selected AVPs. Further in vitro and in vivo studies must be conducted to authenticate the anti-COVID-19 potential of these AVP candidates.

Supplementary Materials

The following are available online at https://www.mdpi.com/article/10.3390/antibiotics11101318/s1, Table S1: Total peptide sequences of the putative peptidome of rice bran (Oryza sativa) hydrolyzed by pepsin, trypsin, and chymotrypsin.

Author Contributions

Conceptualization, P.C.; methodology, T.K., N.M., N.H. and P.P.; software, P.C.; validation, T.K., N.M., P.P., N.H. and T.N.; formal analysis, T.K., N.M., P.P. and N.H.; investigation, P.C. and T.N.; resources, P.C.; data curation, P.C.; writing—original draft preparation, T.K., N.M., P.P., N.H. and P.C.; writing—review and editing, P.C.; visualization, T.K., N.M., P.P., N.H. and P.C.; supervision, P.C.; project administration, P.C. and T.N.; funding acquisition, P.C. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the Department of Zoology, Faculty of Science, Kasetsart University, and Mahidol Wittayanusorn School, Thailand.

Acknowledgments

The authors would like to thank the Mahidol Wittayanusorn School and Department of Zoology, Faculty of Science, Kasetsart University, for support for all research facilities and equipment.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Wang, M.-Y.; Zhao, R.; Gao, L.-J.; Gao, X.-F.; Wang, D.-P.; Cao, J.-M. SARS-CoV-2: Structure, biology, and structure-based therapeutics development. Front. Cell. Infect. Microbiol. 2020, 10, 587269. [Google Scholar] [CrossRef] [PubMed]
  2. Mittal, A.; Manjunath, K.; Ranjan, R.K.; Kaushik, S.; Kumar, S.; Verma, V. COVID-19 pandemic: Insights into structure, function, and hACE2 receptor recognition by SARS-CoV-2. PLoS Pathog. 2020, 16, e1008762. [Google Scholar] [CrossRef] [PubMed]
  3. Bundó, M.; Montesinos, L.; Izquierdo, E.; Campo, S.; Mieulet, D.; Guiderdoni, E.; Rossignol, M.; Badosa, E.; Montesinos, E.; Segundo, B.S.; et al. Production of cecropin A antimicrobial peptide in rice seed endosperm. BMC Plant Biol. 2014, 14, 102. [Google Scholar] [CrossRef] [PubMed]
  4. Li, X.; Okita, T.W. Accumulation of Prolamines and Glutelins during Rice Seed Development: A Quantitative Evaluation. Plant Cell Physiol. 1993, 34, 385–390. [Google Scholar]
  5. Kawakatsu, T.; Yamamoto, M.P.; Hirose, S.; Yano, M.; Takaiwa, F. Characterization of a new rice glutelin gene GluD-1 expressed in the starchy endosperm. J. Exp. Bot. 2008, 59, 4233–4245. [Google Scholar]
  6. Cagampang, G.B.; Perdon, A.A.; Juliano, B.O. Changes in salt-soluble proteins of rice during grain development. Phytochemistry 1976, 15, 1425–1429. [Google Scholar] [CrossRef]
  7. Wang, M.; Hettiarachchy, N.S.; Qi, M.; Burks, W.; Siebenmorgen, T. Preparation and Functional Properties of Rice Bran Protein. Isolate. J. Agric. Food Chem. 1999, 47, 411–416. [Google Scholar] [CrossRef]
  8. Cicero, A.F.; Derosa, G. Rice bran and its main components: Potential role in the management of coronary risk factors. Curr. Top. Nutraceutical Res. 2005, 3, 29–46. [Google Scholar]
  9. Zaky, A.A.; Abd El-Aty, A.; Ma, A.; Jia, Y. An overview on antioxidant peptides from rice bran proteins: Extraction, identification, and applications. Crit. Rev. Food Sci. Nutr. 2022, 62, 1350–1362. [Google Scholar] [CrossRef]
  10. Phongthai, S.; Rawdkuen, S. Fractionation and characterization of antioxidant peptides from rice bran protein hydrolysates stimulated by in vitro gastrointestinal digestion. Cereal Chem. 2020, 97, 316–325. [Google Scholar] [CrossRef]
  11. Ngamsuk, S.; Hsu, J.-L.; Huang, T.-C.; Suwannaporn, P. Ultrasonication of Milky Stage Rice Milk with Bioactive Peptides from Rice Bran: Its Bioactivities and Absorption. Food Bioprocess Technol. 2020, 13, 462–474. [Google Scholar] [CrossRef]
  12. Gasymov, O.K.; Celik, S.; Agaeva, G.; Akyuz, S.; Kecel-Gunduz, S.; Qocayev, N.M.; Ozel, A.E.; Agaeva, U.; Bakhishova, M.; Aliyev, J.A. Evaluation of anti-cancer and anti-covid-19 properties of cationic pentapeptide Glu-Gln-Arg-Pro-Arg, from rice bran protein and its d-isomer analogs through molecular docking simulations. J. Mol. Graph. Model. 2021, 108, 107999. [Google Scholar] [CrossRef]
  13. BioMedTech. Bioinformatics. Available online: https://www.blockdit.com/posts/5db93b2bb6a0111f3f4a722b (accessed on 5 June 2020).
  14. Chang, K.Y.; Yang, J.-R. Analysis and Prediction of Highly Effective Antiviral Peptides Based on Random Forests. PLoS ONE 2013, 8, e70166. [Google Scholar] [CrossRef] [Green Version]
  15. Skalickova, S.; Heger, Z.; Krejcova, L.; Pekarik, V.; Bastl, K.; Janda, J.; Kostolansky, F.; Vareckova, E.; Zitka, O.; Adam, V.; et al. Perspective of Use of Antiviral Peptides against Influenza Virus. Viruses 2015, 7, 5428–5442. [Google Scholar] [CrossRef] [PubMed]
  16. Ahmed, A.; Siman-Tov, G.; Hall, G.; Bhalla, N.; Narayanan, A. Human Antimicrobial Peptides as Therapeutics for Viral Infections. Viruses 2019, 11, 704. [Google Scholar] [CrossRef] [PubMed]
  17. Nyanguile, O. Peptide Antiviral Strategies as an Alternative to Treat Lower Respiratory Viral Infections. Front. Immunol. 2019, 10, 1366. [Google Scholar] [CrossRef]
  18. Sala, A.; Ardizzoni, A.; Ciociola, T.; Magliani, W.; Conti, S.; Blasi, E.; Cermelli, C. Antiviral Activity of Synthetic Peptides Derived from Physiological Proteins. Intervirology 2018, 61, 166–173. [Google Scholar] [CrossRef]
  19. Vilas Boas, L.C.P.; Campos, M.L.; Berlanda, R.L.A.; de Carvalho Neves, N.; Franco, O.L. Antiviral peptides as promising therapeutic drugs. Cell Mol. Life Sci. 2019, 76, 3525–3542. [Google Scholar] [CrossRef]
  20. Shoombuatong, W.; Schaduangrat, N.; Nantasenamat, C. Unraveling the bioactivity of anticancer peptides as deduced from machine learning. EXCLI J. 2018, 17, 734–752. [Google Scholar]
  21. Chen, W.; Ding, H.; Feng, P.; Lin, H.; Chou, K.-C. iACP: A sequence-based tool for identifying anticancer peptides. Oncotarget 2016, 7, 16895–16909. [Google Scholar] [CrossRef]
  22. Schaduangrat, N.; Nantasenamat, C.; Prachayasittikul, V.; Shoombuatong, W. Meta-iAVP: A sequence-based meta-predictor for improving the prediction of antiviral peptides using effective feature representation. Int. J. Mol. Sci. 2019, 20, 5743. [Google Scholar]
  23. Bulet, P.; Stocklin, R.; Menin, L. Anti-microbial peptides: From invertebrates to vertebrates. Immunol. Rev. 2004, 198, 169–184. [Google Scholar] [CrossRef] [PubMed]
  24. Oren, Z.; Shai, Y. Mode of action of linear amphipathic α-helical antimicrobial peptides. Pept. Sci. 1998, 47, 451–463. [Google Scholar]
  25. Li, S.; Guo, S.; Li, F.; Xiang, J. Functional Diversity of Anti-Lipopolysaccharide Factor Isoforms in Shrimp and Their Characters Related to Antiviral Activity. Mar. Drugs 2015, 13, 2602–2616. [Google Scholar] [CrossRef] [PubMed]
  26. Scott, M.G.; Hancock, R.E.W. Cationic antimicrobial peptides and their multifunctional role in the immune system. Crit. Rev. Immunol. 2000, 20, 24. [Google Scholar] [CrossRef]
  27. Wang, C.-K.; Shih, L.-Y.; Chang, K.Y. Large-Scale Analysis of Antimicrobial Activities in Relation to Amphipathicity and Charge Reveals Novel Characterization of Antimicrobial Peptides. Molecules 2017, 22, 2037. [Google Scholar]
  28. Chowdhury, A.S.; Reehl, S.M.; Kehn-Hall, K.; Bishop, B.; Webb-Robertson, B.-J.M. Better understanding and prediction of antiviral peptides through primary and secondary structure feature importance. Sci. Rep. 2020, 10, 1–8. [Google Scholar]
  29. Timmons, P.B.; Hewage, C.M. ENNAVIA is a novel method which employs neural networks for antiviral and anti-coronavirus activity prediction for therapeutic peptides. Brief. Bioinform. 2021, 22, bbab258. [Google Scholar]
  30. Sterri, S.H.; Fonnum, F. CHAPTER 68—Role of Carboxylesterases in Therapeutic Intervention of Nerve Gas Poisoning. In Handbook of Toxicology of Chemical Warfare Agents; Gupta, R.C., Ed.; Academic Press: San Diego, CA, USA, 2009; pp. 1033–1040. [Google Scholar]
  31. Fakih, T.M. Dermaseptin-based antiviral peptides to prevent COVID-19 through in silico molecular docking studies against SARS-Cov-2 spike protein. Pharm. Sci. Res. 2020, 7, 8. [Google Scholar]
  32. Jeffrey, G.A.; Jeffrey, G.A. An Introduction to Hydrogen Bonding; Oxford University Press: New York, NY, USA, 1997. [Google Scholar]
  33. El-Demerdash, A.; Al-Karmalawy, A.A.; Abdel-Aziz, T.M.; Elhady, S.S.; Darwish, K.M.; Hassan, A.H. Investigating the structure–activity relationship of marine natural polyketides as promising SARS-CoV-2 main protease inhibitors. RSC Adv. 2021, 11, 31339–31363. [Google Scholar] [CrossRef]
  34. Kumar, N.; Sood, D.; van der Spek, P.J.; Sharma, H.S.; Chandra, R. Molecular binding mechanism and pharmacology comparative analysis of noscapine for repurposing against SARS-CoV-2 protease. J. Proteome Res. 2020, 19, 4678–4689. [Google Scholar] [CrossRef] [PubMed]
  35. Robles-Loaiza, A.A.; Pinos-Tamayo, E.A.; Mendes, B.; Ortega-Pila, J.A.; Proaño-Bolaños, C.; Plisson, F.; Teixeira, C.; Gomes, P.; Almeida, J.R. Traditional and Computational Screening of Non-Toxic Peptides and Approaches to Improving Selectivity. Pharmaceuticals 2022, 15, 323. [Google Scholar] [CrossRef] [PubMed]
  36. Dimitrov, I.; Bangov, I.; Flower, D.R.; Doytchinova, I. AllerTOP v. 2—A server for in silico prediction of allergens. J. Mol. Model. 2014, 20, 1–6. [Google Scholar] [CrossRef] [PubMed]
  37. Recio, C.; Maione, F.; Iqbal, A.J.; Mascolo, N.; De Feo, V. The potential therapeutic application of peptides and peptidomimetics in cardiovascular disease. Front. Pharmacol. 2017, 7, 526. [Google Scholar] [CrossRef] [PubMed]
  38. Cook, Q.S.; Burks, A.W. Peptide and Recombinant Allergen Vaccines for Food Allergy. Clin. Rev. Allergy Immunol. 2018, 55, 162–171. [Google Scholar] [CrossRef]
  39. Silva, A.F.; Bastos, E.L.; Torres, M.D.T.; Costa-Da-Silva, A.L.; Ioshino, R.S.; Capurro, M.L.; Alves, F.L.; Miranda, A.; Vieira, R.D.F.F.; Oliveira, V.X. Antiplasmodial activity study of angiotensin II via Ala scan analogs. J. Pept. Sci. 2014, 20, 640–648. [Google Scholar] [CrossRef]
  40. Osorio, D.; Rondón-Villarreal, P.; Torres, R. Peptides: A package for data mining of antimicrobial peptides. Small 2015, 12, 44–444. [Google Scholar] [CrossRef]
  41. IDimitrov, I.; Naneva, L.; Doytchinova, I.; Bangov, I. AllergenFP: Allergenicity prediction by descriptor fingerprints. Bioinformatics 2013, 30, 846–851. [Google Scholar] [CrossRef]
  42. Dimitrov, I.; Doytchinova, I. An alignment-independent platform for allergenicity prediction. In Immunoinformatics; Springer: New York, NY, USA, 2020; pp. 147–153. [Google Scholar]
  43. Bhattacharya, M.; Sharma, A.R.; Patra, P.; Ghosh, P.; Sharma, G.; Patra, B.C.; Ghosh, P.; Sharma, G.; Patra, B.C.; Saha, R.P.; et al. A SARS-CoV-2 vaccine candidate: In-silico cloning and validation. Inform. Med. Unlocked 2020, 20, 100394. [Google Scholar] [CrossRef]
  44. Liu, Y.Q.; Strappe, P.; Shang, W.T.; Zhou, Z.K. Functional peptides derived from rice bran proteins. Crit. Rev. Food Sci. Nutr. 2017, 59, 349–356. [Google Scholar] [CrossRef]
  45. Wattanasiritham, L.; Theerakulkait, C.; Wickramasekara, S.; Maier, C.S.; Stevens, J.F. Isolation and identification of antioxidant peptides from enzymatically hydrolyzed rice bran protein. Food Chem. 2016, 192, 156–162. [Google Scholar] [CrossRef] [PubMed]
  46. Wang JHu, J.; Cui, J.; Bai, X.; Du, Y.; Miyaguchi, Y.; Lin, B. Purification and identification of a ACE inhibitory peptide from oyster proteins hydrolysate and the antihypertensive effect of hydrolysate in spontaneously hypertensive rats. Food Chem. 2008, 111, 302–308. [Google Scholar] [CrossRef] [PubMed]
  47. Taniguchi, M.; Kameda, M.; Namae, T.; Ochiai, A.; Saitoh, E.; Tanaka, T. Identification and characterization of multifunctional cationic peptides derived from peptic hydrolysates of rice bran protein. J. Funct. Foods 2017, 34, 287–296. [Google Scholar] [CrossRef]
  48. Ochiai, A.; Tanaka, S.; Tanaka, T.; Taniguchi, M. Rice Bran Protein as a Potent Source of Antimelanogenic Peptides with Tyrosinase Inhibitory Activity. J. Nat. Prod. 2016, 79, 2545–2551. [Google Scholar] [CrossRef] [PubMed]
  49. Gibb, S. Cleaver: Cleavage of Polypeptide Sequences, R package version 1.34.1; 2022. Available online: https://www.bioconductor.org/packages/release/bioc/html/cleaver.html (accessed on 12 January 2022).
  50. Thakur, N.; Qureshi, A.; Kumar, M. AVPpred: Collection and prediction of highly effective antiviral peptides. Nucleic Acids Res. 2012, 40, W199–W204. [Google Scholar] [CrossRef]
  51. Chung, C.-R.; Kuo, T.-R.; Wu, L.-C.; Lee, T.-Y.; Horng, J.-T. Characterization and identification of antimicrobial peptides with different functional activities. Briefings Bioinform. 2019, 21, 1098–1114. [Google Scholar] [CrossRef]
  52. Oliveros, J.; Venny, C. An Interactive Tool for Comparing Lists with Venn’s Diagrams; BioinfoGP, CNB-CSIC; Spanish National Biotechnology Centre: Madrid, Spain, 2007; Available online: https://bioinfogp.cnb.csic.es/tools/venny/ (accessed on 12 January 2022).
  53. Zhou, X.; Zhong, F.; Lin, C.; Hu, X.; Zhang, Y.; Xiong, B.; Yin, X.; Fu, J.; He, W.; Duan, J.; et al. Structure of SARS-CoV-2 main protease in the apo state. Sci. China Life Sci. 2021, 64, 656–659. [Google Scholar] [CrossRef]
  54. Pettersen, E.F.; Goddard, T.D.; Huang, C.C.; Couch, G.S.; Greenblatt, D.M.; Meng, E.C.; Ferrin, T.E. UCSF Chimera? A visualization system for exploratory research and analysis. J. Comput. Chem. 2004, 25, 1605–1612. [Google Scholar] [CrossRef]
  55. Mathew, O.K.; Sowdhamini, R. PIMA: Protein-Protein interactions in Macromolecular Assembly—A web server for its Analysis and Visualization. Bioinformation 2016, 12, 9–11. [Google Scholar] [CrossRef]
  56. Xue, L.C.; Rodrigues, J.P.; Kastritis, P.L.; Bonvin, A.M.; Vangone, A. PRODIGY: A web server for predicting the binding affinity of protein–protein complexes. Bioinformatics 2016, 32, 3676–3678. [Google Scholar]
  57. Yang, J.-F.; Wang, F.; Chen, Y.-Z.; Hao, G.-F.; Yang, G.-F. LARMD: Integration of bioinformatic resources to profile ligand-driven protein dynamics with a case on the activation of estrogen receptor. Brief. Bioinform. 2020, 21, 2206–2218. [Google Scholar] [CrossRef] [PubMed]
Antibiotics 11 01318 g001 550
Figure 1. The workflow of the bioinformatic virtual screening for antiviral peptide (AVP) candidates and the in silico analysis of the SARS-CoV-2 main protease inhibition using molecular docking simulation.
Figure 1. The workflow of the bioinformatic virtual screening for antiviral peptide (AVP) candidates and the in silico analysis of the SARS-CoV-2 main protease inhibition using molecular docking simulation.
Antibiotics 11 01318 g001
Antibiotics 11 01318 g002 550
Figure 2. Percentage of the peptides’ properties with respect to their (a) length, (b) hydrophobicity, and (c) net charge based on 71 putative AVPs.
Figure 2. Percentage of the peptides’ properties with respect to their (a) length, (b) hydrophobicity, and (c) net charge based on 71 putative AVPs.
Antibiotics 11 01318 g002
Antibiotics 11 01318 g003 550
Figure 3. Compositional analysis represents the preferences between the significant AVPs and non-AVPs.
Figure 3. Compositional analysis represents the preferences between the significant AVPs and non-AVPs.
Antibiotics 11 01318 g003
Antibiotics 11 01318 g004 550
Figure 4. Comparative molecular docking of the ten top-ranked AVPs (AVP1 to AVP10) on the crystal structure of the SARS-CoV-2 main protease in the apo state (PDB ID: 7C2Q) shown with (A) and without (B) the enzyme structure. The structure of the SARS-CoV-2 main protease is shaded in gold, and the peptides are labeled with different colors.
Figure 4. Comparative molecular docking of the ten top-ranked AVPs (AVP1 to AVP10) on the crystal structure of the SARS-CoV-2 main protease in the apo state (PDB ID: 7C2Q) shown with (A) and without (B) the enzyme structure. The structure of the SARS-CoV-2 main protease is shaded in gold, and the peptides are labeled with different colors.
Antibiotics 11 01318 g004
Antibiotics 11 01318 g005 550
Figure 5. Molecular docking of the ten top-ranked AVPs (AVP1 (A), AVP2 (B), AVP3 (C), AVP4 (D), AVP5 (E), AVP6 (F), AVP7 (G), AVP8 (H), AVP9 (I), and AVP10 (J)) to the crystal structure of the SARS-CoV-2 main protease in the apo state (PDB ID: 7C2Q). The structure of the SARS-CoV-2 main protease is shaded in gold, and the peptide sequences are colored as labeled above. The hydrogen bonds are shown as red lines.
Figure 5. Molecular docking of the ten top-ranked AVPs (AVP1 (A), AVP2 (B), AVP3 (C), AVP4 (D), AVP5 (E), AVP6 (F), AVP7 (G), AVP8 (H), AVP9 (I), and AVP10 (J)) to the crystal structure of the SARS-CoV-2 main protease in the apo state (PDB ID: 7C2Q). The structure of the SARS-CoV-2 main protease is shaded in gold, and the peptide sequences are colored as labeled above. The hydrogen bonds are shown as red lines.
Antibiotics 11 01318 g005
Antibiotics 11 01318 g006 550
Figure 6. Molecular dynamics simulation of the Mpro–AVP4 complex; radius of teg gyrus plot (A), fraction of native contact analysis (B), and statistics of the hydrogen bonds (C).
Figure 6. Molecular dynamics simulation of the Mpro–AVP4 complex; radius of teg gyrus plot (A), fraction of native contact analysis (B), and statistics of the hydrogen bonds (C).
Antibiotics 11 01318 g006
Table
Table 1. Predictive probability scores or SVM scores from antiviral peptide prediction tools.
Table 1. Predictive probability scores or SVM scores from antiviral peptide prediction tools.
Peptide IDSecondary StructuresSequences
(from Protein/Cut by)		LengthAVPpredMeta-iAVPAMP
Fun		ENNAVIA
M3*M4*A*B*C*D*
AVP1Antibiotics 11 01318 i001SWCRCSALNHMVGGIY
(albumin/pepsin)		1655.5672.370.510.390.471.000.791.00
AVP2Antibiotics 11 01318 i002HQASSLLRGIKNY
(globulin/pepsin)		1354.8224.410.980.490.970.740.570.23
AVP3Antibiotics 11 01318 i003VVFSALLLIIVSVLAATATMADHHK
(albumin/trypsin)		2552.4025.590.540.591.000.001.000.01
AVP4Antibiotics 11 01318 i004QQHSIVATPFWQPATFQLINNQVMQQQCCQQLR
(prolamin/trypsin)		3337.0664.030.060.650.680.000.990.16
AVP5Antibiotics 11 01318 i005IIFVFALLAIVACNASAR
(prolamin/pepsin/trypsin)		1856.1233.780.990.431.000.000.750.00
AVP6Antibiotics 11 01318 i006PILSLVQMSAVKVNLY
(glutelin/pepsin)		1646.6255.770.320.530.100.000.800.32
AVP7Antibiotics 11 01318 i007GGHGPHWPLPPF
(globulin/pepsin)		1242.0425.901.000.510.740.910.000.08
AVP8Antibiotics 11 01318 i008VIALPAGVAHWCY
(glutelin/pepsin)		1341.6762.940.540.460.500.310.350.28
AVP9Antibiotics 11 01318 i009ATILLLLAAVLFAAAAAASGEDR
(globulin/trypsin)		2347.3962.950.740.430.960.000.400.00
AVP10Antibiotics 11 01318 i010QQVGVVY
(prolamin/chymotrypsin)		750.4643.090.990.660.050.200.200.45
Note: For the AVPpred program; M3* = composition of the SVM-based models (Model 3), M4* = physiochemical property SVM-based models (Model 4). For the ENNAVIA program, the neural network models based on the antiviral vs. non-antiviral datasets (A*), antiviral vs. random datasets (B*), anti-coronavirus vs. non-antiviral datasets (C*), and anti-coronavirus vs. random datasets (D*).
Table
Table 2. Calculated physicochemical property scores of the selected AVP candidates (AVP1 to AVP10) by ToxinPred.
Table 2. Calculated physicochemical property scores of the selected AVP candidates (AVP1 to AVP10) by ToxinPred.
Peptide IDHydrophobicitySteric HindranceSidebulkHydropathicityAmphipathicityHydrophilicityNet HydrogenChargepIMol wt
AVP1−0.010.60.60.340.24−0.710.691.58.381797.33
AVP2−0.220.580.58−0.620.68−0.131.082.510.011486.89
AVP30.170.560.561.520.26−0.730.3617.262621.56
AVP4−0.150.610.61−0.450.5−0.5611.58.43916.01
AVP50.230.620.622.10.14−0.940.3918.61892.59
AVP60.090.620.621.050.31−0.740.5618.941775.43
AVP70.080.430.43−0.690.24−0.720.2517.261298.65
AVP80.250.540.541.320.11−1.180.230.57.061399.86
AVP90.150.580.581.450.16−0.450.35−14.382228.94
AVP100.060.690.690.560.36−0.910.7105.88792.01
Table
Table 3. List of hydrogen bonds observed from the molecular docking of ten peptides (AVP1 to AVP10) to the crystal structure of the SARS-CoV-2 main protease in the apo state (PDB ID: 7C2Q).
Table 3. List of hydrogen bonds observed from the molecular docking of ten peptides (AVP1 to AVP10) to the crystal structure of the SARS-CoV-2 main protease in the apo state (PDB ID: 7C2Q).
PeptidesPeptide
Residues		Protease
Residues		Distance
(Å)		PeptidesPeptide
Residues		Protease
Residues		Distance
(Å)
AVP1Arg4Gln1802.485AVP4Arg33Thr192.285
Ser6Thr1782.124(continue)Arg33Thr242.335
Ala7Glu1541.971 Arg33Gln612.041
His10Cys422.978 Arg33Thr192.018
Val12Gly1312.124AVP5Ile1Ala1812.404
AVP2Ser4Gln1801.976 Asn14Thr242.085
Leu6Glu1541.912AVP6Pro1Gln1771.917
Arg8His1512.036 Ile2Glu1542.013
Arg8His1602.047 Val6Asn1112.428
Ile10Thr242.078 Gln7Thr241.902
Tyr13His392.778 Ala10Thr242.463
AVP3Ala5Thr1782.410 Lys12Gly212.137
Leu6Glu1541.875 Lys12Asn571.939
Ser12Asn1112.330 Tyr16Asn552.536
Ala15Thr242.045 Tyr16Val691.995
Thr19Thr222.017AVP7Phe12Gly1312.526
AVP4Ser4Val2862.047 Phe12Cys1331.924
Thr8Pro1562.031 Phe12His1522.091
Trp11Ser2742.279AVP8Val1Asp2091.898
Gln12Asp2092.021 Val1Thr1571.975
Pro13Arg1232.196 Cys12Thr233.071
Thr15Ala1822.471AVP9Ala1Arp2092.000
Thr15Gly1831.938 Ala1Thr1571.861
Asn20Gln1802.430 Thr2Gly2712.093
Gln22Gln1802.057 Thr2Arp1851.972
Gln22Gln1772.493 Val10Glu1541.894
Gln22Thr1782.523 Ala13Thr242.390
Gln22Asp1752.064 Ser19Asn1112.178
Val23Glu1541.935 Glu21Thr242.242
Gln25Ser1322.372 Arp22Gln612.349
Gln25Glu1542.035 Arg23Gln172.163
Gln26Cys422.026 Arg23Asn642.379
Cys28Glu1543.828AVP10Gln1Ala1792.554
Gln30Thr241.993 Val5Glu1541.868
Gln31Thr222.426 Tyr7Ser1322.501
Gln31Thr232.093 Tyr7Glu1542.111
Table
Table 4. Calculated binding affinity (ΔG), dissociation constant (Kd), and binding energy scores from the molecular docking of ten peptides (AVP1 to AVP10) to the SARS-CoV-2 main protease based on the PROGIDY and PIMA web servers. Molecular docking scores, H-bonds, electrostatic energy, and van der Waals energy (H-bond Ener., Elec. Ener., and VDW. Ener.) are presented in kJ/mol.
Table 4. Calculated binding affinity (ΔG), dissociation constant (Kd), and binding energy scores from the molecular docking of ten peptides (AVP1 to AVP10) to the SARS-CoV-2 main protease based on the PROGIDY and PIMA web servers. Molecular docking scores, H-bonds, electrostatic energy, and van der Waals energy (H-bond Ener., Elec. Ener., and VDW. Ener.) are presented in kJ/mol.
Protein-Peptide ComplexΔG (kcal/mol)Kd (M)
at 25.0 °C		H-Bond Ener. (kJ/mol)Elec. Ener. (kJ/mol)VDW. Ener. (kJ/mol)Molecular Docking Score (kJ/mol)
Mpro–AVP1−11.25.8E−09−23.624.48−206.06−225.20
Mpro–AVP2−10.51.9E−08−31.672.43−196.55−225.79
Mpro–AVP3−11.09.3E−09−27.370.00−153.91−181.28
Mpro–AVP4−14.52.2E−11−37.560.00−325.48−363.04
Mpro–AVP5−9.95.7E−08−3.190.00−160.49−163.68
Mpro–AVP6−11.09.1E−09−42.274.91−187.82−225.19
Mpro–AVP7−8.29.3E−07−23.680.00−133.58−157.26
Mpro–AVP8−10.04.5E−08−6.39−13.14−107.89−127.42
Mpro–AVP9−12.65.9E−10−23.330.00−219.64−242.97
Mpro–AVP10−9.41.3E−07−25.320.00−109.52−134.84
Table
Table 5. Predictive SVM and QM scores from ToxinPred.
Table 5. Predictive SVM and QM scores from ToxinPred.
Peptide IDSVM MethodQuantitative Matrix (QM) Method
Model A*Model B*Model C*Model D*Model E*Model F*Model G*Model H*
AVP1−0.22−0.22−0.28−0.2819.709.060.33−0.91
AVP2−0.94−0.94−1.16−1.16−14.50−12.95−3.82−2.13
AVP3−1.40−1.40−1.77−1.77−27.60−46.75−4.69−7.79
AVP4−1.07−1.07−1.14−1.14−1.90−17.13−1.64−1.91
AVP5−0.68−0.68−0.68−0.68−8.50−13.44−2.30−2.32
AVP6−1.31−1.31−2.03−2.03−15.70−23.15−0.32−1.35
AVP7−0.75−0.75−1.30−1.3012.5022.191.281.85
AVP8−0.19−0.19−0.61−0.617.709.28−1.360.43
AVP9−1.69−1.69−1.38−1.38−42.50−52.14−9.10−10.11
AVP10−1.45−1.45−1.18−1.18−16.00−7.16−0.91−0.47
Note: Moldel A = SVM (Swiss-Prot) based, Model B = SVM (Swiss-Prot) + Motif based, Model C = SVM (TrEMBL) based, Model D = SVM (TrEMBL) + Motif based, Model E = Monopeptide (Swiss-Prot), Model F = Monopeptide (TrEMBL), Model G = Dipeptide (Swiss-Prot), and Model H = Dipeptide (TrEMBL).
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Harnkit, N.; Khongsonthi, T.; Masuwan, N.; Prasartkul, P.; Noikaew, T.; Chumnanpuen, P. Virtual Screening for SARS-CoV-2 Main Protease Inhibitory Peptides from the Putative Hydrolyzed Peptidome of Rice Bran. Antibiotics 2022, 11, 1318. https://doi.org/10.3390/antibiotics11101318

AMA Style

Harnkit N, Khongsonthi T, Masuwan N, Prasartkul P, Noikaew T, Chumnanpuen P. Virtual Screening for SARS-CoV-2 Main Protease Inhibitory Peptides from the Putative Hydrolyzed Peptidome of Rice Bran. Antibiotics. 2022; 11(10):1318. https://doi.org/10.3390/antibiotics11101318

Chicago/Turabian Style

Harnkit, Nathaphat, Thanakamol Khongsonthi, Noprada Masuwan, Pornpinit Prasartkul, Tipanart Noikaew, and Pramote Chumnanpuen. 2022. "Virtual Screening for SARS-CoV-2 Main Protease Inhibitory Peptides from the Putative Hydrolyzed Peptidome of Rice Bran" Antibiotics 11, no. 10: 1318. https://doi.org/10.3390/antibiotics11101318

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop