Next Article in Journal
The 3Rs Principle in Animal Experimentation: A Legal Review of the State of the Art in Europe and the Case in Italy
Previous Article in Journal
Harmonization of Biobank Education for Biobank Technicians: Identification of Learning Objectives
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

PR-1-Like Protein as a Potential Target for the Identification of Fusarium oxysporum: An In Silico Approach

1
Environmental Biotechnology Laboratory (EBL), Department of Biotechnology, University of the Western Cape, Cape Town 7535, South Africa
2
Bioinformatics Research Group, Department of Biotechnology, University of the Western Cape, Cape Town 7535, South Africa
*
Author to whom correspondence should be addressed.
Submission received: 24 December 2020 / Revised: 12 March 2021 / Accepted: 16 March 2021 / Published: 14 April 2021

Abstract

:
Fusarium oxysporum remains one of the leading causes of economic losses and poor crop yields; its detection is strained due to its presentation in various morphological and physiological forms. This research work sought to identify novel biomarkers for the detection of Fusarium oxysporum using in silico approaches. Experimentally validated anti-Fusarium oxysporum antimicrobial peptides (AMPs) were used to construct a profile against Fusarium oxysporum. The performance and physicochemical parameters of these peptides were predicted. The gene for the Fusarium oxysporum receptor protein PR-1-like Protein, Fpr1, was identified and translated. The resulting protein model from the translation was then validated. The anti-Fusarium oxysporum AMPs and Fusarium oxysporum receptor protein 3-D structures were characterized, and their docking interaction analyses were carried out. The HMMER in silico tool identified novel anti-Fusarium oxysporum antimicrobial peptides with good performance in terms of accuracy, sensitivity, and specificity. These AMPs also displayed good physicochemical properties and bound with greater affinity to Fusarium oxysporum protein receptor PR-1-like Protein. The tendency of these AMPs to precisely detect Fusarium oxysporum PR-1-like Protein, Fpr1, would justify their use for the identification of the fungus. This study would enhance and facilitate the identification of Fusarium oxysporum to reduce problems associated with poor crop yield, economic losses, and decreased nutritional values of plants to keep up with the growing population.

1. Introduction

Fusarium oxysporum is a significant threat to agricultural production. Due to its considerable variation of morphological and physiological makeup resulting from an anamorphic species complex, it tends to escape detection [1]. This fungal pathogen is common globally in soils with the tendency to grow saprophytically or colonize plants. Its economic importance ranges from decreased crop yield, reduced nutritional and market value of farm produce to plants’ reduced resistance to the harsh environmental conditions [2]. The pathogen achieves this by blocking the plant’s water-conducting xylem tissues and subsequently producing germinating spores in the host [3]. The consequence of the aforementioned challenges is a negative effect on storage to an off-season period, causing scarcity for the ever-increasing population [4]. Hence, there is a need to prevent its menace for food abundance and security to meet the demand of our growing population.
The pathogenic strains of Fusarium oxysporum may cause infection such as severe vascular wilts and root rot diseases to not only Phaseolus vulgaris but also other plant hosts such as tomato, banana, cotton and legumes [5]. It is also being reported as an emerging human pathogen for immunocompromised patients [6]. Despite this tendency to infect different plant hosts, isolated Fusarium oxysporum strains only infect very few plant species during inoculation [6]. This inconsistency between field and laboratory conditions limits the fungal pathogen’s handling rate; hence, fast and dependable detection of this microbe is imperative with the end goal of appropriate and timely infection management measures.
Several mechanisms have been described that allow Fusarium oxysporum to recognize and penetrate a host and subdue the innate defenses and nutrient components of the host. The overall effect of the virulence factors produced determines both the infectious potential and severity of the disease caused [7]. Plants have the same innate immunity with the tendency to recognize pathogen-associated molecular patterns (PAMPs) or the presence of pathogenesis-related 1 (PR-1) proteins [8]. Of all the PR proteins, Fpr1 is the most highly expressed class with at least 10% of total Protein in infected hosts [7]. A secreted PR-1-like protein, Fpr1, in Fusarium oxysporum has been functionally characterized and its function is required for full virulence [9]. Despite the widespread use of PR-1-like Protein in Fusarium oxysporum as a virulence factor, its biochemical function and biological importance remain elusive.
Several methods are available to detect Fusarium oxysporum strains that are less accurate due to their similarity to other species. One detection method is based on the characteristic shape and size of macro-and microconidia with the presence and absence of chlamydospores [10]. There are also reports of observation of colony appearances, pigmentation, and growth rates on agar media, giving false positive results [11]. These methods’ reliability is questionable in terms of specificity as they can give false positive results with other Fusarium species [12]. Polymerase chain reaction with restriction fragment length polymorphism (PCR-RFLP) of the intergenic spacer (IGS) region is a gold standard and an emerging technique used for the identification of Fusarium oxysporum and has proven to be highly dependable for the differentiation of strains at the intra-specific level [13]. However, the difficulty of the continuous use of this method lies in the fact that Fusarium oxysporum has a cell wall that impedes its efficient lysis and liberation of DNA. However, this can lead to false negative PCR results, creating demand for more sensitive and accurate diagnostics [14].
Antimicrobial peptides (AMPs) are produced as part of the innate immune response, and they exhibit a broad spectrum of activities against pathogens [13]. They have several compensatory characteristics that make them unique, including cationic charge, hydrophobicity, and diverse structural forms [15]. They establish their activities by binding to membranes and cell walls, causing a non-enzymatic disruption with selectivity owing to the membrane composition of different microbes and cell types [16]. Application of AMPs has been reported recently for the early diagnosis of HIV using bioinformatics [17] with its molecular validation in a lateral flow device for Point-of-Care (POC) detection [18]. Thus, due to the AMP compensatory mechanism of actions highlighted above, the false negative results from existing methods such as PCR may be eliminated. Therefore, this present research work aims to use novel AMPs for the accurate detection of Fusarium oxysporum in Phaseolus vulgaris.

2. Materials and Methods

2.1. Data Collection

The anti-Fusarium oxysporum AMPs were retrieved from antimicrobial peptide databases such as Antimicrobial Peptides Database (APD3) [19,20] and Collection of Antimicrobial Peptides (CAMP) [21]. Thereafter, literature mining was carried out to confirm that all the recovered AMPs were either experimentally validated or predicted. Duplicate experimentally validated AMPs were then removed from the list using the Cluster Database at High Identity with Tolerance (CD-HIT) [22].

2.2. Division of Data into Training and Testing Datasets

The screened, experimentally validated AMPs were randomly partitioned into two subsets: 75% of the data were utilized as the training partition (to build each profile), while the remaining 25% was used for testing (including optimization/calibration of profiles).

2.3. Construction of AMPs Profiles

The Hidden Markov Models (HMMER) algorithm version 3.8 [23] was utilized to build specific pathogen-targeted profiles using the constructed datasets. The task was carried out on the terminal of the Ubuntu operating system version 12.04; (Canonical Ltd., London, UK) with the command line used for building the profile written in the flow chart (i) (Figure 1):
For the initial step, the training datasets of each target class were arranged by utilizing the Clustalo alignment tool [24].

2.4. Independent Profile Testing

The autonomous query of the profiles was performed in a step called “Query profiles”. The testing data were queried against each target profile utilizing the command line as stated above in the flow chart (iii) (Figure 1) with an E-value threshold of 0.05.

2.5. Performance Measurement of Each Profile

