Next Article in Journal
Antibacterial, Antifungal and Anticancer Activities of Compounds Produced by Newly Isolated Streptomyces Strains from the Szczelina Chochołowska Cave (Tatra Mountains, Poland)
Next Article in Special Issue
First Global Report of Plasmid-Mediated mcr-1 and Extended-Spectrum Beta-Lactamase-Producing Escherichia coli from Sheep in Portugal
Previous Article in Journal
Efficacy of Cathelicidin LL-37 in an MRSA Wound Infection Mouse Model
Previous Article in Special Issue
Antibiotic Resistance: From Pig to Meat
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Antibiotics
Volume 10
Issue 10
10.3390/antibiotics10101211
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

Targeting PirAvp and PirBvp Toxins of Vibrio parahaemolyticus with Oilseed Peptides: An In Silico Approach

by
Joe-Hui Ong
1,
Wey-Lim Wong
2,3,
Fai-Chu Wong
1,3 and
Tsun-Thai Chai
1,3,*
1
Department of Chemical Science, Faculty of Science, Universiti Tunku Abdul Rahman, Jalan Universiti, Bandar Barat, Kampar 31900, Perak, Malaysia
2
Department of Biological Science, Faculty of Science, Universiti Tunku Abdul Rahman, Jalan Universiti, Bandar Barat, Kampar 31900, Perak, Malaysia
3
Center for Agriculture and Food Research, Universiti Tunku Abdul Rahman, Jalan Universiti, Bandar Barat, Kampar 31900, Perak, Malaysia
*
Author to whom correspondence should be addressed.
Antibiotics 2021, 10(10), 1211; https://doi.org/10.3390/antibiotics10101211
Submission received: 5 September 2021 / Revised: 28 September 2021 / Accepted: 28 September 2021 / Published: 5 October 2021
(This article belongs to the Special Issue Antibiotic Resistance in Terrestrial and Aquatic Farmed Animals: An Urgent One Health Issue for Sustainability)
Browse Figures
Versions Notes

Abstract

:
Acute hepatopancreatic necrosis disease (AHPND), caused by PirAvp- and PirBvp-releasing Vibrio parahaemolyticus strains, has resulted in massive mortality in shrimp aquaculture. Excessive use of antibiotics for AHPND management has led to antibiotic resistance, highlighting the urgency to search for alternatives. Using an in silico approach, we aimed to discover PirAvp/PirBvp-binding peptides from oilseed meals as alternatives to antibiotics. To search for peptides that remain intact in the shrimp digestive tract, and therefore would be available for toxin binding, we focused on peptides released from tryptic hydrolysis of 37 major proteins from seeds of hemp, pumpkin, rape, sesame, and sunflower. This yielded 809 peptides. Further screening led to 24 peptides predicted as being non-toxic to shrimp, fish, and humans, with thermal stability and low water solubility. Molecular docking on the 24 peptides revealed six dual-target peptides capable of binding to key regions responsible for complex formation on both PirAvp and PirBvp. The peptides (ISYVVQGMGISGR, LTFVVHGHALMGK, QSLGVPPQLGNACNLDNLDVLQPTETIK, ISTINSQTLPILSQLR, PQFLVGASSILR, and VQVVNHMGQK) are 1139–2977 Da in mass and 10–28 residues in length. Such peptides are potential candidates for the future development of peptide-based anti-AHPND agents which potentially mitigate V. parahaemolyticus pathogenesis by intercepting PirAvp/PirBvp complex formation.

Graphical Abstract

1. Introduction

The shrimp is one of the most valuable aquaculture species consumed by people from many parts of the world. However, intensive shrimp culture and improper management have led to the outbreak of several diseases, such as acute hepatopancreatic necrosis disease (AHPND), hepatopancreatic microsporidiosis, infectious hypodermal and hematopoietic necrosis, and white spot disease [1,2,3]. AHPND, commonly known as early mortality syndrome, has caused huge loss in the whiteleg shrimp (Litopenaeus vannamei) industry. This bacterial shrimp disease has resulted in mass shrimp mortality and low shrimp production [3,4]. Administration of antibiotics in shrimp culture has been most widely used in the past in controlling the AHPND infection caused by Vibrio species [4]. Excessive use of antibiotics in shrimp culture has created antibiotic resistance in the pathogenic Vibrio. For instance, AHPND-causing Vibrio campbellii from China was found to carry multiple antibiotic resistance genes [5]. In Mexico, several Vibrio parahaemolyticus AHPND strains were reported to have the tetB gene coding tetracycline resistance [6]. The spread of antibiotic resistance genes among the pathogens may lead to the failure of AHPND control if the management of AHPND mainly depends on antibiotics. Moreover, the contamination of aquatic ecosystems by antibiotics is of immense concern to environmental and human health [7]. Thus, there is an urgent need to develop other anti-AHPND agents as alternatives to antibiotics.
V. parahaemolyticus is a halophilic, Gram-negative bacterium that harbors the pVA1 plasmid and expresses the Photorhabdus insect-related (Pir) binary toxin PirAvp/PirBvp. These two toxins, the virulence factors of AHPND, are reported to form a binary complex, destructing the hepatopancreatic tissues of shrimp [8,9]. Upon infection, AHPND-causing bacteria colonize the stomach, and subsequently release the toxic PirAvp and PirBvp proteins in the hepatopancreas of shrimp [10]. Since PirAvp and PirBvp are both secreted into the digestive tract of shrimp, the proteins can be easily targeted and neutralized by orally administrated compounds or medicated feed [11]. Medicated feed is a common and convenient mode of drug administration in aquaculture [12]. In shrimp farming, oxytetracycline was used to control Vibrio infection through medicated feed [13]. Several AHPND measures based on the antimicrobial properties of probiotics [14], prebiotics [15], plant extracts, or essential oils [16] have been applied in shrimp feeds to prevent outbreaks of AHPND. However, the modes of action of these antimicrobial agents are complicated, and not specific in eliminating the causative agents.
The use of peptides as biofunctional ingredients in shrimp aquaculture for disease control has been highlighted. For instance, the application of bioactive peptides to partially replace fish meal in the diet of shrimp was found to promote growth performance and strengthen immunity [17]. Oilseed meals, which are good sources of bioactive peptides [18,19,20], have been shown to contribute to improved shrimp health when used as shrimp feed ingredients [21,22]. Although the aforementioned studies [21,22] did not investigate how oilseed meals specifically affect the incidence of AHPND, the incorporation of peptides derived from Bacillus subtilis E20-fermented soybean meal into shrimp feed protected shrimp against AHPND [23]. Hence, the application of oilseed-derived bioactive peptides as toxin-neutralizing agents for the inhibition of AHPND could be a promising anti-AHPND strategy. Research in this area could provide valuable information for the formulation of a novel anti-AHPND shrimp feed fortified with oilseed-derived peptides. With the availability of the crystal structures of PirAvp and PirBvp and knowledge of the structural characteristics of the PirAvp/PirBvp complex, structure-based drug design, as recommended by [11,24], is a feasible and efficient approach to discover potential anti-AHPND agents. Specifically, in silico screening for inhibitors of PirAvp/PirBvp complex formation was proposed to be a promising strategy [24]. Therefore, this study adopted an in silico approach to discover oilseed-derived peptides that can bind to the key regions responsible for PirAvp/PirBvp complex formation. The binding of such peptides to the key regions is anticipated to preclude PirAvp/PirBvp complex formation. Major proteins of five oilseed meals, namely Cannabis sativa (hemp), Cucurbita maxima (pumpkin), Brassica napus (rape), Sesamum indicum (sesame), and Helianthus annuus (sunflower) were investigated as sources of PirAvp- and PirBvp- binding peptides following tryptic digestion. Screening for non-toxicity, low water solubility, and high thermal stability was also conducted to increase the chance of discovering peptides that can be readily incorporated as shrimp feed ingredients.

2. Results and Discussion

2.1. In Silico Tryptic Hydrolysis of Oilseed Meal Proteins

