Title: | Prediction of Antimicrobial Peptides |
---|---|
Description: | Predicts antimicrobial peptides using random forests trained on the n-gram encoded peptides (10.3390/ijms21124310). The implemented algorithm can be accessed from both the command line and shiny-based GUI. The AmpGram model is too large for CRAN and it has to be downloaded separately from the repository: <https://github.com/michbur/AmpGramModel>. |
Authors: | Michal Burdukiewicz [cre, aut] , Katarzyna Sidorczuk [ctb] , Filip Pietluch [ctb] , Dominik Rafacz [aut] , Stefan Roediger [ctb] , Jaroslaw Chilimoniuk [ctb] |
Maintainer: | Michal Burdukiewicz <[email protected]> |
License: | GPL-3 |
Version: | 1.1 |
Built: | 2024-12-19 04:38:24 UTC |
Source: | https://github.com/michbur/ampgram |
Antimicrobial peptides (AMPs) are ancient and evolutionarily conserved molecules widespread in all living organisms that participate in host defence and/or microbial competition. Due to their positive charge, hydrophobicity and amphipathicity, they preferentially disrupt negatively-charged bacterial membranes. AMPs are considered an important alternative to traditional antibiotics, especially in times when the latter are drastically losing their effectiveness. Therefore, efficient computational tools for AMP prediction are essential to identify the best AMP candidates without undertaking expensive experimental studies. AmpGram is our novel tool for predicting AMPs based on the stacked random forests and n-gram analysis, able to successfully predict antimicrobial peptides in proteomes.
AmpGram is available as R function (predict.ampgram_model
) or
shiny GUI (AmpGram_gui
).
AmpGram requires the external package, AmpGramModel, which
contains models necessary to perform the prediction. The model
can be installed using install_AmpGramModel
Maintainer: Michal Burdukiewicz <[email protected]>
Burdukiewicz M, Sidorczuk K, Rafacz D, Pietluch F, Chilimoniuk J, Roediger S, Gagat P. (2020) AmpGram: a proteome screening tool for prediction and design of antimicrobial peptides. (submitted)
Launches graphical user interface that predicts presence of antimicrobial peptides.
AmpGram_gui()
AmpGram_gui()
No return value, called for side effects.
Any ad-blocking software may cause malfunctions.
Predictions made with the AmpGram methods.
A list of length three: random forest for 10-mer analysis, random forest for predictions of AMPs, and a vector of important n-grams.
Function gets sequences recognized as antimicrobial peptides and returns as data.frame.
get_AMPs(x)
get_AMPs(x)
x |
AmpGram predictions for a single protein |
a data.frame with sequences recognized as antimicrobial peptides (AMPs). It consists of two columns:
amino acid sequence of a 10-mer (subsequence of an analyzed peptide) predicted as AMP.
Probability with which a 10-mer is recognized as AMP.
data(AmpGram_predictions) get_AMPs(AmpGram_predictions[[2]])
data(AmpGram_predictions) get_AMPs(AmpGram_predictions[[2]])
Installs AmpGramModel package containing model required for prediction of antimicrobial peptides. Due to large size of our model and file size limit on CRAN, it needs to be stored in the external repository. See readme for more information or in case of installation problems.
install_AmpGramModel()
install_AmpGramModel()
Convert predictions to data.frame Return predictions as data.frame
pred2df(x)
pred2df(x)
x |
results of prediction as produced by |
a data.frame with two columns and number of rows corresponding to the number of peptides/proteins in the results of prediction. Columns contain following information:
Name of an analyzed sequence
Probability that a protein/peptide possesses antimicrobial activity. It assumes values from 0 (non-AMP) to 1 (AMP).
Row names contain sequence name and decision if a peptide/protein is classified
as AMP (TRUE
) or non-AMP (FALSE
).
data(AmpGram_predictions) pred2df(AmpGram_predictions)
data(AmpGram_predictions) pred2df(AmpGram_predictions)
Recognizes antimicrobial peptides using the AmpGram algorithm.
## S3 method for class 'ampgram_model' predict(object, newdata, ...)
## S3 method for class 'ampgram_model' predict(object, newdata, ...)
object |
|
newdata |
|
... |
further arguments passed to or from other methods. |
AmpGram requires the external package, AmpGramModel, which
contains models necessary to perform the prediction. The model
can be installed using install_AmpGramModel
.
Predictions for each protein are stored in objects of class
single_ampgram_pred
. It consists of three elements:
Character vector of amino acid sequence of an analyzed peptide/protein
Numeric vector of predictions for each 10-mer (subsequence of 10 amino acids) of a sequence. Prediction value indicates probability that a 10-mer possesses antimicrobial activity and ranges from 0 (non-AMP) to 1 (AMP).
Named numeric vector of a single prediction value for
a whole peptide/protein. Its value corresponds to the probability that a
peptide/protein exhibits antimicrobial activity. It assumes name TRUE
if probability is equal or greater than 0.5, i.e. peptide/protein is classified
as antimicrobial (AMP), and FALSE
if probability is less that 0.5,
i.e. peptide/protein is classified as non-antimicrobial (non-AMP).
list
of objects of class single_ampgram_pred
. Each object
of this class contains analyzed sequence, values of predictions for 10-mers and
result of the prediction for the whole peptide/protein.
Read sequence data saved in text file.
read_txt(connection)
read_txt(connection)
connection |
a |
The input file should contain one or more amino acid sequences separated by empty line(s).
a list of sequences.
(sequences <- read_txt(system.file("AmpGram/prots.txt", package = "AmpGram")))
(sequences <- read_txt(system.file("AmpGram/prots.txt", package = "AmpGram")))