The statistical measures were carried out utilizing sensitivity, specificity, accuracy, and Matthews Correlation Coefficient as parameters. The parameters utilized are as described below where TP indicates true positive, TN indicates true negative, FP indicates false positive, and FN indicates false negative:
Percentage sensitivity of the anti-Fusarium oxysporum AMPs against a specific pathogen (testing sets) effectively predicted as anti-Fusarium oxysporum AMPs (positive). The equation of the sensitivity is written below as (1):
Sensitivity = ( T P T P + F N ) × 100
Percentage specificity of the non-anti-Fusarium oxysporum AMPs (negative sets) effectively predicted as non-anti-Fusarium oxysporum AMPs (negative). The equation of the specificity is written below as (2):
Specificity = ( T N T N + F P ) × 100
Percentage accuracy of the effectively predicted peptides (anti-Fusarium oxysporum AMPs and non-anti-Fusarium oxysporum AMPs). The equation of the accuracy is written below as (3):
Accuracy = ( T P + T N T P + F P + T N + F N ) × 100
Matthew’s correlation coefficient (MCC) measures the sensitivity and specificity. MCC = 0 is an indication of absolutely random prediction, while MCC = 1 means perfect prediction. See the Equation (4) as below:
MCC = ( ( T P × T N ) ( F N × F P ) ( T P + F N ) × ( T N + F P ) × ( T P + F P ) × ( T N + F N ) )

2.6. Novel Putative Anti-Fusarium oxysporum AMPs Identification