To discover peptides that can potentially block the formation of the PirAvp/PirBvp toxin complex, we focused on peptides resulting from the in silico trypsin digestion of oilseed proteins. The in vitro trypsin digestion of feedstuffs is comparable to that of in vivo whiteleg shrimp digestion in terms of protein digestibility [25]. Thus, in silico tryptic digestion can be used to simulate the in vivo shrimp digestion of oilseed proteins and predict the peptide fragments that are released in the shrimp digestive tract. Peptides resulting from tryptic digestion also represent peptides that can be fed to the shrimp directly and remain undegraded in the shrimp digestive tract. Either way, such peptides could remain intact in the shrimp digestive tract, and would hence be available for neutralizing the toxicity of the pathogenic V. parahaemolyticus strains. In silico trypsin hydrolysis released 809 peptide fragments from 37 oilseed meal proteins (Table 1). Rape and hemp seed proteins released 272 and 247 peptide fragments, respectively, accounting for about a 4.5-fold greater than those from the sunflower proteins (Table 1). Regardless of oilseed type, in silico trypsin hydrolysis generated more peptide fragments from the globulin-type proteins than from the albumin-type. For instance, 93% of the hemp peptides resulting from tryptic hydrolysis was generated from the globulin-type proteins (edestin 1–3 and 7S vicilin-like protein) (Table 1). A similar trend was observed in other studies of the in silico hydrolysis of oilseed proteins [26,27]. Overall, our results suggest that the globulin-type proteins of oilseeds are prominent sources of peptides that could remain intact in the shrimp digestive tract.

2.2. Toxicity of Peptides

A bottleneck in the success of drugs used in aquaculture is their toxicity [7]. Peptides to be adopted as a potential anti-AHPND strategy should ideally be non-toxic to shrimp, other aquatic organisms, or human consumers. Therefore, by using the online admetSAR and ToxinPred tools, we screened the 809 peptides liberated from oilseed proteins by in silico hydrolysis for potential toxicity to crustaceans, fish, and humans. Remarkably, all the peptides were predicted as non-toxic to crustaceans (Table 2). Our results suggest that the five oilseed meals appear safe as protein ingredients in shrimp feed. This is in agreement with the report of good survival rate (>87%) in shrimp fed 25–100 g/kg rapeseed meal [28]. In other growth trials, the inclusion of 20–40% sesame meal and 1–10% sunflower cake to shrimp diets also led to good survival rates and health conditions [29,30]. Our results show that 25–35% of the peptides released from the five oilseeds were potentially toxic to fish (Table 2). Our results concur with previous observations of adverse effects in fish fed oilseed-containing diets at high inclusion levels. For instance, fish fed a diet containing 25% sesame meal as a replacement for soybean meal displayed comparable growth performance and survival rate to the control group; but a diet containing 75–100% sesame meal reduced growth performance and survival [31]. Besides a decline in weight gain, the specific growth rate and survival rate of fish was observed when 25–100% of soybean meal was substituted with pumpkin seed meal in fish feed [32]. Our results indicate that although all of the oilseed-derived peptides we investigated in this study were likely non-toxic to shrimp, some of them may be toxic to fish. The drugs used in aquaculture may enter environment, where they contaminate wild seafood [7]. Hence, in this study, only the 580 peptides that were predicted as being non-toxic to fish were considered in our subsequent analysis when searching for potential anti-AHPND peptides. Meanwhile, our ToxinPred analysis revealed that 96–100% of tryptic peptides generated from the oilseed proteins were non-toxic to human. This consideration is important as residual drugs in treated aquatic organisms, including shrimp, can pose health issues to people who consume aquatic food products [7]. Altogether, our toxicity screening narrowed down the initial pool of 809 peptides to a set of 572 oilseed-derived peptides which are potentially non-toxic to shrimp, fish and humans. Such peptides, 463 derived from globulin- and 109 from albumin-type proteins of the five oilseeds, were promising peptides for subsequent search of anti-AHPND agents in this study.

2.3. Water Solubility and Thermal Stability of Peptides

In shrimp farming, feed ingredients with high water solubility may be rapidly leached and lost on exposure to water, thus reducing their availability for shrimp ingestion [33]. Thus, in our search for potential anti-AHPND peptides as in-feed medication for shrimp, we further screened the 572 non-toxic oilseed-derived peptides for water solubility. This revealed 59 low-water-solubility peptides (data not shown), each made of 25% to 100% hydrophobic amino acids (e.g., Ala, Val, and Phe). This is in agreement with the finding that the nonpolar portion of peptides accounts for their low solubility in aqueous solutions [34]. In both feed processing and aquafarming, highly water-soluble compounds are prone to leaching. For instance, water-soluble crystalline amino acids were reported to leach out rapidly from shrimp feed, compromising the nutritional quality of shrimp diets [33]. Another study also pointed out that the leaching of water-soluble protein hydrolysates or peptides from feed may weaken the responses of fish to diets supplemented with such protein hydrolysates or peptides [35]. To minimize the potential issue of leaching, only the 59 low-water-solubility oilseed peptides were selected for further analysis in our search for in-feed anti-AHPND peptides. The 59 peptides included some repetitive sequences. For example, two identical sequences (i.e., GVLYK) could be found among the non-toxic hemp-derived peptides. By excluding the repetitive sequences, the 59 low-water-solubility peptides were consolidated into a set of 49 peptides with unique sequences, as shown in Table 3.
Thermal treatment is an important part of aquaculture feed production [36]. Hence, anti-AHPND peptides to be incorporated into shrimp feed should be heat-stable. In this study, the thermal stability of the 49 low-water-solubility oilseed peptides were predicted based on their aliphatic index and compared with lantibiotic nisin-A (Table 3). The aliphatic index, computed based on the percentage of aliphatic amino acids (Ala, Ile, Leu, and Val), is correlated with thermal stability to proteins [37]. Based on their aliphatic indices, 24 of the 49 peptides were likely thermally more stable than nisin-A (Table 3). This is notable as nisin-A is a heat-stable peptide that retains its antimicrobial activity after thermal treatment at 121 °C for 15 min; it has also been used as food preservative for decades [38]. The thermal stability predicted for the 24 oilseed-derived peptides could be attributed, at least in part, to a relatively high percentage of aliphatic amino acids in their sequences. For instance, PI and PVV, the sesame- and hemp-derived peptides predicted as having the highest aliphatic indices, are composed of 50% and 67% of aliphatic amino acids, respectively (Table 3). Taken together, our analysis led to 24 oilseed-derived peptides with favorable properties namely non-toxicity, low water solubility, and high thermal stability.

2.4. Molecular Docking Analysis

