NCTU_Formosa 2017 had completed a Peptide prediction model which can predict peptides for new function. In the model, they used Scoring Card Method (SCM) for machine learning. This year, NCTU_Formosa 2018 continued to use the same method for predicting antimicrobial peptide, in order to seek more candidates for our project.
Team:NCTU Formosa/Dry Lab/Peptide Prediction
Prediction Process
We started from the Uniprot to pick out the peptides with the key word of antimicrobial. After deleting peptides with similarity over 50% with each other, there were totally 425 sequences left. These data will be taken as positive data for our model.
For negative data, we chose the sequences which are not antimicrobial peptides and the sequences length should be around 30 to 300 amino acids. The amount of these random selection is the same as the amount of positive datasets. Therefore, the ratio between positive and negative would be 1:1.
Next, all the datasets were mixed together, 2/3 for training and 1/3 for testing. After training, we got a threshold and the score of each amino acid sequence. Then we can calculate the score of any peptide through this scoring card.
At last, comparing the score of peptides with threshold, we can easily determine whether the unknown peptide might have the function of inhibiting bacteria’s growth.
Figure 1: Scoring card of antimicrobial peptide.
The more blue the dipeptide the less likely it is antimicrobial, the more red the dipeptide, the more likely it is antimicrobial.
Figure 2: The distribution between positive and negative datasets.
The smaller the overlap area, the more precise the peptide prediction model is.
Figure 3: The receiver operating characteristic curve of the peptide prediction model.
Result
Finally, our peptide pridiction model help us to filter out six potential antimicrobial peptide in the table below. In the future, with this model, we can promote our system to further usage, such as pathogentic disease and so on.
Bacteriocin |
Score |
|---|---|
Leucocyclicin Q |
438.38 |
Enterocin B |
464.36 |
Enterocin 96 |
464.06 |
Lacticin Z |
450.92 |
Bovicin HJ50 |
459.87 |
Durancin TW-49M |
478.97 |
Threshold |
431.65 |
If you are interesting in our model, please click the icon to find out how them work in GitHub!