Query of the proteome sequences were carried out by the respective profiles using the list of all proteome sequences collected from the Ensembl database (http://www.ensembl.org/index.html, accessed on 22 December 2019) [25] and the UniProt database (http://www.uniprot.org/, accessed on 23 December 2019) [26]. An E-value cut-off was set to 0.05 for the discovery of putative anti-Fusarium oxysporum AMPs. The accomplishment of this task was done using “hmmsearch” module of the HMMER software with the command line employed stated in the flow chart above (iv) (Figure 1).
Specific FOTrainings.hmm in the profile, target class query.txt representing the species scanned against the profile and resultfile.txt is the output file acquired after testing the species against the constructed Fusarium oxysporum (FO) profile.

2.7. Identification of Receptors

The gene for the receptor, PR-1-like protein, Fpr1, was identified for Fusarium oxysporum (isolate 4287 PR-1-like protein) and collected from the National Center for Biotechnology Information (NCBI) (https://0-www-ncbi-nlm-nih-gov.brum.beds.ac.uk/, accessed on 26 December 2019) [27], through literature mining. Thereafter, curation was performed to verify that the retrieved Fusarium oxysporum gene was complete. Thereafter, the translate tool of ExPAsy (https://web.expasy.org/translate/, accessed on 27 December 2019) [28] was used to translate the reading frame of the coding portion of the gene into protein. BLAST analysis was then performed using the UniProt interface (https://www.uniprot.org/help/uniprotkb, accessed on 23 January 2020) [26] for further assurance of specificity such that the PR-1-like protein of interest was specific for Fusarium oxysporum.

2.8. Evaluation of the Protein Receptor Model

The quality of the receptor PR-1-like protein, Fpr1, was analyzed using BIOVIA (https://www.3ds.com/products-services/biovia/, accessed on 27 January 2020) [29] online software that incorporates modeler tool (https://salilab.org/modeller, accessed on 27 January 2020) [30] by computing the validation/authentication of its residue regions during binding interaction.

2.9. Physicochemical Properties of the Putative Anti-Fusarium oxysporum AMPs and the Fusarium oxysporum Fpr1 Protein

Physicochemical properties of the anti-Fusarium oxysporum AMPs were calculated using the prediction interface of BACTIBASE (http://bactibase.pfba-lab-tun.org/physicochem, accessed on 31 January 2020) [31,32], DBAASP (https://dbaasp.org/, accessed on 31 January 2020) [33], and APD3 (https://wangapd3.com/main.php, accessed on 28 February 2020) [19,34] and the receptor PR-1-like protein, Fpr1, was carried out using ProtParam tool (http://web.expasy.org/protparam/, accessed on 28 March 2020) from the ExPAsy server [35] using the amino acid sequences of the PR-1-like protein, Fpr1, as input.

2.10. Structure Predictions of the Putative Anti-Fusarium oxysporum AMPs and Fusarium oxysporum Proteins

The I-TASSER (Iterative Threading ASSembly Refinement) server, which is an example of a de novo method of peptide or protein structure prediction, was used to generate the putative anti-Fusarium oxysporum AMPs as well as the Fusarium oxysporum PR-1-like protein, Fpr1, structures [36]. In brief, the prediction was performed by uploading each sequence onto the I-TASSER website [37]. PyMOL (Version 1.3), (Schrödinger, Inc., New York, NY, USA) was then used to visualize the 3-D structures of the AMPs and the protein receptor [38].

2.11. Interaction Analysis of the Putative Anti-Fusarium oxysporum AMPs and Fusarium oxysporum Protein

The PatchDock 1.3 web-server that enables the docking of the protein-small ligand molecule, available at http://bioinfo3d.cs.tau.ac.il/PatchDock/ (accessed on 31 March 2020) was used for the docking of the anti-Fusarium oxysporum AMPs to the Fusarium oxysporum PR-1-like protein, Fpr1 [39]. In brief, the PDB files generated from the I-TASSER for the 3-D structures of the anti-Fusarium oxysporum putative AMPs and the Fusarium oxysporum protein receptor were uploaded onto the PatchDock server. The complex formation with the interaction analysis between the anti-Fusarium oxysporum putative AMPs and the PR-1-like protein receptor was achieved using RasMol 2.7.5 Software (NextMove Software Ltd., Cambridge Science Park, UK) [40]. Subsequently, binding energy scores of the complex formed between the AMPs and the receptor protein were computed using HDock server (http://hdock.phys.hust.edu.cn/, accessed on 3 March 2021) [41].

3. Results

3.1. Data Collection

Experimentally validated anti-Fusarium oxyporum AMPs were collected from different databases—literature mining revealed that CAMP, APD3, DBAASP, and BACTIBASE had 2, 32, and 6 experimentally validated anti-Fusarium oxysporum antimicrobial peptides, respectively. After duplicate removal, a final list of 32 anti-Fusarium oxysporum AMPs was generated.

3.2. Profile Construction

The first step in the profile creation was the random partitioning of the experimentally validated AMPs (Table 1). HMMER was then used to cluster, build and search new AMPs with diagnostic relevance against Fusarium oxysporum.

3.3. Testing and Performance Measurement of the Profile

The profile was tested against a positive dataset which represented about a quarter of the dataset, from which the training dataset used for the construction of the profile was derived. In addition, the trained profile was scanned against a negative control dataset, made up of random fragments of 17236 neuropeptides, which had no recorded anti-Fusarium oxysporum activity (Table 2). The profile discriminated against the negative dataset, with only six of its eight positive datasets being a true positive. Thus, the purpose for dividing the AMPs into training and testing datasets was to ascertain the robustness and discriminatory power of the profile built by HMMER [17]. The performance result of the profile also showed that it was specific, accurate, and sensitive with a significant Matthews correlation coefficient (MCC).

3.4. Proteome Sequence Database Query and Discovery of Anti-Fusarium oxysporum AMPs

Scanning of the profile was carried out to identify novel anti-Fusarium oxysporum AMP sequences that adhered to the 0.05 E-value cut-off. This yielded 12 AMPs across all proteomes scanned that matched the profile (Table 3).

3.5. Receptor Identification

Fusarium oxysporum PR-1-like protein, Fpr1, was used as a receptor to serve as targets for the novel antimicrobial peptides for its detection in a plant host. The PR-1-like protein gene (Fpr1) was identified for Fusarium oxysporum from the National Centre for Bioinformatics Institute (NCBI) database. It was translated using the ExPAsy translate tool using the coding unit of the gene. It was projected that this PR-1-like protein is potentially relevant in detecting Fusarium oxysporum because of its compensatory advantages ranging from production in very high concentration to accessibility for detection [42].

3.6. 3-D Model Structure Validation

BIOVIA, an online tool for verification of structure using modeler Ramachandran plot [29], was used to validate and evaluate the protein model’s quality (Figure 2). The model structure of the PR-1-like protein has 91.8% residues in the most favored region, 7.6% residues in an additional allowed region, 0.2% residues in a generously allowed region, and 0.4% residues in the disallowed regions.

3.7. Physicochemical Analysis of the Anti-Fusarium oxysporum AMPs and Fusarium oxysporum PR-1-Like Protein

Physicochemical features such as molecular weight amino acid composition, hydrophobicity, Boman index, net charge, isoelectric potential, and half-life were used to evaluate the anti-Fusarium oxysporum AMPs. From Table 4, BOMK-1 to 9, including 11 and 12, had glycine as a common amino acid. BOMK-7 and 12 had, in addition to glycine, alanine, and cysteine, respectively. BOMK-10 and 13 had cysteine. All the AMPs had significant hydrophobicity values between 32 and 42, with the lowest value observed for BOMK-9. The hydrophobicity result was the percentage of the total hydrophobic amino acids in the peptides as calculated from APD3 and BACTIBASE. All the AMPs had positive charges with the exception of BOMK-11, 12, and 13, which had negative and neutral charges, respectively. The isoelectric point for the AMPs was between 3.75 and 8.70, while the Boman Index was observed to be between 0.26 and 2.04. Lastly, the AMPs had significant half-lives with the lowest value observed for BOMK-10 (1.3 h).
In Table 5, PR-1-like Protein, Fpr1, had common amino acid alanine, a molecular weight of 95,472.44 Da, significant hydrophobicity, an isoelectric point of 10.00, and a half-life of 1.1 h in mammals with aliphatic and instability indices of 82.72 and 71.04 respectively.

3.8. Structure Prediction and Docking

The structure of the anti-Fusarium oxysporum AMPs was predicted using certain parameters such as C score, TM score, and RSMD as indicators (Table 6). All AMPs had significant C score values, TM scores, and RSMD, where the C score between −5 and 2 indicates structural prediction with high confidence. All the protein and AMP structures were predicted with high confidence because they had existing templates for their database validation. The TM scores for the AMPs and the receptor protein were >0.5, indicating correct topology for the anti-Fusarium oxysporum AMPs and Fusarium oxysporum protein, while the RSMD for these AMPs and receptor protein were between 2 and 4 Å indicating good prediction except for BOMK-4 and BOMK-11, indicating ideal predictions.
Representative output images from the I-TASSER server after predicting the 3-D structures of the anti-Fusarium oxysporum AMPs (ligands) and the protein receptors are indicated in Figure 3.

3.9. Protein-Peptide Interaction between Anti-Fusarium oxysporum and Fusarium oxysporum Fpr1

The docking results of the complex between the putative AMPs and PR-1-like protein, Fpr1, is displayed in Table 7. All the AMPs showed good binding affinity to PR-1-like protein, Fpr1, greater than 8741 [43]. It was observed that BOMK-10, 12, 6, and 8 had the highest binding scores with the lowest observed for BOMK-11 and 9, respectively, using PatchDock. These binding geometry scores are indicators of high detection of the Fusarium oxysporum Frp1. Using binding energy scores from the HDock server, it was observed that all the anti-Fusarium oxysporum displayed high binding energy with BOMK-7 having the highest binding energy, followed by BOMK-3 and -5. The complex formation occurred at the Fusarium oxysporum PR-1-like protein, Fpr1, the most favored regions.
The structural complex of the docking results between the PR-1-like protein, Fpr1, of Fusarium oxysporum, and anti-Fusarium oxysporum AMPs downloaded as PDB files and visualized using RasMol software is shown in Figure 4 in which the blue represents Fpr1 and the anti-Fusarium oxysporum AMPs are shown in red. BOMK-1 and -2 bound at the same position, BOMK-3 and -9 bound at the same position, BOMK-4 and -7 bound at the same position, BOMK-5 and -8 bound at the same position, while BOMK-6 bound at the different orientation of the PR-1-like protein from others.

4. Discussion

4.1. Data Retrieval and Profile Construction of the Anti-Fusarium oxysporum AMPs

Anti-Fusarium oxysporum AMPs that have been experimentally approved were recovered from different databases since they have been demonstrated to possess activities against Fusarium oxysporum by utilizing the agar dilution or broth microdilution techniques with the minimum inhibitory concentration (MIC) assay [44]. The rundown of anti-Fusarium oxysporum AMPs from the databases was recovered after eliminating duplicates to take into account specific species/pathogen profile creation.
The training dataset was made up of ¾ of the retrieved peptides required to prepare the algorithm to test whether the functionally critical amino acid consensus is preserved. After this, multiple alignments were created utilizing HMMER, which keeps the profile from being sensitive to little misalignments and report significant E-values. This allows the tendency to capture sequence diversity since the AMPs were obtained from various life forms [45]. Clusters by HMMER likewise permit a minimum measure of closeness between all peptides.

4.2. Testing of the Profiles

The profile constructed utilizing the training dataset was applied against the held out, positive testing dataset to assess the trained model’s ability to recognize and distinguish this subset of AMPs. Since experimentally confirmed AMPs were utilized, the assumption will be that the profiles developed should have the ability to identify different sequences with similar action and reject those that have no anti-Fusarium oxysporum activity. The utilization of a negative dataset (neuropeptides) was done to affirm whether the prepared profiles would discriminate non-anti-Fusarium oxysporum peptides. The utilization of random sequences as a negative dataset is a regularly utilized method [46].
The assessment of the autonomous profile testing was accomplished utilizing the TP, FP, TN, and FN measures as inputs to the sensitivity, specificity, accuracy, and MCC descriptive statistics. The HMMER E-value cut-off was set to 0.05 to improve the discovery capacity of the profile between the true positive anti-Fusarium oxysporum AMP and false negative anti-Fusarium oxysporum AMPs. The FO (against Fusarium oxysporum) profile had six of its eight positive datasets as true positive, bringing about high sensitivity. The MCC is considered to give the best performance estimation of profiles since it provides the best connection between sensitivity, specificity, and accuracy [47]. The high specificity implies the profile had no comparable capacity with different profiles. The accuracy result was exceptionally high for all the profiles demonstrating the elimination of errors by invalidating misclassified AMPs from both positive and negative datasets. The MCC value “0.5 to 1” relates to correct prediction, while “0” focuses on an irregular forecast. MCC is considered the most robust estimation for assessing the prediction of profile performance. Along these lines, the FO profile shows the right prediction.
HMMER utilized a default E-value of 0.05 for each hit viewed as a true positive. The anti-Fusarium oxysporum profile yielded true positives with E-values lower than 0.05 showing that there was just a 5% possibility that the hit was false or arbitrary. This outcome concurs with the work of Bhadra, Yan [48] where performance was analyzed in terms of accuracy, specificity, sensitivity, and Matthews Correlation Coefficient (MCC) utilizing benchmark datasets as information sources.

4.3. Proteome Sequence Database Query and Discovery of Anti-Fusarium oxysporum AMPs

The discovery stage was to identify novel AMPs with the ability to detect anti-Fusarium oxysporum in infected plant tissues. This was carried out to discover the AMPs with the same signature/motif as the input sequences. A final list of twelve AMPs was identified, and the AMPs were categorized according to their E-values, with those having the smallest E-values considered the most probable putative anti-Fusarium oxysporum AMPs. There was a very small likelihood that these peptides were incorrectly predicted to be anti-Fusarium oxysporum AMPs.

4.4. Receptor Identification

PR-1-like protein, Fpr1, is a protein of Fusarium oxysporum used for proteolytic processing and activation of secreted effectors by fungal and plant host proteases (Avr4). It is a well-characterized type of PR-1 like protein in humans that has been associated with rudimentary biological processes such as cancer, reproduction, and immune response, which are inferred indirectly based on gene expression, localization in specific cell types (glioma or sperm cells), or in response to certain stimuli (pathogen attack) rather than by firm genetic evidence [49]. The highly specific role of PR-1-like protein, Fpr1, during fungus-host interaction makes it a promising target for Fusarium oxysporum detection.
From Table 4, PR-1-like protein, Fpr1, of Fusarium oxysporum is a moderately stable protein. The 3-D model structure validation using BIOVIA also supports its use because of its high quality in terms of the distribution of amino acid residues (Table 1), and thus, this justifies its use for the detection of this fungus [50,51].

4.5. Physicochemical Analysis

The physicochemical properties of the putative AMPs were resolved by utilizing APD and BACTIBASE to guarantee that the distinguished sequences adjust to other AMPs dependent on the qualities estimated. The hydrophobicity result, which was lower than 30% is not an ideal physicochemical parameter [52]. Peptides with higher hydrophobicity would penetrate further into the cell’s hydrophobic center to exert their antimicrobial effects through several mechanisms exhibited by the peptides [53]. All the anti-Fusarium oxysporum AMPs which were positively charged demonstrated congruity of ideal AMPs with improved antimicrobial activity. Notwithstanding, the absence of the positive charge in the net charge of BOMK-10, 11, and 12 does not imply a lack of antimicrobial activity since some negatively charged AMPs have quite recently been accounted for. For example, the surfactant-related anionic peptide in the APD3 database (AP00528) with a net charge of −5 has an anti-bacterial action, and maximin H5 with a charge between −1 and −7 has a bacterial growth restraint action against Listeria monocytogenes [54]. The range of isoelectric values of the AMPs between 3.75 and 8.70 shows characteristic solubility properties for the AMPs in acid and alkaline media despite the variability of charges [55]. The isoelectric point (pI) of peptides is a component of individual amino acids in both original structures. A negative Boman index is said to be related to a more hydrophobic peptide, demonstrating a high protein binding potential, while a more hydrophilic peptide will, in general, have a more positive index [56]. In any case, the propensity of certain peptides to be positive in their Boman index values has been associated with the capacity to identify HIV in a lateral flow device [17].
The physicochemical parameters of the PR-1-like protein, Fpr1 (Table 5) indicate that it is an ideal candidate for the identification of Fusarium oxysporum in terms of stability (as indicated by the instability index), with alanine, valine, isoleucine, and leucine being the most abundant contributors to the aliphatic side chains resulting in an increased thermo-stability (with alanine being the most abundant).

4.6. Structure Prediction and Docking Interaction Analysis of the Putative Anti-Fusarium oxysporum and Fusarium oxysporum PR-1-Like Protein

The structure prediction of the AMPs and the Fusarium oxysporum protein receptor was analyzed in Table 6. The C-score is a certainty score for assessing the nature of anticipated models by I-TASSER. Its assessment depends on the significance of threading template arrangements and the combination parameters of the structure assembly simulations, which are regularly in the scope of −5 to 2. A C-score inside this scope of values connotes a model with high certainty [57]. The prediction of the models of the anti-Fusarium oxysporum AMPs and the Fusarium oxysporum Fpr1 had high confidence in terms of the templates used for their prediction.
On the other hand, TM-score is a proposed scale for estimating the basic convergence/similarity between two structures [58]. A TM-score of >0.5 indicates a correct topology model, and a TM-score of <0.17 means a random similarity. All the AMPs, including the receptor protein, had the correct topology without arbitrary similarity to any other models.
Even though there is certifiably not a characterized RMSD value for 3-D structure prediction, an RMSD estimation of 2–4 Å is viewed as acceptable, and an RMSD of ≤1 Å is considered to be ideal. All models had ideal qualities for RMSD. The results similarly showed that all the AMPs different secondary structures, including α-helices, parallel β-sheet, anti-parallel β-sheet, extended, and loop conformational structures. The outcomes observed associate with the various structural conformations displayed by known AMPs. Examples of known AMPs and their structures include tachyplesin from horseshoe crabs and bovine lactoferricin, which have beta-sheet structures [59]; magainin analog and melittin having alpha-helical conformations [60]. Consequently, the peptides can be considered bona fide AMPs. In any case, the AMPs identified in this study are thought to be putative anti-Fusarium oxysporum peptides because of the absence of experimental proof for these molecules at present.
Utilizing the binding geometry scores in Table 7, all the putative AMPs indicated a huge binding affinity to the PR-1-like protein of Fusarium oxysporum. The AMPs also displayed high binding energy scores with PR-1-like protein, Fpr1, with the AMPs having the most noteworthy inclination to identify the fungus. PatchDock and HDock servers use the scoring function as provided in this research to sample ligands’ conformations on the protein receptor [39,41]. The HDock server, for instance, uses a flexible receptor molecular docking approach to estimate and assess the non-bonded (electrostatic and van der Waals) interactions utilizing the classical force-field-based scoring function [41]. The utilization of HMMER for the discovery of putative AMPs in this research can be used to identify Fusarium oxysporum in plants by utilizing PR-1-like protein, Fpr1, as a target under high sensitivity, specificity, and accuracy.

5. Conclusions

This research identified novel AMPs for the potential diagnosis of Fusarium oxysporum using HMMER in silico technology, where 12 anti-Fusarium oxysporum AMPs were generated. The putative anti-Fusarium oxysporum AMPs showed conformity to other known AMPs in terms of their physicochemical characteristics. This diagnostic system’s primary goal is to ease the search and identify a standard reference for a biomarker for early detection of the fungus to solve the current problem, which leads to the reduction of crop yield, market value, and nutritional value of crop plants, including Phaseolus vulgaris. AMPs have demonstrated incredible promise in evading the downsides related to the current diagnostic systems of this fungus. This research work could be pursued for molecular validation through the binding of these AMPs with the PR-1-like protein, Fpr1, using an “on/off” binding experiment in an LFD setting to develop a prototype with these specific AMPs conjugated to gold nanoparticles (AuNPs) to accurately and sensitively detect the fungal pathogen within plant samples.

Author Contributions

Conceived and designed the experiments: M.K., O.O.B.; analyzed the data: A.G., O.O.B.; and wrote the paper: O.O.B. All authors have read and agreed to the published version of the manuscript.

Funding

National Research Foundation, South Africa (120712). M.K. and the APC is funded by GrainSA (GB0200066).

Data Availability Statement

All data generated or analyzed during this study are included in this published article.

Acknowledgments

The authors are very grateful to the Department of Biotechnology for the provision of an enabling environment to perform this research. Special appreciation to National Research Foundation (NRF), South Africa, and South African National Zakat Fund (SANZAF).

Conflicts of Interest

The authors declare that there is no conflict of interest whatsoever.

References

  1. Michielse, C.B.; Rep, M. Pathogen profile update: Fusarium oxysporum. Mol. Plant Pathol. 2009, 10, 311–324. [Google Scholar] [CrossRef] [PubMed]
  2. Koike, S.T.; Gordon, T.R. Management of Fusarium wilt of strawberry. Crop. Prot. 2015, 73, 67–72. [Google Scholar] [CrossRef]
  3. Ploetz, R.C.; Herbert, J.; Sebasigari, K.; Hernandez, J.H.; Pegg, K.G.; Ventura, J.A.; Mayato, L.S. Importance of Fusarium Wilt in Different Banana Growing Regions; APS Press: St. Paul, MN, USA, 1990; pp. 9–26. [Google Scholar]
  4. Jacobson, D.; Gordon, T. Fusarium oxysporum f. sp. melonis: A case study of diversity within a forma specialis. Phytopathologyogy 1991, 81, 1064–1067. [Google Scholar]
  5. Ploetz, R.C. Fusarium Wilt of Banana. Phytopathologyogy 2015, 105, 1512–1521. [Google Scholar] [CrossRef] [Green Version]
  6. Ordoñez, N.; García-Bastidas, F.; Laghari, H.B.; Akkary, M.Y.; Harfouche, E.N.; Al Awar, B.N.; Kema, G.H.J. First Report of Fusarium oxysporum f. sp. cubense Tropical Race 4 Causing Panama Disease in Cavendish Bananas in Pakistan and Lebanon. Plant Dis. 2016, 100, 209. [Google Scholar] [CrossRef]
  7. Lu, S.; Edwards, M.C. Molecular Characterization and Functional Analysis of PR-1-Like Proteins Identified from the Wheat Head Blight Fungus Fusarium graminearum. Phytopathology 2018, 108, 510–520. [Google Scholar] [CrossRef] [Green Version]
  8. Peng, Y.; Van Wersch, R.; Zhang, Y. Convergent and Divergent Signaling in PAMP-Triggered Immunity and Effector-Triggered Immunity. Mol. Plant Microbe Interact. 2018, 31, 403–409. [Google Scholar] [CrossRef] [Green Version]
  9. Baysal, O.; Karaaslan, C.; Siragusa, M.; Alessandro, R.; Carimi, F.; De Pasquale, F.; Da Silva, J.A.T. Molecular markers reflect differentiation of Fusarium oxysporum forma speciales on tomato and forma on eggplant. Biochem. Syst. Ecol. 2013, 47, 139–147. [Google Scholar] [CrossRef]
  10. Gopinath, V.; Velusamy, P. Extracellular biosynthesis of silver nanoparticles using Bacillus sp. GP-23 and evaluation of their antifungal activity towards Fusarium oxysporum. Spectrochim. Acta Part A Mol. Biomol. Spectrosc. 2013, 106, 170–174. [Google Scholar] [CrossRef]
  11. Henry, P.M.; Kirkpatrick, S.C.; Islas, C.M.; Pastrana, A.M.; Yoshisato, J.A.; Koike, S.T.; Daugovish, O.; Gordon, T.R. The Population of Fusarium oxysporum f. sp. fragariae, Cause of Fusarium Wilt of Strawberry, in California. Plant Dis. 2017, 101, 550–556. [Google Scholar] [CrossRef] [Green Version]
  12. Lievens, B.; Rep, M.; Thomma, B.P. Recent developments in the molecular discrimination of formae speciales of fusarium oxysporum. Pest. Manag. Sci. 2008, 64, 781–788. [Google Scholar] [CrossRef]
  13. García-Bastidas, F.; Ordóñez, N.; Konkol, J.; Al-Qasim, M.; Naser, Z.; Abdelwali, M.; Salem, N.; Waalwijk, C.; Ploetz, R.C.; Kema, G.H.J. First Report of Fusarium oxysporum f. sp. cubense Tropical Race 4 Associated with Panama Disease of Banana outside Southeast Asia. Plant Dis. 2014, 98, 694. [Google Scholar] [CrossRef] [PubMed]
  14. Lin, Y.-H.; Su, C.-C.; Chao, C.-P.; Chen, C.-Y.; Chang, C.-J.; Huang, J.-W.; Chang, P.-F.L. A molecular diagnosis method using real-time PCR for quantification and detection of Fusarium oxysporum f. sp. cubense race 4. Eur. J. Plant Pathol. 2012, 135, 395–405. [Google Scholar] [CrossRef]
  15. Tassanakajon, A.; Rimphanitchayakit, V.; Visetnan, S.; Amparyup, P.; Somboonwiwat, K.; Charoensapsri, W.; Tang, S. Shrimp humoral responses against pathogens: Antimicrobial peptides and melanization. Dev. Comp. Immunol. 2018, 80, 81–93. [Google Scholar] [CrossRef] [PubMed]
  16. Vogel, H.; Müller, A.; Heckel, D.G.; Gutzeit, H.; Vilcinskas, A. Nutritional immunology: Diversification and diet-dependent expression of antimicrobial peptides in the black soldier fly Hermetia illucens. Dev. Comp. Immunol. 2018, 78, 141–148. [Google Scholar] [CrossRef]
  17. Tincho, M.; Gabere, M.N.; Pretorius, A. In Silico Identification and Molecular Validation of Putative Antimicrobial Peptides for HIV Therapy. J. AIDS Clin. Res. 2016, 7. [Google Scholar] [CrossRef] [Green Version]
  18. Williams, M.E.; Tincho, M.; Gabere, M.; Uys, A.; Pretorius, M.M.A.A. Molecular Validation of Putative Antimicrobial Peptides for Improved Human Immunodeficiency Virus Diagnostics via HIV Protein p24. J. AIDS Clin. Res. 2016, 7, 2. [Google Scholar] [CrossRef] [Green Version]
  19. Wang, Z.; Wang, G. APD: The antimicrobial peptide database. Nucleic Acids Res. 2004, 32 (Suppl. 1), D590–D592. [Google Scholar] [CrossRef] [Green Version]
  20. Wang, J.T.L.; Zaki, M.J.; Toivonen, H.; Shasha, D. Data Mining in Bioinformatics; Springer: London, UK, 2005; pp. 3–8. [Google Scholar]
  21. Thomas, S.; Karnik, S.; Barai, R.S.; Jayaraman, V.K.; Idicula-Thomas, S. CAMP: A useful resource for research on antimicrobial peptides. Nucleic Acids Res. 2010, 38 (Suppl. 1), D774–D780. [Google Scholar] [CrossRef] [Green Version]
  22. Li, W.; Godzik, A. Cd-hit: A fast program for clustering and comparing large sets of protein or nucleotide sequences. Bioinformatics 2006, 22, 1658–1659. [Google Scholar] [CrossRef] [Green Version]
  23. Roddy, J. Reducing False Sequence Annotation Due to Alignment Overextension. In Proceedings of the Conference on Undergraduate Research (UMCUR), Missoula, MT, USA, 27 April 2018. [Google Scholar]
  24. Fang, Q. Predicting Functional Alterations Caused By Non-Synonymous Variants in CHO Using Models Based on Phylogenetic Tree and Evolutionary Preservation. Ph.D. Thesis, University of Sheffield, Sheffield, UK, 2018. [Google Scholar]
  25. Hubbard, T.; Barker, D.; Birney, E.; Cameron, G.; Chen, Y.; Clark, L.; Cox, T.; Cuff, J.; Curwen, V.; Down, T.; et al. The Ensembl genome database project. Nucleic Acids Res. 2002, 30, 38–41. [Google Scholar] [CrossRef] [Green Version]
  26. Consortium, U. UniProt: A hub for protein information. Nucleic Acids Res. 2015, 43, D204–D212. [Google Scholar] [CrossRef]
  27. Pruitt, K.D.; Tatusova, T.; Maglott, D.R. NCBI reference sequence (RefSeq): A curated non-redundant sequence database of genomes, transcripts and proteins. Nucleic Acids Res. 2005, 33 (Suppl. 1), D501–D504. [Google Scholar]
  28. Artimo, P.; Jonnalagedda, M.; Arnold, K.; Baratin, D.; Csardi, G.; De Castro, E.; Duvaud, S.; Flegel, V.; Fortier, A.; Gasteiger, E.; et al. ExPASy: SIB bioinformatics resource portal. Nucleic Acids Res. 2012, 40, W597–W603. [Google Scholar] [CrossRef]
  29. Kemmish, H.; Fasnacht, M.; Yan, L. Fully automated antibody structure prediction using BIOVIA tools: Validation study. PLoS ONE 2017, 12, e0177923. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  30. Šali, A.; Potterton, L.; Yuan, F.; van Vlijmen, H.; Karplus, M. Evaluation of comparative protein modeling by MODELLER. Proteins Struct. Funct. Bioinf. 1995, 23, 318–326. [Google Scholar] [CrossRef] [PubMed]
  31. Hammami, R.; Zouhir, A.; Ben Hamida, J.; Fliss, I. BACTIBASE: A new web-accessible database for bacteriocin characterization. BMC Microbiol. 2007, 7, 89. [Google Scholar] [CrossRef] [Green Version]
  32. Hammami, R.; Zouhir, A.; Le Lay, C.; Ben Hamida, J.; Fliss, I. BACTIBASE second release: A database and tool platform for bacteriocin characterization. BMC Microbiol. 2010, 10, 22. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  33. Pirtskhalava, M.; Gabrielian, A.; Cruz, P.; Griggs, H.L.; Squires, R.B.; Hurt, D.E.; Grigolava, M.; Chubinidze, M.; Gogoladze, G.; Vishnepolsky, B.; et al. DBAASP v.2: An enhanced database of structure and antimicrobial/cytotoxic activity of natural and synthetic peptides. Nucleic Acids Res. 2016, 44, D1104–D1112. [Google Scholar] [CrossRef]
  34. Wang, G.; Li, X.; Wang, Z. APD2: The updated antimicrobial peptide database and its application in peptide design. Nucleic Acids Res. 2008, 37 (Suppl. 1), D933–D937. [Google Scholar] [CrossRef] [Green Version]
  35. Kyte, J.; Doolittle, R.F. A simple method for displaying the hydropathic character of a protein. J. Mol. Biol. 1982, 157, 105–132. [Google Scholar] [CrossRef] [Green Version]
  36. Eswar, N.; Sali, A. Protein structure modeling. In From Molecules to Medicines; Springer: Berlin/Heidelberg, Germany, 2009; pp. 139–151. [Google Scholar]
  37. Roy, A.; Kucukural, A.; Zhang, Y. I-TASSER: A unified platform for automated protein structure and function prediction. Nat. Protoc. 2010, 5, 725–738. [Google Scholar] [CrossRef] [Green Version]
  38. De Lano, W.L. PyMOL. Available online: http://www.pymol.org (accessed on 31 March 2020).
  39. Schneidman-Duhovny, D.; Inbar, Y.; Nussinov, R.; Wolfson, H.J. PatchDock and SymmDock: Servers for rigid and symmetric docking. Nucleic Acids Res. 2005, 33 (Suppl. 2), W363–W367. [Google Scholar] [CrossRef] [Green Version]
  40. Sayle, R.A.; Milner-White, E.J. RASMOL: Biomolecular graphics for all. Trends Biochem. Sci. 1995, 20, 374–376. [Google Scholar] [CrossRef]
  41. Yan, Y.; Tao, H.; He, J.; Huang, S.-Y. The HDOCK server for integrated protein–protein docking. Nat. Protoc. 2020, 15, 1829–1852. [Google Scholar] [CrossRef] [PubMed]
  42. Perez-Nadales, E.; Di Pietro, A. The transmembrane protein Sho1 cooperates with the mucin Msb2 to regulate invasive growth and plant infection in Fusarium oxysporum. Mol. Plant Pathol. 2015, 16, 593–603. [Google Scholar] [CrossRef] [PubMed]
  43. Gujjula, K.R. Prediction and Comparison of HIV-1 Protease Inhibitor Binding Energies by Various Molecular Docking Methods. Ph.D. Thesis, National Institute of Technology Rourkela, Orissa, India, 2008. [Google Scholar]
  44. Kim, B.C.; Kim, Y.S.; Chung, J.Y.; Jurng, J.S.; Song, M.Y. Single-Stranded Nucleic Acid Aptamers Specifically Binding to Klebsiella Pneumoniae and Method for Detecting K. pneumonia Using the Same. U.S. Patent No. 9,284,550, 15 March 2016. [Google Scholar]
  45. Madera, M.; Gough, J. A comparison of profile hidden Markov model procedures for remote homology detection. Nucleic Acids Res. 2002, 30, 4321–4328. [Google Scholar] [CrossRef] [PubMed]
  46. Sharma, A.; Kapoor, P.; Gautam, A.; Chaudhary, K.; Kumar, R.; Chauhan, J.S.; Tyagi, A.; Raghava, G.P.S. Computational approach for designing tumor homing peptides. Sci. Rep. 2013, 3, srep01607. [Google Scholar] [CrossRef]
  47. Liu, L.; Oza, S.; Hogan, D.; Perin, J.; Rudan, I.; Lawn, E.J.; Cousens, S.; Mathers, C.; Black, E.R. Global, regional, and national causes of child mortality in 2000–13, with projections to inform post-2015 priorities: An updated systematic analysis. Lancet 2015, 385, 430–440. [Google Scholar] [CrossRef]
  48. Bhadra, P.; Yan, J.; Li, J.; Fong, S.; Siu, S.W.I. AmPEP: Sequence-based prediction of antimicrobial peptides using distribution patterns of amino acid properties and random forest. Sci. Rep. 2018, 8, 1–10. [Google Scholar] [CrossRef] [Green Version]
  49. Silva, C.M.; Carneiro, F.; O’Neill, A.; Cabral, J.S.M.; Guebitz, G.; Cavaco-Paulo, A. Cutinase? A new tool for biomodification of synthetic fibers. J. Polym. Sci. Part A Polym. Chem. 2005, 43, 2448–2450. [Google Scholar] [CrossRef] [Green Version]
  50. Alexandrakis, G.; Jalali, S.; Gloor, P. Diagnosis of Fusarium keratitis in an animal model using the polymerase chain reaction. Br. J. Ophthalmol. 1998, 82, 306–311. [Google Scholar] [CrossRef] [PubMed]
  51. Kim, Y.-H.; Ahn, J.-Y.; Moon, S.-H.; Lee, J. Biodegradation and detoxification of organophosphate insecticide, malathion by Fusarium oxysporum f. sp. pisi cutinase. Chemosphere 2005, 60, 1349–1355. [Google Scholar] [CrossRef]
  52. Chen, Y.; Guarnieri, M.T.; Vasil, A.I.; Vasil, M.L.; Mant, C.T.; Hodges, R.S. Role of Peptide Hydrophobicity in the Mechanism of Action of α-Helical Antimicrobial Peptides. Antimicrob. Agents Chemother. 2006, 51, 1398–1406. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  53. Tachi, T.; Epand, R.F.; Epand, R.M.; Matsuzaki, K. Position-dependent hydrophobicity of the antimicrobial magainin peptide affects the mode of peptide-lipid interactions and selective toxicity. Biochemistry 2002, 41, 10723–10731. [Google Scholar] [CrossRef] [PubMed]
  54. Jodoin, J.; Hincke, M.T. Histone H5 is a potent Antimicrobial Agent and a template for novel Antimicrobial Peptides. Sci. Rep. 2018, 8, 1–15. [Google Scholar] [CrossRef]
  55. Abraham, T.; Prenner, E.J.; Lewis, R.N.; Mant, C.T.; Keller, S.; Hodges, R.S.; McElhaney, R.N. Structure–activity relationships of the antimicrobial peptide gramicidin S and its analogs: Aqueous solubility, self-association, conformation, antimicrobial activity and interaction with model lipid membranes. Biochim. Biophys. Acta BBA Biomembr. 2014, 1838, 1420–1429. [Google Scholar] [CrossRef] [Green Version]
  56. Bakare, O.O.; Keyster, M.; Pretorius, A. Identification of biomarkers for the accurate and sensitive diagnosis of three bacterial pneumonia pathogens using in silico approaches. BMC Mol. Cell Biol. 2020, 21, 1–13. [Google Scholar] [CrossRef] [PubMed]
  57. Zhang, Y. I-TASSER server for protein 3D structure prediction. BMC Bioinform. 2008, 9, 1–8. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  58. Zhang, Y.; Skolnick, J. Scoring function for automated assessment of protein structure template quality. Proteins Struct. Funct. Bioinform. 2004, 57, 702–710. [Google Scholar] [CrossRef]
  59. Bacalum, M.; Radu, M. Cationic antimicrobial peptides cytotoxicity on mammalian cells: An analysis using therapeutic index integrative concept. Int. J. Peptide Res. Ther. 2015, 21, 47–55. [Google Scholar] [CrossRef]
  60. Mihajlovic, M.; Lazaridis, T. Charge distribution and imperfect amphipathicity affect pore formation by antimicrobial peptides. Biochim. Biophys. Acta Biomembr. 2012, 1818, 1274–1283. [Google Scholar] [CrossRef] [PubMed] [Green Version]
Figure 1. Flow chart of the HMMER command lines. The command line (i) in the flow chart above used “Clustalo” module of the HMMER software for the multiple alignment and GCG postscript output for the graphical printing of the AMPs. The command line hmmbuildin (ii) built the aligned sequences in (i) to enhance the construction of the profile by showing common motifs/signatures within the profile. The command line hmmsearchin (iii) evaluated the performance of the resulting constructed profile in (ii) by querying it on independent datasets. The command line (iv) allowed the identification of the anti-Fusarium oxysporum AMPs.
Figure 1. Flow chart of the HMMER command lines. The command line (i) in the flow chart above used “Clustalo” module of the HMMER software for the multiple alignment and GCG postscript output for the graphical printing of the AMPs. The command line hmmbuildin (ii) built the aligned sequences in (i) to enhance the construction of the profile by showing common motifs/signatures within the profile. The command line hmmsearchin (iii) evaluated the performance of the resulting constructed profile in (ii) by querying it on independent datasets. The command line (iv) allowed the identification of the anti-Fusarium oxysporum AMPs.
Biotech 10 00008 g001
Figure 2. BIOVIA result for modeled PR-1-like Protein (Fpr1) of Fusarium oxysporum using the generated model from I-TASSER. Residues in most favored regions (upper first box on the left), residues in additional allowed regions (lower box on the left), and residues in generously allowed regions (upper box on the right).
Figure 2. BIOVIA result for modeled PR-1-like Protein (Fpr1) of Fusarium oxysporum using the generated model from I-TASSER. Residues in most favored regions (upper first box on the left), residues in additional allowed regions (lower box on the left), and residues in generously allowed regions (upper box on the right).
Biotech 10 00008 g002
Figure 3. 3D structures of the anti-Fusarium oxysporum AMPs and Fusarium oxysporum PR-1-like protein, Fpr1, as determined by I-TASSER. 3D structures of anti-Fusarium oxysporum AMPs (al) represented in red and Fusarium oxysporum Fpr1 (m) represented in blue.
Figure 3. 3D structures of the anti-Fusarium oxysporum AMPs and Fusarium oxysporum PR-1-like protein, Fpr1, as determined by I-TASSER. 3D structures of anti-Fusarium oxysporum AMPs (al) represented in red and Fusarium oxysporum Fpr1 (m) represented in blue.
Biotech 10 00008 g003
Figure 4. Docking interaction analysis of the anti-Fusarium oxysporum AMPs with Fusarium oxysporum PR-1 like protein, Fpr1. The blue color represents PR-1 like protein (Fpr1), the red is the AMPs.
Figure 4. Docking interaction analysis of the anti-Fusarium oxysporum AMPs with Fusarium oxysporum PR-1 like protein, Fpr1. The blue color represents PR-1 like protein (Fpr1), the red is the AMPs.
Biotech 10 00008 g004
Table 1. Profile creation by HMMER.
Table 1. Profile creation by HMMER.
ProfilesTraining DatasetsTesting DatasetsTotal
FO24832
Table 2. Independent testing of the profile.
Table 2. Independent testing of the profile.
True PositiveFalse NegativeTrue NegativeFalse Positive
6217,2360
Sensitivity (%)Specificity (%)Accuracy (%)MCC
7510099.990.87
Table 3. Discovery of anti-Fusarium oxysporum AMPs.
Table 3. Discovery of anti-Fusarium oxysporum AMPs.
OrganismNameAMPsNumber of Amino Acid ResiduesBit ScoresE Values
Selaginella moellendorffiiBOMK-1AlaTrpAlaGlyProGlyCysAsnAsnArgLeu----------ValGlyAlaSerGlnHisGlyGlyTyrSerPheAlaTyrGlnGlyGlnThrAlaAlaAlaTyrAsnThrAlaAsnCysArgGlyValAlaHisThrArgPheSerSerLysGlyGluCysLysSerGlySerValGlnAspCysSerGlyPheGlyTrpArgSerIlePheIleGlnCys8035.35 × 10−8
Selaginella moellendorffiiBOMK-2TrpAlaGlyProGlyCysAsnAsnArgLeuGlu----------GlyAlaSerGlnHisGlyGlyTyrSerPheAlaTyrGlnGlyGlnThrAlaAlaAlaTyrAsnThrAlaAsnCysGlnGlyValAlaHisThrArgPheSerArgLysGlyGluCysLysSerGlySerValGlnAspCysSerGlyPheGlyTrpAsnSerPhePheIleGlnCys8032.83.3 × 10−7
Selaginella moellendorffiiBOMK-3ThrTrpAlaGlyProGlyCysAsnAsnArgLeu----------ValGlyAlaSerGlnHisGlyGlyTyrSerPheGlyTyrGlnGlyGlnThrAlaAlaAlaTyrAsnThrAlaAsnCysGlnGlyValAlaHisThrArgPheSerArgLysGlyGluCysLysSerGlySerValGlnAspCysSerGlyPheGlyTrpAsnSerPhePheIleGlnCys8032.73.6 × 10−7
Selaginella moellendorffiiBOMK-4AlaTrpAlaGlyProGlyCysAsnAsnValLeu----------ValArgAlaSerGlnHisGlyGlyTyrSerPheValTyrGlnGlyGlnThrAlaAlaAlaTyrAsnThrAlaAsnCysArgGlyValAlaHisThrArgPheSerArgLysGlyGluCysLysSerGlySerValGlnAspCysSerGlyPheGlyTrpAsnSerPhePheIleGlnCys8031.11.1 × 10−6
Selaginella moellendorffiiBOMK-5ThrTrpAlaGlyProGlyCysAsnAsnGlnArg----------ValGlyAlaSerGlnHisGlyGlyTyrSerPheGlyTyrGlnGlyGlnThrAlaAlaAlaTyrAsnThrAlaAsnCysGlnGlyValAlaGlnThrArgPheSerAlaLysGlyGluCysLysSerGlySerValGlnAspCysSerGlyPheGlyTrpAsnSerPhePheIleGlnCys8027.61.4 × 10−5
Selaginella moellendorffiiBOMK-6TrpAlaGlyProGlyCysAsnAsnTrpLeuGlu----------AlaSerGlnHisGlyGlyTyrSerValAlaTyrLeuGlyHisAlaAlaAlaAlaTyrAsnThrAlaAsnCysGlnGlyValAlaGlnArgTrpPheArgArgLysGlyHisCysSerSerGlyCysAlaSerGluCysGluGlyPheArgTrpAsnSerIlePheIleGlnCysSerSer8026.43.6 × 10−5
Selaginella moellendorffiiBOMK-7TrpAlaGlyProGlyGlyAsnAsnArgLeuGlu----------AlaSerGlnHisGlyGlyTyrSerValValTyrLeuGlyHisAlaAlaAlaAlaTyrAsnThrAlaAsnCysGlnGlyValAlaGlnArgTrpPheArgArgLysGlyHisCysSerSerGlyCysAlaSerGluCysGluGlyPheArgTrpAsnSerIlePheIleGlnCysSerSer8025.00.00011
Setaria italicBOMK-8ThrSerTrpAlaGlyProGlyCysSerGlyGln----------AsnLeuGlnPheTyrAspGlyGlnGluLysSerTyrGlnGlyGlnThrAlaArgLeuTyrThrGluThrGlyCysAlaGlyThrSerTyrLeuValPheGluAspThrGlnAlaCysGlySerGlyCysAlaSerGluCysGluAspPheGlyTrpArgSerIle7521.80.00073
Oryza sativumBOMK-9LysIleGlnValGluAlaLysSerCysCysProGly----------TyrAsnSerCysArgPheAlaGlyGlySerArgAspThrCysAlaLysLeuSerGlyCysLysIleValCysAspGlyAsnCysLysProProTyr5423.50.00079
Zea maysBOMK-10GlyGlyHisProAspGlyAlaIleProCysGlyGlu----------ValPheGlyCysArgGlyTrpGlyTyrCysGlu3319.80.0037
Solanum lycopersicumBOMK-11AlaGlnGlnCysGlyIleGlnAlaGlyGlyAla----------PheGlyTyrCysGlyThrThrAlaThrAlaTyrCysGlyProGlyCysGlnSerGlnCys4116.00.026
Arabidopsis thalianaBOMK-12ValGlnGluTyrGlyCysProAsnCysLysArg----------GlyGluLeuValMetGluCysAsnLys3017.70.034
BOMK1–12: Anti-Fusarium oxysporum AMPs, “-” means specific amino acid residues which will be made available on request.
Table 4. Physicochemical properties of the anti-Fusarium oxysporum AMPs.
Table 4. Physicochemical properties of the anti-Fusarium oxysporum AMPs.
AMPMol. Mass (Da)Common Amino AcidsHydrophobicity (%)Isoelectric PointBoman Index (Kcal/mol)ChargeHalf-Life (h)
BOMK-18651.30G358.701.89+54.4
BOMK-28531.35G358.281.69+32.8
BOMK-38618.44G338.281.71+37.2
BOMK-48700.36G378.491.70+44.4
BOMK-58509.84G348.281.51+37.2
BOMK-68900.35GA408.401.63+42.8
BOMK-78852.32G378.701.81+52.8
BOMK-88121.76G323.881.42−57.2
BOMK-95796.17C358.701.72+51.3
BOMK-103376.98G424.420.26−230
BOMK-114042.34GC413.750.45−14.4
BOMK-123537.55C407.082.040100
BOMK1–12: Anti-Fusarium oxysporum AMPs.
Table 5. Physicochemical of the PR-1-like protein, Fpr1, of Fusarium oxysporum.
Table 5. Physicochemical of the PR-1-like protein, Fpr1, of Fusarium oxysporum.
ReceptorM. wt (Da)Common Amino AcidHydrophobicity
(%)
Isoelectric PointInstability
Index
Aliphatic
Index
Half-Life
(Hours)
PR-1-like protein95,472.44SLP4010.0071.0482.721.1
Table 6. Structure prediction of the anti-Fusarium oxysporum AMPs and PR-1-like protein, Fpr1, from I-TASSER.
Table 6. Structure prediction of the anti-Fusarium oxysporum AMPs and PR-1-like protein, Fpr1, from I-TASSER.
AMPsC ScoreTM ScoreRSMD (Å)
BOMK-10.840.83 ± 0.082.0 ± 1.6
BOMK-20.490.78 ± 0.102.6 ± 1.9
BOMK-30.830.83 ± 0.082.0 ± 1.6
BOMK-40.880.83 ± 0.081.9 ± 1.6
BOMK-50.740.81 ± 0.092.2 ± 1.7
BOMK-60.770.82 ± 0.092.1 ± 1.7
BOMK-70.730.81 ± 0.092.3 ± 1.8
BOMK-80.670.80 ± 0.092.2 ± 1.7
BOMK-90.420.66 ± 0.133.6 ± 2.5
BOMK-100.860.61 ± 0.143.5 ± 2.4
BOMK-110.580.79 ± 0.091.4 ± 1.3
BOMK-12−0.750.62 ± 0.143.3 ± 2.3
PR-1-like protein−1.440.54 ± 0.159.1 ± 4.6
Table 7. PATCHDOCK results for each AMPs and PR-1-like protein with the geometry binding affinity.
Table 7. PATCHDOCK results for each AMPs and PR-1-like protein with the geometry binding affinity.
AMPsPathdock Geometry Binding Affinity ScoresHDock Binding Energy Scores (Kcal/mol)
BOMK−112,828−254.09
BOMK−212,976−238.30
BOMK−313,776−263.43
BOMK−412,652−244.88
BOMK−513,688−263.63
BOMK−614,334−242.95
BOMK−712,776−273.34
BOMK−814,016−213.09
BOMK−911,468−239.40
BOMK−1014,806−230.06
BOMK−119958−220.11
BOMK−1214,420−236.68
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Bakare, O.O.; Gokul, A.; Keyster, M. PR-1-Like Protein as a Potential Target for the Identification of Fusarium oxysporum: An In Silico Approach. BioTech 2021, 10, 8. https://0-doi-org.brum.beds.ac.uk/10.3390/biotech10020008

AMA Style

Bakare OO, Gokul A, Keyster M. PR-1-Like Protein as a Potential Target for the Identification of Fusarium oxysporum: An In Silico Approach. BioTech. 2021; 10(2):8. https://0-doi-org.brum.beds.ac.uk/10.3390/biotech10020008

Chicago/Turabian Style

Bakare, Olalekan Olanrewaju, Arun Gokul, and Marshall Keyster. 2021. "PR-1-Like Protein as a Potential Target for the Identification of Fusarium oxysporum: An In Silico Approach" BioTech 10, no. 2: 8. https://0-doi-org.brum.beds.ac.uk/10.3390/biotech10020008

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