Interception of the complex formation between PirAvp and PirBvp was recommended as a promising anti-AHPND strategy [24]. To discover peptides that could preclude the formation of PirAvp/PirBvp complex, the 24 aforementioned peptides were each docked onto PirAvp and PirBvp. When docking onto PirAvp, 18 peptides produced similar or more negative docking scores (−144.315 to −194.881) than the six PirBvp regions (−142.536 to −174.899) which we used for comparison (Table 4). Hence, when compared with PirBvp, the 18 peptides potentially bound more stably to PirAvp. This implies the potential of the 18 peptides to intercept the binding between PirAvp-PirBvp. Intermolecular interactions between the 18 peptides and PirAvp were analyzed by using LigPlot+. The 18 peptides could each form interactions with 4–10 residues within the two PirBvp-binding regions of PirAvp (Table 4 and Table S1). For example, sunflower-derived VIQNLPNQCDLEVQQCTTCTG could bind to one residue (Val23) in the 15-WTVEPNGGVTEVDSKHTPIIPEVGRS-40 region, and to five residues (Thr52, Gln54, Ser71, Gln75, and Arg76) in the 52-TIQYQWGAPFMAGGWKVAKSHVVQRDET-79 region of PirAvp. Overall, most of the predicted peptide-PirAvp binding involved the 52-TIQYQWGAPFMAGGWKVAKSHVVQRDET-79 region of PirAvp. This region was also responsible for most of the predicted binding between the six PirBvp regions and PirAvp (Supplementary Materials Table S1). Overall, hydrophobic interactions formed the majority of the peptide-PirAvp interactions. A similar trend was observed when the six PirBvp regions were docked onto PirAvp (Table 4). In particular, sesame-derived ISTINSQTLPILSQLR stood out among the 18 peptides due to its largest number of interactions with residues in 52-TIQYQWGAPFMAGGWKVAKSHVVQRDET-79, including Lys67 and Lys70 (Table 4 and Table S1). The two Lys residues were found to be localized in the dimeric interface between PirAvp and PirBvp. The binding of peptides to the two Lys residues may block the cross-linking between PirAvp with PirBvp [24]. A graphical representation of the intermolecular interactions between ISTINSQTLPILSQLR and PirAvp is shown in Figure 1.
When the 24 non-toxic, low-water-solubility, and thermally-stable oilseed peptides were docked onto PirBvp, only six produced more negative docking scores, suggesting more stable binding, than the PirBvp-binding region 15-WTVEPNGGVTEVDSKHTPIIPEVGRS-40 of PirAvp. By contrast, all of the 24 peptides could not form a more stable binding to PirBvp when compared with the PirAvp region 52-TIQYQWGAPFMAGGWKVAKSHVVQRDET-79 (Table 5). A LigPlot+ analysis also revealed that the six best-scored peptides could bind to 2–4 residues within only two of the six PirAvp-binding regions of PirBvp, namely 386-FVVGENSGKPSVRLQL-401 and 426-YELFHPDEF-434. For instance, rape-derived QSLGVPPQLGNACNLDNLDVLQPTETIK, predicted to form the largest number of interactions with PirBvp, could bind to Val397, Arg398, and Gln400 in the 386-FVVGENSGKPSVRLQL-401 region, and could bind to Phe429 in the 426-YELFHPDEF-434 region of PirBvp. Figure 2 shows the intermolecular interactions contributing to QSLGVPPQLGNACNLDNLDVLQPTETIK-PirBvp binding. Similarly, limited interactions between the two PirAvp regions 15-WTVEPNGGVTEVDSKHTPIIPEVGRS-40 and 52-TIQYQWGAPFMAGGWKVAKSHVVQRDET-79 with PirBvp were also observed (Table 5 and Table S2). Interestingly, the six oilseed peptides having predicted interactions with PirBvp (Table 5) were also among the 18 best-scoring PirAvp-binding peptides presented in Table 4. The six peptides ranged between 1139 Da and 2977 Da in mass, and 10–28 residues in length, mostly having a net charge of +1. Modeling with PEP-FOLD 3.5 revealed the 3D structures of all the peptides except the rape-derived ISYVVQGMGISGR to be comprised of helical elements (Figure 3). The six peptides were derived from globulin-type proteins from rape, sesame, and hemp. Our findings suggest that the globulin-type proteins of these three seeds are potential sources of peptides that can bind to both PirAvp and PirBvp. As observed in oilseed peptide-PirAvp interactions (Table 4), hydrophobic interactions contributed to the majority of the interactions between the six peptides and PirBvp (Table 5). Taken together, our results suggest that hydrophobic interactions are critical in facilitating stable binding between the six oilseed peptides and the proteins of both PirAvp and PirBvp.
The virulence factors of AHPND have been shown to be both PirAvp and PirBvp [9,39]. Hence, targeting only one of the proteins may not result in full protection against AHPND. For example, the binding of the chicken egg yolk immunoglobulin (IgY) specific for a recombinant PirBvp-like protein (anti-PirBvp-IgY) to the PirBvp toxin still allowed AHPND infection to progress, resulting in a 14% survival rate of AHPND-challenged shrimp [40]. On the other hand, the binding of the IgY specific for the recombinant PirAvp-like protein (anti-PirAvp-IgY) to PirAvp gave rise to only an 86% survival rate in AHPND-challenged shrimp [40]. Therefore, we speculate that dual-target oilseed peptides that can bind to both PirAvp and PirBvp may be a relatively promising approach for intercepting PirAvp/PirBvp complex formation, thus mitigating AHPND occurrence in shrimp. Based on our results, the six oilseed peptides that can bind to both PirAvp and PirBvp (ISYVVQGMGISGR, LTFVVHGHALMGK, QSLGVPPQLGNACNLDNLDVLQPTETIK, ISTINSQTLPILSQLR, PQFLVGASSILR, and VQVVNHMGQK) deserve attention in future research. In this theoretical investigation, we considered peptides undegraded by trypsin to be bioavailable in the shrimp digestive tract. Considering the complexity of factors governing bioavailability, investigation into the in vivo availability of the six aforementioned peptides should be the first priority in future research. Next, testing of the in vivo effectiveness of these peptides individually and in combination in live shrimp, delivered as ingredients incorporated into shrimp feed, would be of great interest. Besides validating the properties predicted in this study (absence of toxicity, low water solubility, and thermal stability) and confirming the anti-AHPND potency of the peptides, such studies may also reveal whether the peptides could function synergistically. An interesting recent development is the novel discovery by Almanza-Martínez et al. [41] of an α-amylase-like protein serving as a binding receptor to PirBvp in the digestive tract tissue of the whiteleg shrimp. The specific role of the α-amylase-like protein in the mechanism of action of PirBvp remains to be established, as are the residues participating in the α-amylase-PirBvp binding interface. Once such information has been elucidated, it would be interesting to investigate, through molecular docking and/or in vitro experiments, whether the oilseed peptides reported in this study can also intercept the α-amylase-PirBvp interaction. Meanwhile, the potential off-target interactions of the six selected oilseed peptides should also be addressed to ensure that they do not disrupt other normal biological processes, causing adverse side-effects in the shrimp.

3. Materials and Methods

3.1. In Silico Hydrolysis of Oilseed Meal Proteins

A total of 37 major proteins found in five oilseed meals, namely those of hemp [42,43], pumpkin [44], rape [45], sesame [20], and sunflower [46] were chosen for this study. The major proteins chosen can be classified into two types, globulins and albumins, which account for more than 78% of seed protein constituents [43,44,45,46,47]. The sequence of each protein was downloaded from the UniProt Knowledgebase (https://www.uniprot.org/) [48] on 5 May 2021. All 37 proteins were hydrolyzed by trypsin (EC 3.4.21.4) using the “Enzyme(s) action” tool on the BIOPEP-UWM bioactive peptides database (http://www.uwm.edu.pl/biochemia/index.php/en/biopep) [49] on 23 May 2021.

3.2. Toxicity of Peptides

Aquatic toxicity of peptides towards crustaceans and fish was analyzed by using admetSAR (http://lmmd.ecust.edu.cn/admetsar2) [50]. Peptide sequences were converted into the molecular input line entry system (SMILES) format using BIOPEP-UWM (http://www.uwm.edu.pl/biochemia/index.php/en/biopep) [49] before analysis with admetSAR. Toxicity of peptides to humans was predicted with ToxinPred (https://webs.iiitd.edu.in/raghava/toxinpred/multi_submit.php) [51]. Support vector machine (SVM) and the default SVM threshold of 0.0 were chosen for ToxinPred toxicity prediction. Peptides with an SVM score < 0.0 were predicted as non-toxic. The aforementioned analyses were performed between 25 May and 31 August 2021.

3.3. Water Solubility and Thermal Stability of Peptides

Water solubility of peptides were analyzed using an online peptide property calculator (https://pepcalc.com/). Aliphatic index, an indicator of thermal stability, was predicted by using the Peptides Package in R (https://rdrr.io/snippets/) [37]. The two online tools were accessed between 26 to 31 August 2021.

3.4. Molecular Docking Analysis

Oilseed peptides predicted as non-toxic to crustaceans, fish and humans, with low water solubility and high thermal stability, were chosen for molecular docking analysis. The two target proteins from V. parahaemolyticus, namely PirAvp (Protein Data Bank (PDB) code: 3X0T) and PirBvp (PDB code: 3X0U) [9] were used as receptors. The docking of peptides onto chain A of PirAvp and of PirBvp was performed by submitting PDB ID: Chain ID information and peptide sequences to HPEPDOCK (http://huanglab.phys.hust.edu.cn/hpepdock/) [52]. For comparison, two PirAvp regions (15-WTVEPNGGVTEVDSKHTPIIPEVGRS-40 and 52-TIQYQWGAPFMAGGWKVAKSHVVQRDET-79), and six PirBvp regions (214-WADNDSYNNANQD-226, 290-DEIPQPLKPNM-300, 322-YNRVGRLKL-330, 386-FVVGENSGKPSVRLQL-401, 409-MLADQEGSDKVAA-421, and 426-YELFHPDEF-434) reported to be involved in PirAvp/PirBvp complex formation [24], were used. These sequences were docked onto each corresponding target protein as described above. All docking analyses were carried out between 27 May 2021 and 31 August 2021. The best (most negative) docking score of each peptide was tabulated. Three-dimensional (3D) structures of the protein-peptide docked models were visualized 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). Two-dimensional (2D) protein-peptide interaction diagrams were generated with LigPlot+ v.2.2 [53,54], through which intermolecular interactions were visualized and analyzed. PEP-FOLD 3.5 (https://bioserv.rpbs.univ-paris-diderot.fr/services/PEP-FOLD3/) [55,56,57] was used to model the 3D structures of selected oilseed peptides, setting the number of simulations to 200 and with model sorting based on sOPEP energy. PepDraw (http://pepdraw.com/) was used to compute the molecular masses and net charges of the selected peptides. PEP-FOLD 3.5 and PepDraw were accessed on 1 September 2021.

4. Conclusions

Our in silico study found the globulin-type proteins of the oilseeds investigated to be promising sources of PirAvp- and PirBvp- binding peptides. Upon selection for non-toxicity against crustaceans, fish, and humans, as well as thermal stability and low water solubility, six peptides that could bind to both PirAvp and PirBvp with stability comparable or superior to the binding between the key regions in the native PirAvp/PirBvp complex, were identified. These six peptides could potentially play a role in intercepting the formation of the toxic PirAvp/PirBvp complex, thus serving as promising anti-AHPND agents. The findings of this study can be used to guide the future discovery and development of anti-AHPND peptides from oilseed proteins.

Supplementary Materials

The following are available online at https://www.mdpi.com/article/10.3390/antibiotics10101211/s1, Table S1: Number of residues in the two PirBvp-binding regions of PirAvp protein that are involved in interactions with 18 chosen oilseed peptides, Table S2: Number of residues in the six PirAvp-binding regions of PirBvp protein that are involved in interactions with six chosen oilseed peptides.

Author Contributions

Conceptualization, W.-L.W. and T.-T.C.; methodology and software, F.-C.W. and T.-T.C.; formal analysis, investigation, data curation, project administration, J.-H.O.; writing—original draft preparation, J.-H.O. and T.-T.C.; writing—review and editing, W.-L.W. and F.-C.W.; supervision, W.-L.W., F.-C.W. and T.-T.C. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Data Availability Statement

The data presented in this study are available on request from the corresponding author.

Acknowledgments

We thank Hongbin Yang and Chaofeng Lou for their technical assistance in our admetSAR analysis.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Tandel, G.M.; John, K.R.; Rosalind George, M.; Prince Jeyaseelan, M.J. Current status of viral diseases in Indian shrimp aquaculture. Acta Virol. 2017, 61, 131–137. [Google Scholar] [CrossRef]
  2. Kmmari, S.; Rathlavath, S.; Pillai, D.; Rajesh, G. Hepatopancreatic microsporidiasis (HPM) in shrimp culture: A review. Int. J. Curr. Microbiol. App. Sci. 2018, 7, 3208–3215. [Google Scholar] [CrossRef] [Green Version]
  3. Zorriehzahra, M.J.; Banaederakhshan, R. Early mortality syndrome (EMS) as new emerging threat in shrimp industry. Adv. Anim. Vet. Sci. 2015, 3, 64–72. [Google Scholar] [CrossRef] [Green Version]
  4. Santos, H.M.; Tsai, C.-Y.; Maquiling, K.R.A.; Tayo, L.L.; Mariatulqabtiah, A.R.; Lee, C.-W.; Chuang, K.P. Diagnosis and potential treatments for acute hepatopancreatic necrosis disease (AHPND): A review. Aquac. Int. 2020, 28, 169–185. [Google Scholar] [CrossRef]
  5. Dong, X.; Bi, D.; Wang, H.; Zou, P.; Xie, G.; Wan, X.; Yang, Q.; Zhu, Y.; Chen, M.; Guo, C.; et al. pirABvp-bearing Vibrio parahaemolyticus and Vibrio campbellii pathogens isolated from the same AHPND-affected pond possess highly similar pathogenic plasmids. Front. Microbiol. 2017, 8, 1859. [Google Scholar] [CrossRef]
  6. Han, J.E.; Mohney, L.L.; Tang, K.F.; Pantoja, C.R.; Lightner, D.V. Plasmid mediated tetracycline resistance of Vibrio parahaemolyticus associated with acute hepatopancreatic necrosis disease (AHPND) in shrimps. Aquac. Rep. 2015, 2, 17–21. [Google Scholar] [CrossRef] [Green Version]
  7. Okocha, R.C.; Olatoye, I.O.; Adedeji, O.B. Food safety impacts of antimicrobial use and their residues in aquaculture. Public Health Rev. 2018, 39, 21. [Google Scholar] [CrossRef] [PubMed]
  8. Han, J.E.; Tang, K.F.J.; Tran, L.H.; Lightner, D.V. Photorhabdus insect-related (Pir) toxin-like genes in a plasmid of Vibrio parahaemolyticus, the causative agent of acute hepatopancreatic necrosis disease (AHPND) of shrimp. Dis. Aquat. Org. 2015, 113, 33–40. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  9. Lee, C.-T.; Chen, I.-T.; Yang, Y.-T.; Ko, T.-P.; Huang, Y.-T.; Huang, J.-Y.; Huang, M.-F.; Lin, S.-J.; Chen, C.-Y.; Lin, S.-S.; et al. The opportunistic marine pathogen Vibrio parahaemolyticus becomes virulent by acquiring a plasmid that expresses a deadly toxin. Proc. Natl. Acad. Sci. USA 2015, 112, 10798–10803. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  10. Tran, L.; Nunan, L.; Redman, R.M.; Mohney, L.L.; Pantoja, C.R.; Fitzsimmons, K.; Lightner, D.V. Determination of the infectious nature of the agent of acute hepatopancreatic necrosis syndrome affecting penaeid shrimp. Dis. Aquat. Org. 2013, 105, 45–55. [Google Scholar] [CrossRef]
  11. Lin, S.-J.; Hsu, K.-C.; Wang, H.-C. Structural insights into the cytotoxic mechanism of Vibrio parahaemolyticus PirAvp and PirBvp toxins. Mar. Drugs 2017, 15, 373. [Google Scholar] [CrossRef] [Green Version]
  12. Singh, M.; Singh, P. Drugs and chemicals applied in aquaculture industry: A review of commercial availability, recommended dosage and mode of application. J. Entomol. Zool. Stud. 2018, 6, 903–907. [Google Scholar]
  13. Gómez-Jimenez, S.; Espinosa-Plascencia, A.; Valenzuela-Villa, F.; Bermúdez-Almada, M.D.C. Oxytetracycline (OTC) accumulation and elimination in hemolymph, muscle and hepatopancreas of white shrimp Litopenaeus vannamei following an OTC-feed therapeutic treatment. Aquaculture 2008, 274, 24–29. [Google Scholar] [CrossRef]
  14. Pinoargote, G.; Flores, G.; Cooper, K.; Ravishankar, S. Effects on survival and bacterial community composition of the aquaculture water and gastrointestinal tract of shrimp (Litopenaeus vannamei) exposed to probiotic treatments after an induced infection of acute hepatopancreatic necrosis disease. Aquac. Res. 2018, 49, 3270–3288. [Google Scholar] [CrossRef] [Green Version]
  15. Fuandila, N.N.; Widanarni, W.; Yuhana, M. Growth performance and immune response of prebiotic honey fed pacific white shrimp Litopenaeus vannamei to Vibrio parahaemolyticus infection. J. Appl. Aquac. 2020, 32, 221–235. [Google Scholar] [CrossRef]
  16. Jha, R.K.; Babikian, Y.H.; Babikian, H.Y.; Khoa, L.V.; Wisoyo, D.; Srisombat, S.; Jiaravanon, B. Efficacy of natural herbal formulation against acute hepatopancreatic necrosis disease (AHPND) causing Vibrio parahaemolyticus in Penaeus vannamei. Open J. Vet. Med. 2017, 2, 1–6. [Google Scholar] [CrossRef]
  17. Wang, G.; Yu, E.; Li, Z.; Yu, D.; Wang, H.; Gong, W. Effect of bioactive peptides (BPs) on the development of Pacific white shrimp (Litopenaeus vannamei Boone, 1931). J. Ocean Univ. China 2016, 15, 495–501. [Google Scholar] [CrossRef]
  18. Tang, C.-H.; Wang, X.-S.; Yang, X.-Q. Enzymatic hydrolysis of hemp (Cannabis sativa L.) protein isolate by various proteases and antioxidant properties of the resulting hydrolysates. Food Chem. 2009, 114, 1484–1490. [Google Scholar] [CrossRef]
  19. Zhou, C.; Yu, X.; Qin, X.; Ma, H.; Yagoub, A.E.A.; Hu, J. Hydrolysis of rapeseed meal protein under simulated duodenum digestion: Kinetic modeling and antioxidant activity. LWT 2016, 68, 523–531. [Google Scholar] [CrossRef]
  20. Lu, X.; Zhang, L.; Sun, Q.; Song, G.; Huang, J. Extraction, identification and structure-activity relationship of antioxidant peptides from sesame (Sesamum indicum L.) protein hydrolysate. Food Res. Int. 2019, 116, 707–716. [Google Scholar] [CrossRef]
  21. Vaghei, R.G.; Abolhasani, M.H.; Ghorbani, R.; Matinfar, A. Production of soybean meal-based feed and its effect on growth performance of western white shrimp. Iran. J. Fish. Sci. 2017, 16, 578–586. [Google Scholar]
  22. Shao, J.; Zhao, W.; Han, S.; Chen, Y.; Wang, B.; Wang, L. Partial replacement of fishmeal by fermented soybean meal in diets for juvenile white shrimp (Litopenaeus vannamei). Aquac. Nutr. 2019, 25, 145–153. [Google Scholar] [CrossRef] [Green Version]
  23. Cheng, A.-C.; Lin, H.-L.; Shiu, Y.-L.; Tyan, Y.-C.; Liu, C.-H. Isolation and characterization of antimicrobial peptides derived from Bacillus subtilis E20-fermented soybean meal and its use for preventing Vibrio infection in shrimp aquaculture. Fish. Shellfish Immunol. 2017, 67, 270–279. [Google Scholar] [CrossRef]
  24. Lin, S.-J.; Chen, Y.-F.; Hsu, K.-C.; Chen, Y.-L.; Ko, T.-P.; Lo, C.-F.; Wang, H.-C.; Wang, H.-C. Structural insights to the heterotetrameric interaction between the Vibrio parahaemolyticus PirAvp and PirBvp toxins and activation of the Cry-Like pore-forming domain. Toxins 2019, 11, 233. [Google Scholar] [CrossRef] [Green Version]
  25. Lazo, J.P.; Romaire, R.P.; Reigh, R.C. Evaluation of three in vitro enzyme assays for estimating protein digestibility in the pacific white shrimp Penaeus vannamei. J. World Aquac. Soc. 1998, 29, 441–450. [Google Scholar] [CrossRef]
  26. Duan, X.; Zhang, M.; Chen, F. Prediction and analysis of antimicrobial peptides from rapeseed protein using in silico approach. J. Food Biochem. 2021, 45, e13598. [Google Scholar] [CrossRef] [PubMed]
  27. Wong, F.-C.; Ong, J.-H.; Chai, T.-T. SARS-CoV-2 spike protein-, main protease- and papain-like-protease-targeting peptides from seed proteins following gastrointestinal digestion: An in silico study. Phytomed. Plus 2021, 1, 100016. [Google Scholar] [CrossRef]
  28. Jannathulla, R.; Dayal, J.S.; Ambasankar, K.; Yuvapushpa, R.; Kumar, J.A.; Muralidhar, M. Evaluation of fungal fermented rapeseed meal as a fishmeal substitute in the diet of Penaeus vannamei. J. Coast. Res. 2019, 86, 82–89. [Google Scholar] [CrossRef]
  29. Dayal, J.S.; Rajaram, V.; Ambasankar, K.; Ali, S.A. Sunflower oil cake as a replacement for fish meal in feeds of tiger shrimp, Penaeus monodon reared in tanks and in net cages. Indian J. Geo-Mar. Sci. 2011, 40, 460–470. [Google Scholar]
  30. Bae, J.; Hamidoghli, A.; Djaballah, M.S.; Maamri, S.; Hamdi, A.; Souffi, I.; Farris, N.W.; Bai, S.C. Effects of three different dietary plant protein sources as fishmeal replacers in juvenile whiteleg shrimp, Litopenaeus vannamei. Fish. Aquat. Sci. 2020, 23, 2. [Google Scholar] [CrossRef] [Green Version]
  31. Jimoh, W.; Aroyehun, H.T. Evaluation of cooked and mechanically defatted sesame (Sesamum indicum) seed meal as a replacer for soybean meal in the diet of African catfish (Clarias gariepinus). Turk. J. Fish. Aquat. Sci. 2011, 11, 185–190. [Google Scholar] [CrossRef]
  32. Omosowone, O.O.; Ogunrinde, A.M. Effect of partial or total supplementation of soybean meal with fluted pumpkin (Telfairia occidentalis) seed meal in the diet of hybrid Catfish (Heteroclarias) fingerlings. Food Sci. Nutr. 2018, 6, 1190–1195. [Google Scholar] [CrossRef]
  33. Simon, C.J.; Truong, H.; Habilay, N.; Hines, B. Feeding behaviour and bioavailability of essential amino acids in shrimp Penaeus monodon fed fresh and leached fishmeal and fishmeal-free diets. Animals 2021, 11, 847. [Google Scholar] [CrossRef]
  34. Sarma, R.; Wong, K.-Y.; Lynch, G.C.; Pettitt, B.M. Peptide solubility limits: Backbone and side-chain interactions. J. Phys. Chem. B. 2018, 122, 3528–3539. [Google Scholar] [CrossRef]
  35. Kvåle, A.; Yúfera, M.; Nygård, E.; Aursland, K.; Harboe, T.; Hamre, K. Leaching properties of three different micropaticulate diets and preference of the diets in cod (Gadus morhua L.) larvae. Aquaculture 2006, 251, 402–415. [Google Scholar] [CrossRef]
  36. Mo, W.Y.; Man, Y.B.; Wong, M.H. Use of food waste, fish waste and food processing waste for China’s aquaculture industry: Needs and challenge. Sci. Total Environ. 2018, 613–614, 635–643. [Google Scholar] [CrossRef]
  37. Osorio, D.; Rondón-Villarreal, P.; Torres, R. Peptides: A package for data mining of antimicrobial peptides. R J. 2015, 7, 4–14. [Google Scholar] [CrossRef]
  38. Delves-Broughton, J.; Weber, G. Nisin, natamycin and other commercial fermentates used in food biopreservation. In Protective Cultures, Antimicrobial Metabolites and Bacteriophages for Food and Beverage Biopreservation; Elsevier: Amsterdam, The Netherlands, 2011; pp. 63–99. [Google Scholar]
  39. Phiwsaiya, K.; Charoensapsri, W.; Taengphu, S.; Dong, H.T.; Sangsuriya, P.; Nguyen, G.T.T.; Pham, H.Q.; Amparyup, P.; Sritunyalucksana, K.; Taengchaiyaphum, S.; et al. A natural Vibrio parahaemolyticus ΔpirAVppirBVp+ mutant kills shrimp but produces neither Pir Vp toxins nor acute hepatopancreatic necrosis disease lesions. Appl. Environ. Microbiol. 2017, 83, e00680-17. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  40. Nakamura, R.; Pedrosa-Gerasmio, I.R.; Alenton, R.R.R.; Nozaki, R.; Kondo, H.; Hirono, I. Anti-PirA-like toxin immunoglobulin (IgY) in feeds passively immunizes shrimp against acute hepatopancreatic necrosis disease. J. Fish. Dis. 2019, 42, 1125–1132. [Google Scholar] [CrossRef]
  41. Almanza-Martínez, N.; Francisco Martínez Díaz, S.; Flores-Ramírez, G.; Zuñiga-Navarrete, F.; Gómez, I.; Cardona-Félix, C.S. An α-amylase-like protein interacts with PirB toxin from Vibrio parahaemolyticus in digestive tract tissue of white shrimp Litopenaeus vannamei. Aquac. Res. 2020, 51, 3910–3914. [Google Scholar] [CrossRef]
  42. Kotecka-Majchrzak, K.; Kasałka-Czarna, N.; Sumara, A.; Fornal, E.; Montowska, M. Multispecies identification of oilseed- and meat-specific proteins and heat-stable peptide markers in food products. Molecules 2021, 26, 1577. [Google Scholar] [CrossRef] [PubMed]
  43. Mamone, G.; Picariello, G.; Ramondo, A.; Nicolai, M.A.; Ferranti, P. Production, digestibility and allergenicity of hemp (Cannabis sativa L.) protein isolates. Food Res. Int. 2019, 115, 562–571. [Google Scholar] [CrossRef] [PubMed]
  44. Rezig, L.; Chibani, F.; Chouaibi, M.; Dalgalarrondo, M.; Hessini, K.; Guéguen, J.; Hamdi, S. Pumpkin (Cucurbita maxima) seed proteins: Sequential extraction processing and fraction characterization. J. Agric. Food Chem. 2013, 61, 7715–7721. [Google Scholar] [CrossRef]
  45. Perera, S.P.; McIntosh, T.C.; Wanasundara, J.P.D. Structural properties of cruciferin and napin of Brassica napus (Canola) show distinct responses to changes in pH and temperature. Plants 2016, 5, 36. [Google Scholar] [CrossRef] [Green Version]
  46. Wildermuth, S.R.; Young, E.E.; Were, L.M. Chlorogenic acid oxidation and its reaction with sunflower proteins to form green-colored complexes. Compr. Rev. Food Sci. Food Saf. 2016, 15, 829–843. [Google Scholar] [CrossRef]
  47. Onsaard, E. Sesame proteins. Int. Food Res. J. 2012, 19, 1287–1295. [Google Scholar]
  48. The UniProt Consortium. UniProt: A worldwide hub of protein knowledge. Nucleic Acids Res. 2019, 47, D506–D515. [Google Scholar] [CrossRef] [Green Version]
  49. Minkiewicz, P.; Iwaniak, A.; Darewicz, M. BIOPEP-UWM database of bioactive peptides: Current opportunities. Int. J. Mol. Sci. 2019, 20, 5978. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  50. Liu, L.; Yang, H.; Cai, Y.; Cao, Q.; Sun, L.; Wang, Z.; Li, W.; Liu, G.; Lee, P.W.; Tang, Y. In silico prediction of chemical aquatic toxicity for marine crustaceans via machine learning. Toxicol. Res. 2019, 8, 341–352. [Google Scholar] [CrossRef]
  51. Gupta, S.; Kapoor, P.; Chaudhary, K.; Gautam, A.; Kumar, R.; Raghava, G.P.S. In silico approach for predicting toxicity of peptides and proteins. PLoS ONE 2013, 8, e73957. [Google Scholar]
  52. Zhou, P.; Jin, B.; Li, H.; Huang, S.-Y. HPEPDOCK: A web server for blind peptide-protein docking based on a hierarchical algorithm. Nucleic Acids Res. 2018, 46, W443–W450. [Google Scholar] [CrossRef]
  53. Wallace, A.C.; Laskowski, R.A.; Thornton, J.M. LIGPLOT: A program to generate schematic diagrams of protein-ligand interactions. Protein Eng. 1995, 8, 127–134. [Google Scholar] [CrossRef] [PubMed]
  54. Laskowski, R.A.; Swindells, M.B. LigPlot+: Multiple ligand-protein interaction diagrams for drug discovery. J. Chem. Inf. Model. 2011, 51, 2778–2786. [Google Scholar] [CrossRef] [PubMed]
  55. Lamiable, A.; Thévenet, P.; Rey, J.; Vavrusa, M.; Derreumaux, P.; Tufféry, P. PEP-FOLD3: Faster de novo structure prediction for linear peptides in solution and in complex. Nucleic Acids Res. 2016, 44, W449–W454. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  56. Shen, Y.; Maupetit, J.; Derreumaux, P.; Tufféry, P. Improved PEP-FOLD approach for peptide and miniprotein structure prediction. J. Chem. Theory Comput. 2014, 10, 4745–4758. [Google Scholar] [CrossRef]
  57. Thévenet, P.; Shen, Y.; Maupetit, J.; Guyon, F.; Derreumaux, P.; Tufféry, P. PEP-FOLD: An updated de novo structure prediction server for both linear and disulfide bonded cyclic peptides. Nucleic Acids Res. 2012, 40, W288–W293. [Google Scholar] [CrossRef] [Green Version]
Antibiotics 10 01211 g001 550
Figure 1. The ISTINSQTLPILSQLR-PirAvp docked model presented in 3D (A,B) and 2D (C) diagrams. In (A,B), PirAvp is shown in the cartoon in maroon; ISTINSQTLPILSQLR is presented in a ball-and-stick style. In the 2D diagram (C), the hydrophobic bonds and hydrogen bonds are displayed in red spoked arcs and green dashed lines, respectively.
Figure 1. The ISTINSQTLPILSQLR-PirAvp docked model presented in 3D (A,B) and 2D (C) diagrams. In (A,B), PirAvp is shown in the cartoon in maroon; ISTINSQTLPILSQLR is presented in a ball-and-stick style. In the 2D diagram (C), the hydrophobic bonds and hydrogen bonds are displayed in red spoked arcs and green dashed lines, respectively.
Antibiotics 10 01211 g001
Antibiotics 10 01211 g002 550
Figure 2. The QSLGVPPQLGNACNLDNLDVLQPTETIK-PirBvp docked model presented in 3D (A,B) and 2D (C) diagrams. In (A,B), PirBvp is shown in the cartoon in green; QSLGVPPQLGNACNLDNLDVLQPTETIK is presented in a ball-and-stick style. In the 2D diagram (C), the hydrophobic bonds, hydrogen bonds, and external bonds are displayed in red spoked arcs, green dashed lines, and purple lines, respectively.
Figure 2. The QSLGVPPQLGNACNLDNLDVLQPTETIK-PirBvp docked model presented in 3D (A,B) and 2D (C) diagrams. In (A,B), PirBvp is shown in the cartoon in green; QSLGVPPQLGNACNLDNLDVLQPTETIK is presented in a ball-and-stick style. In the 2D diagram (C), the hydrophobic bonds, hydrogen bonds, and external bonds are displayed in red spoked arcs, green dashed lines, and purple lines, respectively.
Antibiotics 10 01211 g002
Antibiotics 10 01211 g003 550
Figure 3. The 3D structures of six peptides predicted to bind to the key regions of both PirAvp and PirBvp that are involved in PirAvp/PirBvp complex formation: (A) ISYVVQGMGISGR, (B) LTFVVHGHALMGK, (C) QSLGVPPQLGNACNLDNLDVLQPTETIK, (D) ISTINSQTLPILSQLR, (E) PQFLVGASSILR, and (F) VQVVNHMGQK. The alpha-helical elements are shown in maroon, whereas the beta-strands are shown in light blue. Only the N- and C-terminal residues of each peptide are indicated. Values below each structure represent the molecular mass and the net charge of the peptide, as computed by PepDraw.
Figure 3. The 3D structures of six peptides predicted to bind to the key regions of both PirAvp and PirBvp that are involved in PirAvp/PirBvp complex formation: (A) ISYVVQGMGISGR, (B) LTFVVHGHALMGK, (C) QSLGVPPQLGNACNLDNLDVLQPTETIK, (D) ISTINSQTLPILSQLR, (E) PQFLVGASSILR, and (F) VQVVNHMGQK. The alpha-helical elements are shown in maroon, whereas the beta-strands are shown in light blue. Only the N- and C-terminal residues of each peptide are indicated. Values below each structure represent the molecular mass and the net charge of the peptide, as computed by PepDraw.
Antibiotics 10 01211 g003
Table
Table 1. Numbers of peptide fragments released from major oilseed proteins by in silico trypsin hydrolysis.
Table 1. Numbers of peptide fragments released from major oilseed proteins by in silico trypsin hydrolysis.
SourceProteinUniProt
Accession		Type *Number of ResiduesNumber of Peptides
HempAlbuminA0A219D1L6A11918
Edestin 1A0A090DLH8G48862
Edestin 2A0A090DLI7G46756
Edestin 3A0A219D3H6G46854
7S vicilin-like proteinA0A219D1T7G47257
Pumpkin2S albumin large chainQ39649A6713
11S globulin delta chainP13744G18424
11S globulin gamma chainP13744G27532
RapeCruciferin BnC1 subunit alphaP33523G27724
Cruciferin BnC1 subunit betaP33523G19016
Cruciferin BnC2 subunit alphaP33524G28322
Cruciferin BnC2 subunit betaP33524G19019
Cruciferin CRU1 alpha chainP33525G29622
Cruciferin CRU1 beta chainP33525G19017
Cruciferin CRU4 alpha chainP33522G25420
Cruciferin CRU4 beta chainP33522G18917
Cruciferin subunit alphaP11090G27524
Cruciferin subunit betaP11090G19016
Napin-1A small chainP24565A315
Napin-1A large chainP24565A798
Napin-2 small chainP01090A366
Napin-2 large chainP01090A819
Napin-3 small chainP80208A376
Napin-3 large chainP80208A889
Napin-B small chainP27740A366
Napin-B large chainP27740A8410
Napin embryo-specific small chainP09893A386
Napin embryo-specific large chainP09893A8910
Sesame2S seed storage protein 1 small subunitQ9XHP1A306
2S seed storage protein 1 large subunitQ9XHP1A7010
11S globulin isoform 3Q2XSW7G46854
11S globulin isoform 4Q2XSW6G44946
11S globulin seed storage protein 2 acidic chainQ9XHP0G25631
11S globulin seed storage protein 2 basic chainQ9XHP0G18218
Sunflower2S seed storage proteinP15461A13415
11S globulin seed storage protein G3 acidic chainP19084G28522
11S globulin seed storage protein G3 basic chainP19084G18819
* G, globulin-type; A, albumin-type.
Table
Table 2. In silico prediction of crustacean, fish and human toxicity of oilseed-derived peptides.
Table 2. In silico prediction of crustacean, fish and human toxicity of oilseed-derived peptides.
SourceProtein Type *Total
Number of Peptides		Number of Peptides Predicted as Non-Toxic (NT) and Toxic (T) against Different Organisms
CrustaceanFishHuman
NTTNTTNTT
HempG2292290170592290
A18180162171
PumpkinG565604016560
A13130121121
RapeG1971970115821940
A757506114705
SesameG1491490109401490
A16160151142
SunflowerG414102912400
A15150132141
Total 809809058022979510
* G, globulin-type; A, albumin-type.
Table
Table 3. Aliphatic indices of 49 oilseed-derived peptides predicted to be non-toxic and having low water solubility, in comparison with reference peptide lantibiotic nisin-A.
Table 3. Aliphatic indices of 49 oilseed-derived peptides predicted to be non-toxic and having low water solubility, in comparison with reference peptide lantibiotic nisin-A.
SourcePeptideAliphatic Index
SesamePI195
GLIVMAR167
GHIITVAR146
ISTINSQTLPILSQLR146
IQVVGHK139
GVLYR136
GLQVISPPLQR133
VASA123
QEQFQCAGIVAMR68
MTFVR58
QTFHNIFR49
YWQSLQQHQQHR33
GQHQFGNVFR29
TGGYA20
HCMQWMR0
GSTWQQGQCR0
HempPVV193
GVLYK136
PQFLVGASSILR130
QGQIVTVPQNHAVVK110
ESVILPTSAASPPVK104
LGNLTSYQR87
VQVVNHMGQK87
ATA67
NIPSMCGMQPR35
TTWSWR0
WQSQCQFQR0
SunflowerGHIVNVGQDLQIVR146
VIQNLPNQCDLEVQQCTTCTG83
WVSFK58
GGWSN0
RapeLTFVVHGHALMGK112
QSLGVPPQLGNACNLDNLDVLQPTETIK108
ISYVVQGMGISGR107
TNANAQINTLAGR83
NLPNVCNMK76
TNANAMVSTLAGR75
CSGVSFVR73
ACQQWIR70
ACQQWLHK61
ATSQQFQWIEFK41
QQQGQQMQGQQMQQVISR38
VQGQHGPFQSTR24
QAMQSGGG13
QAMQSGSG13
QAMQPGGGSG10
TMPG0
TMPGPS0
TMPGPSY0
Lantibiotic nisin-AITSISLCTPGCKTGALMGCNMKTATCHCSIHVSK72
Table
Table 4. Docking scores for interactions between the 18 best-scored oilseed peptides and PirAvp, with details on the PirAvp residues involved in different types of interactions.
Table 4. Docking scores for interactions between the 18 best-scored oilseed peptides and PirAvp, with details on the PirAvp residues involved in different types of interactions.
SourcePeptideDocking ScoreInteraction with PirAvp Residues *
Hydrogen BondHydrophobic InteractionSalt BridgeExternal Bond
RapeLTFVVHGHALMGK−194.881Asn87Trp57, Gly58, Pro60, Ala63, Ala69, Lys70, Tyr80, Gln83, Pro85, Asn87, Ala88, Phe89, Tyr90--
ISYVVQGMGISGR−174.661Tyr80(2)Trp57, Gly58, Ala59, Pro60, Val68, Ala69, Lys70, Ser71, His72, Tyr80, His81, Leu82, Gln83, Pro85, Asn87, Phe89, Tyr90--
QSLGVPPQLGNACNLDNLDVLQPTETIK−168.015Gln54, Trp57, Ser71(2), Tyr80, Asn87Thr52, Gln54, Trp57, Val68, Ala69, Lys70, Ser71, His72, Tyr80, His81, Pro85, Asn87, Phe89, Asn99-Asn99(3)
TNANAMVSTLAGR−156.242Ala69, Ser71Pro35, Trp57, Val68, Ala69, Lys70, Ser71, His72, Asn87, Tyr80, Ala88, Tyr90--
TNANAQINTLAGR−151.751Ala69, His72, Tyr80, Asn87(2)Trp57, Ala59, Pro60, Val68, Ala69, Lys70, Ser71, His72, Tyr80, Pro85, Asn87, Phe89, Tyr90--
NLPNVCNMK−147.564Lys29, Arg39(3), Arg84Tyr11, Lys29, Gly38, Arg39, Ser40, Arg84, Asp86, His111, Leu112, Glu113, His114, His115Asp86Arg84
CSGVSFVR−146.934Asn87Trp57, Gly58, Ala59, Pro60, Ala69, Lys70, Pro85, Asp86, Asn87, Ala88, Phe89, Tyr90-Asp86
SesameISTINSQTLPILSQLR−178.623Val68, Ser71, His72, Asn87Thr52, Trp57, Gly58, Lys67, Val68, Ala69, Lys70, Ser71, His72, Val73, Tyr80, Asn87, Ala88, Phe89, Tyr90--
GLQVISPPLQR−171.949Asp10, Lys29, Arg39, Glu113Asp10, His13, Trp15, Asp27, Ser28, Lys29, Gly38, Arg39, Ser40, Arg84, Gln92, Tyr110, His111, Leu112, Glu113, His114, His115Asp10Arg84
GHIITVAR−151.394Gly58, Ser71Trp57, Gly58, Ala59, Pro60, Val68, Ala69, Ser71, His72, Asn87, Phe89, Tyr90--
GLIVMAR−144.315Arg84(3)Lys29, Val37, Gly38, Arg39, Gln83, Arg84, Pro85, Asp86, His111, Leu112, Glu113, His114, His115Asp86(3)Glu113
HempQGQIVTVPQNHAVVK−177.911Tyr80, Asn87, Tyr90Trp57, Pro60, Ala69, Lys70, Ser71, His72, Tyr80, Leu82, Asn87, Ala88, Phe89, Tyr90-Ala69
PQFLVGASSILR−175.973Trp57, Ser71, Asn87Trp57, Lys67, Val68, Lys70, Ser71, His72, Tyr80, His81, Leu82, Gln83, Asn87, Ala88, Phe89-His81
LGNLTSYQR−165.451Trp57, Ser71, Tyr80, Arg84, Asn87Trp57, Val68, Ala69, Lys70, Ser71, His72, Tyr80, Leu82, Gln83, Asp86, Asn87, Ala88, Tyr90Asp86(5)-
VQVVNHMGQK−155.517Ala69, Ser71, Ala88Trp57, Ala59, Pro60, Val68, Ala69, Lys70, Asn87, Ala88, Phe89, Tyr90--
ESVILPTSAASPPVK−155.494Lys70, Asn87Pro35, Trp57, Gly58, Ala59, Pro60, Val68, Ala69, Lys70, Ser71, Tyr80, His81, Leu82, Pro85, Asn87, Ala88, Phe89, Tyr90--
SunflowerGHIVNVGQDLQIVR−175.043Ala69, His72, Asn87Trp57, Pro60, Val68, Ala69, Lys70, His72, Tyr80, Asn87, Ala88, Tyr90Lys70-
VIQNLPNQCDLEVQQCTTCTG−172.423Arg48, Ser71(2), Gln75, Arg76Val23, Arg48, Gly49, Glu50, Thr52, Gln54, Ser71, Gln75, Arg76, Ile97, Asn99, Gly100, Asn101Arg48(2)Thr52
PirBvpYNRVGRLKL−174.899Trp57, Val68, Asn87, Ala88Trp57, Val68, Ala69, Lys70, Ser71, His72, Tyr80, Leu82, Pro85, Asn87, Ala88, Phe89--
WADNDSYNNANQD−170.618His72, Asn87, Ala88Pro35, Ile53, Trp57, Ala69, Lys70, Ser71, His72, Tyr80, Leu82, Pro85, Asn87, Ala88, Phe89Lys70Leu82(2)
FVVGENSGKPSVRLQL−167.625Ser28, Gly38, Gln83, His111, Leu112Ser28, Lys29, His30, Thr31, Ile33, Glu36, Val37, Gly38, Arg39, Ser40, His81, Gln83, Arg84, His111, Leu112, Glu113Glu36Ser28, Ile33
YELFHPDEF−159.953Ser71(3), Tyr80Pro35, Trp57, Lys67, Val68, Ala69, Lys70, Ser71, His72, Tyr80, Asn87, Ala88, Phe89, Tyr90-Val68
DEIPQPLKPNM−148.916Tyr80, Asn87Trp57, Gly58, Ala69, Lys70, Tyr80, Pro85, Asn87, Ala88, Phe89, Tyr90--
MLADQEGSDKVAA−142.536Lys70, Ser71, Asn87(3)Trp57, Gly58, Ala59, Pro60, Val68, Ala69, Lys70, Ser71, His72, Val73, Asn87, Phe89, Tyr90-Ser71
* Residues within two PirAvp regions (15-WTVEPNGGVTEVDSKHTPIIPEVGRS-40 and 52-TIQYQWGAPFMAGGWKVAKSHVVQRDET-79) that were previously reported to bind to PirBvp to form a toxic complex are marked in boldface type. Number in brackets indicates the number of interactions. Six PirBvp regions (322-YNRVGRLKL-330, 214-WADNDSYNNANQD-226, 386-FVVGENSGKPSVRLQL-401, 426-YELFHPDEF-434, 290-DEIPQPLKPNM-300, and 409-MLADQEGSDKVAA-421) that were previously reported to bind to PirAvp to form a toxic complex were analyzed for comparison.
Table
Table 5. Docking scores for interactions between six best-scored oilseed peptides and PirBvp, with details on the PirBvp residues involved in different types of interactions.
Table 5. Docking scores for interactions between six best-scored oilseed peptides and PirBvp, with details on the PirBvp residues involved in different types of interactions.
SourcePeptideDocking ScoreInteraction with PirBvp Residues *
Hydrogen BondHydrophobic InteractionExternal Bond
RapeISYVVQGMGISGR−209.710Tyr35, Glu73, Tyr359Tyr35, Ala36, Ala39, Met40, Phe43, Ile48, Pro49, Asn60, Ile61, Pro64, Asp71, Ile72, Glu73, Gln78, Tyr359, Lys361, Gln400, GLn402, Phe429, Pro431, Thr436-
LTFVVHGHALMGK−208.710Gly53, Asn60Tyr35, Ala39, Met40, Tyr50, Ala51, Gly52, Ser56, Thr57, Asn60, Asn114, Glu118, Tyr359, Phe429, Pro431, Phe434, Gly435, Thr436-
QSLGVPPQLGNACNLDNLDVLQPTETIK−201.070Asn28(2), Tyr35, Gly403Leu25, Asp27, Asn28, Tyr29, Tyr35, Thr66, Pro67, Pro70, Ile72, Asp230, Ala273, Val274, Tyr359, Asp364, Val397, Arg398, Gln400, Glu402, Gly403, His404, Phe429, Thr436Gln400(2)
SesameISTINSQTLPILSQLR−198.924Arg21, Tyr29, Tyr359, Ser381, Arg398(2)Arg21, Tyr29, Tyr35, Met40, Phe43, Ile48, Gln84, Asp85, Glu89, Tyr359, Lys361, Tyr362, Ser381, Arg398, Gln400, Glu402, Thr436-
HempPQFLVGASSILR−198.337Tyr359, Lys361, Phe434Tyr35, Ala36, Ala39, Met40, Ile48, Asn60, Ile61, Pro64, Asn65, Thr66, Ile72, Tyr359, Lys361, Tyr362, Ser381, Asp383, Gly403, Phe429, Phe434, Gly435, Thr436Tyr35
VQVVNHMGQK−196.671-Tyr35, Ala39, Met40, Phe43, Ile48, Pro49, Tyr50, Ala51, Thr57, Asn60, Ile61, Asn65, Thr66, Pro67, Tyr359, Lys361, Gly403, His404, Phe429, Pro431, Thr436Asn65
PirAvpTIQYQWGAPFMAGGWKVAKSHVVQRDET−219.311Asn28, Asn65Asp27, Tyr29, Glu30, Val31, Tyr35, Met40, Ile48, Tyr50, Ala51, Thr57, Asn60, Pro64, Asn65, Thr66, Gln78, Asp81, Arg82, Gln84, Asp85, Tyr362, Ser396, Val397, Arg398, Glu402Val397
WTVEPNGGVTEVDSKHTPIIPEVGRS−191.468Tyr35, Ser396Asp27, Asn28, Tyr29, Val31, Tyr35, Ala36, Ile72, Gln78, Asp81, Asp364, Pro395, Ser396, Val397, Arg398, Phe429-
* Residues within six PirBvp regions (214-WADNDSYNNANQD-226, 290-DEIPQPLKPNM-300, 322-YNRVGRLKL-330, 386-FVVGENSGKPSVRLQL-401, 409-MLADQEGSDKVAA-421, and 426-YELFHPDEF-434) that were previously reported to bind to PirAvp to form a toxic complex are marked in boldface type. Number in brackets indicates the number of interactions. Two PirAvp regions (52-TIQYQWGAPFMAGGWKVAKSHVVQRDET-79 and 15-WTVEPNGGVTEVDSKHTPIIPEVGRS-40) that were previously reported to bind to PirBvp to form a toxic complex were analyzed for comparison.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Ong, J.-H.; Wong, W.-L.; Wong, F.-C.; Chai, T.-T. Targeting PirAvp and PirBvp Toxins of Vibrio parahaemolyticus with Oilseed Peptides: An In Silico Approach. Antibiotics 2021, 10, 1211. https://doi.org/10.3390/antibiotics10101211

AMA Style

Ong J-H, Wong W-L, Wong F-C, Chai T-T. Targeting PirAvp and PirBvp Toxins of Vibrio parahaemolyticus with Oilseed Peptides: An In Silico Approach. Antibiotics. 2021; 10(10):1211. https://doi.org/10.3390/antibiotics10101211

Chicago/Turabian Style

Ong, Joe-Hui, Wey-Lim Wong, Fai-Chu Wong, and Tsun-Thai Chai. 2021. "Targeting PirAvp and PirBvp Toxins of Vibrio parahaemolyticus with Oilseed Peptides: An In Silico Approach" Antibiotics 10, no. 10: 1211. https://doi.org/10.3390/antibiotics10101211

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