BACKGROUND
insecticides are used to control the insects affecting the agricultural crops, parasitizing livestock, as well as to eradicate the pests transmitting dangerous infectious diseases.
however, frequent application of insecticides has resulted in the resurgence of pests and appearance of resistant pest species.
insecticide resistance is the heritable change in the sensitivity of a pest population that is reflected in the repeated failure of a product  to achieve the expected level of control when used according to the level of recommendation for that pest species  <cit> .
several studies have indicated the involvement of multiple genes in conferring the resistance to many insect species .
thus, characterization of these genes is useful to understand the development of resistance and designing new strategies to minimize the development of insecticide resistance  <cit> .
three major mechanisms are involved in insecticide resistance  <cit> :  detoxification of insecticides through alteration in the activities of enzymes like esterase, oxidases or glutathione s-transferases  that prevents the insecticide from reaching to its site of action ,  insensitivity of the insecticide target proteins  <cit> , because of which the insecticide no longer binds to its target  <cit>  and  reduction in insecticide uptake due to decrease in permeability of insect cuticle  <cit> .
though there is evidence of alteration in cuticular penetration, most of the studies have focused and evaluated the target site insensitivity and detoxification of insecticides  mechanisms.
moreover, these two mechanisms have been reported to cover a wide range of resistance levels to almost all available insecticides  <cit> .
the cytochrome p <dig> family of genes in insect play an important role in the detoxification of insecticides resulted in the development of resistance to insecticides  <cit> .
besides, gsts have also been reported to be involved in the detoxification of insecticides  <cit> .
as far as target-based mechanism is concerned, there are three main targets for conventional insecticides viz., gaba -gated chloride ion channel, voltage-gated ion channel and acetylcholinesterases  <cit> .
the gaba receptor is the site of target for cyclodiene  insecticides  <cit> , where the resistance to dieldrin  is conferred by the change of a single amino acid in gaba-gated chloride ion channel encoded by rdl gene  <cit> .
further, knockdown resistance  is one of the major forms of resistance to ddt and pyrethroid insecticides  <cit> , which is associated with mutations in the voltage-gated sodium channel .
acetylcholinesterase  in nerve synapses is the target protein for the insecticides like organophosphorus  and carbamate   <cit> .
the point mutation in the insecticide-binding site of ache has been identified as the cause of insensitivity to these insecticides  <cit> .
the above mentioned works help enable to understand the molecular mechanisms involved in the insecticide resistance.
further, the analysis of bio-molecules involved in this phenomenon has confirmed the importance of single genes in target site resistance and involvement of multi-gene families like cytochrome p <dig> in metabolic resistance  <cit> .
several studies on the effects of mutational changes in target proteins on insecticide resistance aid to the knowledge on the insect proteins involved in this process.
for instance, riveron et al.
<cit>  demonstrated that the single amino acid change  in an up regulated gst gene, gste <dig>  confers high level of metabolic resistance to ddt in the malaria vector anopheles funestus.
in another study, nwane et al.
<cit>  identified that two mutations at position  <dig> of the s <dig> transmembrane segment of domain ii in the voltage gated sodium channel i.e., leucine to a phenylalanine  or to a serine  confers resistance to ddt and pyrethroid insecticides in anopheles gambiae.
in the recent past, several studies have identified species-specific insecticide resistant genes through transcriptome and expression profile analysis.
hsu et al.
<cit>  identified  <dig> p <dig>   <dig> gst,  <dig> coe-related genes in bactrocera doralis, representing three major enzyme families involved in insecticide metabolism and resistance.
in another study,  <dig> p <dig>   <dig> gst and  <dig> ces-specific genes of liposcelis bostrychophila were reported to be involved in insecticide resistance, through transcriptome and differential gene expression analysis  <cit> .
recently, cui et al.
<cit>  identified relevant genes in response to flubendiamide insecticide in asian corn borer , through de novo transcriptome and expression-profile analysis.
though the transcriptome and expression profile analysis is one way of identifying the resistance genes, it is species specific.
moreover the expression profile analysis is expensive as well as time consuming.
thus, development of a computational tool for identifying the resistant genes independent of the species and economically as well would help in augmenting the research related to the identification of insecticide resistant genes.
however, no computational tool is reported till date for the discrimination of insecticide resistant proteins from the proteins that do not confer resistance.
keeping this in view, we propose a computational approach to discriminate the insecticide resistant proteins from non-resistant proteins.
the developed computational approach can be used for identification of the resistant proteins across species as well as with minimum resource .
we have also developed an online prediction server that can be easily used by experimental scientist and researchers to predict an unknown protein sequence as either insecticide-resistant or non-resistant protein.
moreover, computational identification of insecticide resistant proteins will supplement the efforts needed to develop insecticides in targeting the resistance proteins.
methods
collection and processing of data
in this study, protein sequences corresponding to four important groups of insecticide resistant genes viz., cytochrome p <dig>  kdr, rdl and ache were collected from insecticide resistance gene database .
we considered these four categories of genes because they represent important families of insecticide resistant genes which are resistant to commonly used insecticides.
besides, the resistant protein sequences were reported to be involved in two important resistance mechanisms viz., detoxification-based and target-based.
further, target-based resistant proteins are confined to three main targets of insecticides i.e., ache, gaba-gated chloride ion channel and voltage-gated sodium channel.
a total of  <dig> sequences  belonging to  <dig> insect species  were collected.
initially, we removed the sequences having non-standard residues.
then, four positive sets having  <dig>   <dig>   <dig> and  <dig> sequences were prepared, where the maximum pair-wise sequence identities were 40%, 60%, 70% and 90% respectively.
the sequences with more than considered level of pair-wise sequence identity were removed using cdhit  <cit> .
for negative set, protein sequences  of the considered species were collected from the uniprot  database.
for the species acyrthosiphon pisum and tribolium castaneum, only the reviewed sequences were collected, as large number of sequence are present in uniprot for these two species.
on the other hand, all the sequences available for remaining nine species were collected.
after removing the sequences having non-standard residues as well as the identical sequences, a total of  <dig> sequences were obtained.
further, to avoid homologous bias in the negative dataset, sequences with >40% pair-wise identity were removed using cdhit.
finally a dataset with  <dig> sequences was obtained and considered as the negative dataset.
feature generation
protein sequences are the strings of amino acid residues, and hence they need to be mapped onto numeric feature vectors before being used as input in machine learning classifier.
in this study, amino acid composition , di-peptide composition , pseudo amino acid composition , composition-transition-distribution  and auto correlation function  were used to transform the protein sequences into numeric feature vectors.
amino acid composition
aac is a basic feature of protein sequence  <cit> , which is closely associated with its attributes, such as sub-cellular location  <cit> , secondary structure content  <cit>  and domain  <cit> .
aac consists of  <dig> discrete numbers, each of which represents the frequency of the native amino acids in a protein sequence.
based on the aac, each protein sequence was encoded into a 20-dimensional numerical vector.
di-peptide composition
one of the limitations of aac is that it does not take into account the local order information of amino acids in the protein.
on the other hand, dpc, which gives a fixed pattern length of  <dig> , encapsulates the global information about each protein sequence and the order it contains  <cit> .
for any di-peptide, its composition was computed as the ratio of the frequency of that di-peptide to the total number possible di-peptide in the protein sequence.
pseudo amino acid composition
the concept of paac was originally introduced by chou  <cit>  for predicting the protein sub-cellular locations and membrane protein types.
based on the conventional aac, chou proposed a set of discrete numbers to take into account the sequence order effects.
paac has been proven to be an extremely effective feature in many proteins and protein-related systems  <cit> .
the paac for each protein sequence can be represented by a -dimensional vector for d-tier correlation factor.
here, the paac was extracted for 1st-tier correlation only, by which each sequence was transformed into a 21-dimensional numeric vector.
for further details, one can refer to  <cit> .
composition-transition-distribution
the ctd feature was introduced by dubchak et al.
<cit>  for predicting protein folding classes.
thereafter, the ctd feature has been adopted by many researchers for protein function and structure studies  <cit> .
in ctd feature, composition  is the number of amino acids of a particular type divided by the total number of amino acids.
transition  characterizes the frequency percentage with which amino acids of a particular type is followed by other amino acids.
distribution  measures the chain length within which the first 25%, 50%, 75% and 100% of the amino acids of a particular type is located respectively.
based on the ctd feature, each protein sequence of length l was encoded to a l+{l*/2} + -dimension numeric vector.
auto correlation function
sequence autocorrelation-based features assume that the disturbances in each area are systematically related to those in adjacent areas  <cit> .
this concept helps to analyze the dependency among the features of sequences in each location.
autocorrelation features were computed based on the distribution of amino acid properties along the sequence, using all the  <dig> amino acid indices available in aaindex database  <cit> .
in this feature encoding, for an autocorrelation of order n, each sequence was transformed into a numeric vector of length 531*n.
supervised learning technique
for classification purpose we used the support vector machine , which is a nonparametric algorithm developed by vapnik  <cit> .
it is a very promising and popular method for pattern recognition that has been widely used for prediction purpose in the field of bioinformatics .
it is proven to be very efficient in many biological analyses due to their ability to handle noise and large input dataset  <cit> .
a brief description about the working principle of svm is described as follows:
consider a binary classification problem with n samples or input vectors x
i ∈ r
d, , where x
i with class levels y
i ∈ {− <dig>  1} can be considered as the i
th protein or vector defined in d-dimensional space .
in present work,  <dig> refers to resistant class and − <dig> represents non-resistant class.
the objective here is to construct a binary classifier from the available sample  that has less probability of misclassifying future sample .
non-linear svm maps input vectors x
i ' s into high dimensional feature space and constructs an optimal separating hyper-plane  that maximizes the distance between hyper-plane and nearest data points of each class in the space.
mathematically, the hyper-plane is represented as y = sgn, where w represents a weight vector that can map training data in the input space to the outer space and b represents bias.
for a two class problem, it can be formulated as wtxi+b≥1ifyi=1wtxi+b≤−1ifyi=− <dig>
the svm training procedure involves optimization of convex quadratic problem i.e., with lagrangian multipliers α
i ≥  <dig>  maximize∑i=1nαi−12∑i=1n∑j=1nαiαjyiyjkxixj subject to the constraints 0 ≤ α
i ≤ c  and ∑i=1nαiyi= <dig>  where c is the regularization parameter that controls trade-off between margin and classification error.
the x
j ' s are called support vectors only if corresponding α
j >  <dig>  after the svm has been trained, the decision function for classification of query sequence  can be formulated as
fx=sgn∑i=1nyiαikx.xi+b.
the choice of the proper kernel function k is important to train svm model because the power of svm comes from the kernel representation that allows the nonlinear mapping of input space to a higher dimensional feature space.
in this work, four commonly used kernel functions  <cit>  viz., linear , polynomial d), radial basis  and sigmoid ) were used, where r, d, γ > <dig> are the kernel parameters.
validation of the model
cross-validation procedure has been widely accepted for assessing the performance of classifiers  <cit> .
thus, we used the 10-fold cross-validation to assess the performance of our approach.
it was carried out by partitioning the dataset into  <dig> approximately equal-sized sets at random, where nine partitions were used to train the model and the remaining one part was used to assess the model accuracy.
this process was repeated  <dig> times in such a way that each partition was tested once in the model.
performance evaluation
different performance metrics viz., sensitivity , specificity , accuracy , precision  and matthew’s correlation coefficient  were used to measure the accuracy of the developed prediction approach.
the sn, sp, ac, pre and mcc parameters are defined as:sn = tp/, sp = tn/, ac = /, pre = tp/, mcc=tp×tn−fp×fn/tp+fn×tp+fp×tn+fn×tn+fp.
true positive  is the number of resistant proteins correctly predicted as resistant proteins, true negative  is the number of non-resistant proteins correctly predicted as non-resistant proteins, false negative  is the number of resistant proteins incorrectly predicted as non-resistant proteins and false positive  is the number of non- resistant proteins incorrectly predicted as resistant proteins.
besides the above mentioned performance metrics, area under receiving operating characteristic curve   <cit>  was also used to measure the predictive ability.
for given false positive rate  and true positive rate  at different threshold values, the auc-roc was computed as ∑i1−βi.Δα+1/2Δ1−β.Δα, where Δ =  − , Δα = α
i − α
i −  <dig> and i =  <dig> , …, m   <cit> .
a subroutine in r programming language was written to compute the values of these performance metrics.
training and testing datasets
using four positive sets and one negative set , four datasets were prepared that consists of both positive and negative sequences.
here each dataset contains a different positive set and the same negative set .
all the four datasets are highly unbalanced as the number sequences present in one class  is much larger than the other class .
to avoid biasness towards the non-resistant class  while predicting using machine learning classifier like svm, balanced datasets were prepared that consists of same number of sequences from both the classes, where the sequences of the major class were randomly drawn from the available sequences of the major class.
for instance, first balanced dataset contains  <dig> positive and  <dig> negative sequences, where the  <dig> negative sequences were randomly drawn from  <dig> negative sequences.
as the generalized predictive ability cannot be assured based on a single dataset,  <dig> sample sets were prepared, where each sample set consists of same number of positive and negative instances.
further, in each sample set, a 10-fold cross validation procedure was adopted.
the performance metrics were computed by taking average over the  <dig> folds as well as over  <dig> sample sets.
mechanism-based classification
the insecticide resistance mechanism can be broadly categorized into two types, viz., target-based mechanism and detoxification-based mechanism.
the rdl, kdr and ache genes come under target-based and cytochrome p <dig> genes come under detoxification-based mechanism.
to test whether the genes under these two categories are different or not, a binary classification was carried out by employing svm, where  <dig> sequences  from target-based and  <dig> sequences  from detoxification-based category were used.
similar to the classification of resistant and non-resistant proteins,  <dig> sample sets were prepared where each sample set consists of  <dig> sequences from each class.
since, there are  <dig> sequences in the detoxification-based category,  <dig> sequences were randomly drawn each time.
as the number of sequences in each sample is not large, leave-one-out cross validation  technique was adopted for classification of detoxification- and target-based resistant proteins.
here, detoxification-based category was considered as positive class and target-based category as negative class.
comparison with blast algorithm
performance of the proposed approach was also compared with that of blastp  <cit> , psi-blast  <cit>  and delta-blast  <cit> , which are powerful algorithms to detect protein homologs.
further, comparison was made through 10-fold cross validation technique.
for cross validation, offline  blast software was used with blastp, psiblast and deltablast modules/programs in which the training set for each fold of cross validation was defined as the database and sequences of the corresponding test set were used as query.
each query sequence was predicted as the resistant or non-resistant category based on the top hit found in the blast search.
three different e-values i.e.,  <dig> ,  <dig> and  <dig> were used to assess the performance of the blastp, psi-blast and delta-blast.
furthermore, performance of the proposed approach was compared based on best feature set with which higher accuracies were obtained as compared to the other feature sets.
performance evaluation using independent dataset
to assess the generalized predictive ability of the proposed approach, its performance was further tested using an independent test dataset.
the independent dataset was collected based on published literature that includes  <dig> cytochrome p <dig>   <dig> kdr,  <dig> rdl and  <dig> ache proteins.
specifically,  <dig> cytochrome p <dig> genes were reported by hsu et al.
<cit> .
out of  <dig>  we used  <dig> as they are available in ncbi.
similarly,  <dig> kdr,  <dig> rdl and  <dig> ache genes were collected from ncbi, based on the study of zuo et al.
<cit> , wondji et al.
<cit>  and li and han  <cit>  respectively.
sequences of the independent test set are provided in additional file  <dig>
development of prediction server
an online prediction server was developed using html and php, where the combination of best feature set and classifier was used.
a developed r-code was executed in background upon submission of the sequences in fasta format to the server.
the user has to submit the protein sequences having only standard amino acid residues.
this server can be used to predict the likelihood of any unknown protein sequence being an insecticide resistant protein with certain probability.
RESULTS
analysis of amino acid compositions
the composition of amino acids in four different groups of insecticide resistant proteins is shown in fig.
<dig>  it is seen that the proportions of leucine  are higher, whereas the proportions of cystene  and tryptophan  are lower in all the four categories.fig.
<dig> composition of amino acids in all the four categories of insecticide resistant proteins.
it is observed that proportions of leucine are higher, whereas proportions of cystene and tryptophan are lower in all the four categories
analysis of kernel functions
based on a sample dataset consisting of  <dig> positive and  <dig> negative sequences that were drawn randomly from the available positive and negative sequences, performance of svm was analyzed.
roc curves for all the four kernels as well as for all the five feature sets are shown in fig.
2a and the corresponding auc-roc values are shown in bar plots .
from the roc curves it is not clear that which kernel is better, whereas from auc-roc plots it is clear that the values of auc-roc are higher for the rbf kernel, irrespective of the feature set used.
though in rbf kernel the auc-roc for acf feature set is highest, it is difficult to choose the best feature set while other three kernels are taken into account.
therefore, all the feature sets and the rbf kernel were used for further analysis.fig.
<dig>
a roc curves of svm for different kernels and features, b bar plots of corresponding auc-roc values.
it is seen that the auc-roc values are higher for rbf kernel as compared to other kernels
cross-validation performance analysis
for all the four datasets  as well as for all the feature sets, performance metrics averaged over 10-fold as well as  <dig> sample sets are given in table  <dig>  moreover, to analyze the trend in accuracies, performance metrics are also plotted in line graphs .
it is observed that the sensitivities are less as compared to the specificities .
further, higher accuracies are observed for the dataset having resistant proteins with <90% pair-wise sequence identity, whereas lower accuracies are observed for the dataset having resistant proteins with <40% pair-wise sequence identity .
though the specificities are observed almost unchanged, sensitivities are observed to be increased with increase in the percentage of pair-wise sequence identity in the positive dataset .
besides, it is seen that the most of the performance metrics for dpc and ctd feature sets are higher as compared to the other feature sets .
in particular, overall accuracy , mcc  and auc-roc  are observed to be highest for dpc feature set.
since the number of sequences in the positive dataset having sequences with <90% pair-wise sequence identity is larger as compared to the dataset having sequences with <40% pair-wise sequence identity, the former one is used in subsequent analyses.table  <dig> estimates of different performance metrics for svm with rbf kernel in discriminating resistant from non-resistant proteins, under all the feature sets as well as different percentage of sequence identity in the positive dataset
id: maximum percentage of pair-wise sequence identity present in the positive dataset
sn sensitivity, sp specificity, ac accuracy, pre precision, mcc matthew’s correlation coefficient, auc-roc area under roc curves
fig.
<dig> performance metrics of svm with rbf kernel for different feature sets and different percentage of pair-wise sequence identity in the positive set.
it can be seen that the performance metrics are higher for dpc feature set as compared to other feature sets, irrespective of the percentage of sequence identity in the positive dataset
analysis of mechanism-based classification
the values of performance metrics, with regard to classification of resistant proteins involved in target-based mechanism and detoxification-based mechanism, measured over loocv as well as  <dig> sample sets are given in table  <dig>  performance metrics for all the feature sets are observed ≥90% and are found to be highest in case of dpc feature set.
more specifically, overall accuracy for the dpc feature set is observed >97%, with >95% mcc and >97% auc-roc.
though the number of features for aac and paac feature sets are almost same, classification accuracies for aac feature set are seen to be higher than that of paac feature set.
since the sensitivity and specificity are >90%, it is inferred that hardly one sequence is misclassified in each category .table  <dig> estimates of performance metrics for classification of detoxification and target-based resistant proteins, under different feature sets
sn sensitivity, sp specificity, ac accuracy, pre precision, mcc matthew’s correlation coefficient, auc-roc area under roc curves
discriminating target-based resistant proteins from non-resistant proteins
with regard to classification of target-based resistant proteins and non-resistant proteins, performance metrics over loocv and  <dig> sample sets  are given in table  <dig>  the values of performance metrics are observed to be higher for dpc feature set.
specifically, accuracies in terms of all the performance metrics are observed ≥90% for dpc feature set, whereas the values of mcc and auc-roc for rest of the feature sets are observed to be <90%.table  <dig> estimates of performance metrics for discriminating target-based resistant proteins from non-resistant proteins, under different features
sn sensitivity, sp specificity, ac accuracy, pre precision, mcc matthew’s correlation coefficient, auc-roc area under roc curves
discriminating detoxification-based resistant proteins from non-resistant proteins
the classification was also made between  <dig> detoxification-based resistant proteins and  <dig> non-resistant proteins, by using svm with rbf kernel.
performances metrics were computed over  <dig> folds of cross validation as well as  <dig> sample sets  are presented in table  <dig>  it is observed that the accuracies are higher for dpc feature set and lower for aac feature set.
in particular, the values of all the performances metrics for both ctd and dpc feature sets are ≥90% .
barring sensitivity, the values of performance metrics in discriminating the detoxification-based resistant proteins from non-resistant proteins  are higher as compared to that of discriminating target-based resistant proteins from non-resistant proteins .table  <dig> estimates of different performance metrics for discriminating detoxification-based resistant proteins from non-resistant proteins
sn sensitivity, sp specificity, ac accuracy, pre precision, mcc matthew’s correlation coefficient, auc-roc area under roc curves
comparative analysis
for comparing the proposed approach with blast algorithms, we prepared two different datasets.
the first dataset contains  <dig> resistant proteins  and randomly drawn  <dig> non-resistant proteins , and the second dataset contains  <dig> resistant proteins  and randomly drawn  <dig> non-resistant proteins .
furthermore, performance of the proposed approach was compared based on dpc feature set only as higher accuracies were obtained for this feature set as compared to the other feature sets.
in both the datasets, no hits were found for most of the query sequences with e-values  <dig>  and  <dig>  however, hits were found for all the query sequences with e-value  <dig>  therefore, comparison was made based on e-value  <dig> only, and the accuracies averaged over 10-folds are given in table  <dig>  it is observed that the overall accuracies of the proposed approach are ~10% higher than that of blastp, psi-blast and delta-blast, in both datasets .
though, true positive rates  of the blast algorithms are higher, false positive rates  are much lower at the same time.
among the blast algorithms, delta-blast performed better than both blastp and psi-blast, with both the datasets .
barring sensitivity, the proposed approach performed better than blast algorithms in terms of all the performance metrics.
it is further seen that the specificities are higher for the first dataset as compared to the second dataset.table  <dig> performance metrics for the proposed approach, blast, psi-blast and delta-blast, in discriminating the resistant proteins from non-resistant proteins, where the positive dataset consists of <40%  and <90%  pair-wise sequence identity
here, auc-roc values were not computed, as in blast algorithms accuracies are computed based on number of hits
sn sensitivity, sp specificity, ac accuracy, pre precision, mcc matthew’s correlation coefficient
performance analysis based on independent test dataset
both the datasets mentioned in “comparative analysis” section were used to train the model for prediction of the level  of each test sequence.
furthermore, none of the test sequences were present in the training set.
it is observed that  <dig> out of  <dig> are correctly predicted while first dataset is used as training set .
on the other hand, all the  <dig> instances are correctly identified as insecticide resistant proteins with second dataset as training set .
besides, it is seen that most of the sequences are correctly predicted with > <dig>  probabilities irrespective of the training datasets .
more clearly,  <dig> test sequences of cytochrome p <dig> and  <dig> sequences of ache are misclassified in the first training dataset .table  <dig> performance of the proposed approach based on an independent dataset of  <dig> insecticide resistant proteins
fig.
<dig> heat map of the probabilities with which  <dig> test sequences are predicted in two different training datasets.
all the  <dig> sequences are correctly predicted as resistant proteins in the second training dataset, whereas  <dig> are correctly predicted with the first training dataset.
it is further seen that most of the test sequences are correctly predicted with high probabilities
online prediction server: dirprot
a web server dirprot has been developed to discriminate the insecticide resistant proteins from non-resistant proteins.
this server has been trained with the svm  for prediction of insecticide resistant proteins based on dpc feature set.
the web pages showing the execution and results for an example dataset are shown in fig.
5a and b respectively.
help pages are also provided to guide the user regarding generation of features, prediction method and input–output.
the sequences in fasta format along with the annotations and probabilities with which they are predicted as resistance proteins are shown in the result page.
for reproducible research, the trained datasets are also provided in the server.
the prediction server is made freely accessible at http://cabgrid.res.in:8080/dirprot for academic users.fig.
<dig>
a server page of dirprot, b result page after execution with an example dataset.
the result page is displayed in a tabular form, where the last column is the probabilities with which the each sequences are predicted as insecticide-resistant proteins
discussion
extensive use of chemical insecticides has been selecting resistant population of insect species to different insecticides, worldwide  <cit> .
around  <dig> insect species have been reported to resist different insecticides till the end of  <dig>  <cit> .
insecticidal resistance has been associated with the genetic changes in insects.
for instance, a mutation in an insect can alter the behavior, metabolism and physiology by which insect may gain advantage in resisting to different insecticides  <cit> .
most of the earlier studies are dealt with the mutational changes associated with the insecticide resistance.
though insecticide resistance is an important researchable issue, there is no computational tool available for prediction of insecticide resistant proteins.
therefore, we made an attempt to present the first computational approach for prediction of insecticide resistant proteins.
we considered four different categories of insecticide resistant proteins corresponding to four different classes of insecticide resistance genes viz., cytochrome p <dig>  ache, rdl and kdr.
the leucine content was predominantly found in all the four categories of proteins, which has been reported to play an important role in insecticide resistance.
for instance, prince et al.
<cit>  reported that leucine-rich repeat receptor-like kinase “brassinosteroid insensitive1-associated kinase1” contributes to the innate immunity to aphids in arabidopsis.
the valine to leucine  and the leucine to isoleucine mutations  were identified in three pesticide-resistant strains of bed bug   <cit> .
further, the composition of tryptophan which has been reported to present in the active site that interacts with trimethyl-ammonium cationic group of ache was found lowest  <cit> .
hassani et al.
<cit>  described that lysine and tryptophan  are the most reactive residues that play important role in disrupting the function of neuronal sodium channels by ts gamma, which is the most potent neurotoxin in the venom of the brazilian scorpion tityus serrulatus.
for classification of insecticide resistant and non-resistant proteins, initially the sequences were transformed into numeric feature vectors based on different feature generation techniques viz., aac, dpc, paac, acf and ctd.
the encoded numeric vectors were then used as input in binary svm classifier.
prediction accuracies were found to be higher for rbf kernel as compared to the other three kernels of svm.
further, the classification accuracies were found higher for dpc feature set as compared to the other feature sets, which may be due to the fact that in dpc the local ordering of amino acids were taken into account  <cit> .
furthermore, in cross validation analysis , the sensitivity was found to be increased with increase in the percentage of pair-wise sequence identity in the positive dataset.
this may be due to the fact that with increase in the pair-wise sequence identity in the positive dataset, it is less-likely that a positive sequence will be misclassified in the negative dataset.
the accuracy in discriminating the target-based and detoxification-based resistance proteins from non-resistant proteins was also found to be higher.
besides, higher discrimination accuracy was also observed between target-based and detoxification-based resistance proteins.
thus, it can be inferred that the composition of di-peptides are not only different between resistant and non-resistant proteins but also among insecticide resistant proteins involved in different insecticide resistance mechanisms.
the performance of the proposed approach was compared with blast, psi-blast and delta-blast algorithms.
though, prediction was made for three e-values i.e.,  <dig> ,  <dig> and  <dig>  no hits were found for most of the query sequences  for the first two e-values.
thus performance metrics were computed based on e-value  <dig> only, which is also the default e-value in blast algorithms.
though the resistant proteins were predicted with higher accuracy, the specificities were found much lower.
it was also found that the specificities are higher for the first dataset as compared to the second dataset.
one of the possible reasons for this may be that when the pair-wise sequence identity is <40% in the positive class , sequence similarity between the classes will be less.
on the other hand, when the pair-wise sequence identity is <90% in the positive class, sequence similarity between the positive and negative classes will be more by which the likelihood of a sequence of the negative class to be predicted in the positive class will be more and vice versa.
in terms of overall accuracy, the proposed approach outperformed all the three variations of blast algorithm.
among the blast algorithms, delta-blast performed better followed by psi-blast and blast.
the performance of the proposed approach was also assessed using an independent test dataset consisting of  <dig> resistant protein sequences .
out these  <dig> sequences, all were correctly predicted when the pair-wise sequence identity was <90% in the positive dataset of training set, whereas  <dig> were correctly predicted in for the training dataset having positive sequences with <40% pair-wise sequence identity.
nevertheless, the proposed approach achieved higher accuracy for predicting the insecticide resistant proteins.
CONCLUSIONS
this paper presents the first computational approach for predicting the insecticide resistant proteins.
based on this approach, a web server has also been developed that can be easily used by the scientists and researchers to computationally identify the insecticide resistant proteins.
the proposed computational approach is believed to supplement the wet-lab experiments for identifying and targeting the insecticide resistant proteins to develop dynamic and efficient insecticides.
additional file
additional file 1: it contains the list of insect species and corresponding insecticide resistant gene types that were used in this study.
this file also contains the  <dig> insecticide resistant protein sequences that were used as independent dataset.
abbreviations
aacamino acid composition
acaccuracy
acfauto correlation function
acheacetylcholinesterases
auc-rocarea under roc curves
ctdcomposition-trainsition-distribution
dpcdipeptide composition
fnfalse negative
fpfalse positive
gabaγ-amino butyric acid
gstsglutathione s-transferases
kdrknock down resistance
mccmatthew’s correlation coefficient
oshoptimal separating hyper-plane
paacpseudo amino acid composition
preprecision
rbfradial basis function
rdlresistance to dieldrin
rocreceiving operating characteristics
snsensitivity
spspecificity
svmsupport vector machine
tntrue negative
tptrue positive

 insecticides are used to control the insects affecting the agricultural crops, parasitizing livestock, as well as to eradicate the pests transmitting dangerous infectious diseases. however, frequent application of insecticides has resulted in the resurgence of pests and appearance of resistant pest species. insecticide resistance is the heritable change in the sensitivity of a pest population that is reflected in the repeated failure of a product to achieve the expected level of control when used according to the level of recommendation for that pest species. several studies have indicated the involvement of multiple genes in conferring the resistance to many insect species. thus, characterization of these genes is useful to understand the development of resistance and designing new strategies to minimize the development of insecticide resistance. three major mechanisms are involved in insecticide resistance: detoxification of insecticides through alteration in the activities of enzymes like esterase, oxidases or glutathione s-transferases that prevents the insecticide from reaching to its site of action, insensitivity of the insecticide target proteins, because of which the insecticide no longer binds to its target and reduction in insecticide uptake due to decrease in permeability of insect cuticle. though there is evidence of alteration in cuticular penetration, most of the studies have focused and evaluated the target site insensitivity and detoxification of insecticides mechanisms. moreover, these two mechanisms have been reported to cover a wide range of resistance levels to almost all available insecticides. the cytochrome p family of genes in insect play an important role in the detoxification of insecticides resulted in the development of resistance to insecticides. besides, gsts have also been reported to be involved in the detoxification of insecticides. as far as target-based mechanism is concerned, there are three main targets for conventional insecticides viz., gaba -gated chloride ion channel, voltage-gated ion channel and acetylcholinesterases. the gaba receptor is the site of target for cyclodiene insecticides, where the resistance to dieldrin is conferred by the change of a single amino acid in gaba-gated chloride ion channel encoded by rdl gene. further, knockdown resistance is one of the major forms of resistance to ddt and pyrethroid insecticides, which is associated with mutations in the voltage-gated sodium channel. acetylcholinesterase in nerve synapses is the target protein for the insecticides like organophosphorus and carbamate. the point mutation in the insecticide-binding site of ache has been identified as the cause of insensitivity to these insecticides. the above mentioned works help enable to understand the molecular mechanisms involved in the insecticide resistance. further, the analysis of bio-molecules involved in this phenomenon has confirmed the importance of single genes in target site resistance and involvement of multi-gene families like cytochrome p in metabolic resistance. several studies on the effects of mutational changes in target proteins on insecticide resistance aid to the knowledge on the insect proteins involved in this process. for instance, riveron et al. demonstrated that the single amino acid change in an up regulated gst gene, gste confers high level of metabolic resistance to ddt in the malaria vector anopheles funestus. in another study, nwane et al. identified that two mutations at position of the s transmembrane segment of domain ii in the voltage gated sodium channel i.e., leucine to a phenylalanine or to a serine confers resistance to ddt and pyrethroid insecticides in anopheles gambiae. in the recent past, several studies have identified species-specific insecticide resistant genes through transcriptome and expression profile analysis. hsu et al. identified p gst, coe-related genes in bactrocera doralis, representing three major enzyme families involved in insecticide metabolism and resistance. in another study, p gst and ces-specific genes of liposcelis bostrychophila were reported to be involved in insecticide resistance, through transcriptome and differential gene expression analysis. recently, cui et al. identified relevant genes in response to flubendiamide insecticide in asian corn borer, through de novo transcriptome and expression-profile analysis. though the transcriptome and expression profile analysis is one way of identifying the resistance genes, it is species specific. moreover the expression profile analysis is expensive as well as time consuming. thus, development of a computational tool for identifying the resistant genes independent of the species and economically as well would help in augmenting the research related to the identification of insecticide resistant genes. however, no computational tool is reported till date for the discrimination of insecticide resistant proteins from the proteins that do not confer resistance. keeping this in view, we propose a computational approach to discriminate the insecticide resistant proteins from non-resistant proteins. the developed computational approach can be used for identification of the resistant proteins across species as well as with minimum resource. we have also developed an online prediction server that can be easily used by experimental scientist and researchers to predict an unknown protein sequence as either insecticide-resistant or non-resistant protein. moreover, computational identification of insecticide resistant proteins will supplement the efforts needed to develop insecticides in targeting the resistance proteins. methods collection and processing of data in this study, protein sequences corresponding to four important groups of insecticide resistant genes viz., cytochrome p kdr, rdl and ache were collected from insecticide resistance gene database. we considered these four categories of genes because they represent important families of insecticide resistant genes which are resistant to commonly used insecticides. besides, the resistant protein sequences were reported to be involved in two important resistance mechanisms viz., detoxification-based and target-based. further, target-based resistant proteins are confined to three main targets of insecticides i.e., ache, gaba-gated chloride ion channel and voltage-gated sodium channel. a total of sequences belonging to insect species were collected. initially, we removed the sequences having non-standard residues. then, four positive sets having and sequences were prepared, where the maximum pair-wise sequence identities were 40%, 60%, 70% and 90% respectively. the sequences with more than considered level of pair-wise sequence identity were removed using cdhit. for negative set, protein sequences of the considered species were collected from the uniprot database. for the species acyrthosiphon pisum and tribolium castaneum, only the reviewed sequences were collected, as large number of sequence are present in uniprot for these two species. on the other hand, all the sequences available for remaining nine species were collected. after removing the sequences having non-standard residues as well as the identical sequences, a total of sequences were obtained. further, to avoid homologous bias in the negative dataset, sequences with >40% pair-wise identity were removed using cdhit. finally a dataset with sequences was obtained and considered as the negative dataset. feature generation protein sequences are the strings of amino acid residues, and hence they need to be mapped onto numeric feature vectors before being used as input in machine learning classifier. in this study, amino acid composition, di-peptide composition, pseudo amino acid composition, composition-transition-distribution and auto correlation function were used to transform the protein sequences into numeric feature vectors. amino acid composition aac is a basic feature of protein sequence, which is closely associated with its attributes, such as sub-cellular location, secondary structure content and domain. aac consists of discrete numbers, each of which represents the frequency of the native amino acids in a protein sequence. based on the aac, each protein sequence was encoded into a 20-dimensional numerical vector. di-peptide composition one of the limitations of aac is that it does not take into account the local order information of amino acids in the protein. on the other hand, dpc, which gives a fixed pattern length of, encapsulates the global information about each protein sequence and the order it contains. for any di-peptide, its composition was computed as the ratio of the frequency of that di-peptide to the total number possible di-peptide in the protein sequence. pseudo amino acid composition the concept of paac was originally introduced by chou for predicting the protein sub-cellular locations and membrane protein types. based on the conventional aac, chou proposed a set of discrete numbers to take into account the sequence order effects. paac has been proven to be an extremely effective feature in many proteins and protein-related systems. the paac for each protein sequence can be represented by a -dimensional vector for d-tier correlation factor. here, the paac was extracted for 1st-tier correlation only, by which each sequence was transformed into a 21-dimensional numeric vector. for further details, one can refer to. composition-transition-distribution the ctd feature was introduced by dubchak et al. for predicting protein folding classes. thereafter, the ctd feature has been adopted by many researchers for protein function and structure studies. in ctd feature, composition is the number of amino acids of a particular type divided by the total number of amino acids. transition characterizes the frequency percentage with which amino acids of a particular type is followed by other amino acids. distribution measures the chain length within which the first 25%, 50%, 75% and 100% of the amino acids of a particular type is located respectively. based on the ctd feature, each protein sequence of length l was encoded to a l+{l*/2} + -dimension numeric vector. auto correlation function sequence autocorrelation-based features assume that the disturbances in each area are systematically related to those in adjacent areas. this concept helps to analyze the dependency among the features of sequences in each location. autocorrelation features were computed based on the distribution of amino acid properties along the sequence, using all the amino acid indices available in aaindex database. in this feature encoding, for an autocorrelation of order n, each sequence was transformed into a numeric vector of length 531*n. supervised learning technique for classification purpose we used the support vector machine, which is a nonparametric algorithm developed by vapnik. it is a very promising and popular method for pattern recognition that has been widely used for prediction purpose in the field of bioinformatics. it is proven to be very efficient in many biological analyses due to their ability to handle noise and large input dataset. a brief description about the working principle of svm is described as follows: consider a binary classification problem with n samples or input vectors x i ∈ r d,, where x i with class levels y i ∈ {−  1} can be considered as the i th protein or vector defined in d-dimensional space. in present work, refers to resistant class and − represents non-resistant class. the objective here is to construct a binary classifier from the available sample that has less probability of misclassifying future sample. non-linear svm maps input vectors x i' s into high dimensional feature space and constructs an optimal separating hyper-plane that maximizes the distance between hyper-plane and nearest data points of each class in the space. mathematically, the hyper-plane is represented as y = sgn, where w represents a weight vector that can map training data in the input space to the outer space and b represents bias. for a two class problem, it can be formulated as wtxi+b≥1ifyi=1wtxi+b≤−1ifyi=− the svm training procedure involves optimization of convex quadratic problem i.e., with lagrangian multipliers α i ≥  maximize∑i=1nαi−12∑i=1n∑j=1nαiαjyiyjkxixj subject to the constraints 0 ≤ α i ≤ c  and ∑i=1nαiyi= where c is the regularization parameter that controls trade-off between margin and classification error. the x j' s are called support vectors only if corresponding α j >  after the svm has been trained, the decision function for classification of query sequence can be formulated as fx=sgn∑i=1nyiαikx.xi+b. the choice of the proper kernel function k is important to train svm model because the power of svm comes from the kernel representation that allows the nonlinear mapping of input space to a higher dimensional feature space. in this work, four commonly used kernel functions viz., linear, polynomial d), radial basis and sigmoid ) were used, where r, d, γ > are the kernel parameters. validation of the model cross-validation procedure has been widely accepted for assessing the performance of classifiers. thus, we used the 10-fold cross-validation to assess the performance of our approach. it was carried out by partitioning the dataset into approximately equal-sized sets at random, where nine partitions were used to train the model and the remaining one part was used to assess the model accuracy. this process was repeated times in such a way that each partition was tested once in the model. performance evaluation different performance metrics viz., sensitivity, specificity, accuracy, precision and matthew’s correlation coefficient were used to measure the accuracy of the developed prediction approach. the sn, sp, ac, pre and mcc parameters are defined as:sn = tp/, sp = tn/, ac = /, pre = tp/, mcc=tp×tn−fp×fn/tp+fn×tp+fp×tn+fn×tn+fp. true positive is the number of resistant proteins correctly predicted as resistant proteins, true negative is the number of non-resistant proteins correctly predicted as non-resistant proteins, false negative is the number of resistant proteins incorrectly predicted as non-resistant proteins and false positive is the number of non- resistant proteins incorrectly predicted as resistant proteins. besides the above mentioned performance metrics, area under receiving operating characteristic curve was also used to measure the predictive ability. for given false positive rate and true positive rate at different threshold values, the auc-roc was computed as ∑i1−βi.Δα+1/2Δ1−β.Δα, where Δ =  −, Δα = α i − α i −  and i = , …, m. a subroutine in r programming language was written to compute the values of these performance metrics. training and testing datasets using four positive sets and one negative set, four datasets were prepared that consists of both positive and negative sequences. here each dataset contains a different positive set and the same negative set. all the four datasets are highly unbalanced as the number sequences present in one class is much larger than the other class. to avoid biasness towards the non-resistant class while predicting using machine learning classifier like svm, balanced datasets were prepared that consists of same number of sequences from both the classes, where the sequences of the major class were randomly drawn from the available sequences of the major class. for instance, first balanced dataset contains positive and negative sequences, where the negative sequences were randomly drawn from negative sequences. as the generalized predictive ability cannot be assured based on a single dataset, sample sets were prepared, where each sample set consists of same number of positive and negative instances. further, in each sample set, a 10-fold cross validation procedure was adopted. the performance metrics were computed by taking average over the folds as well as over sample sets. mechanism-based classification the insecticide resistance mechanism can be broadly categorized into two types, viz., target-based mechanism and detoxification-based mechanism. the rdl, kdr and ache genes come under target-based and cytochrome p genes come under detoxification-based mechanism. to test whether the genes under these two categories are different or not, a binary classification was carried out by employing svm, where sequences from target-based and sequences from detoxification-based category were used. similar to the classification of resistant and non-resistant proteins, sample sets were prepared where each sample set consists of sequences from each class. since, there are sequences in the detoxification-based category, sequences were randomly drawn each time. as the number of sequences in each sample is not large, leave-one-out cross validation technique was adopted for classification of detoxification- and target-based resistant proteins. here, detoxification-based category was considered as positive class and target-based category as negative class. comparison with blast algorithm performance of the proposed approach was also compared with that of blastp, psi-blast and delta-blast, which are powerful algorithms to detect protein homologs. further, comparison was made through 10-fold cross validation technique. for cross validation, offline blast software was used with blastp, psiblast and deltablast modules/programs in which the training set for each fold of cross validation was defined as the database and sequences of the corresponding test set were used as query. each query sequence was predicted as the resistant or non-resistant category based on the top hit found in the blast search. three different e-values i.e.,, and were used to assess the performance of the blastp, psi-blast and delta-blast. furthermore, performance of the proposed approach was compared based on best feature set with which higher accuracies were obtained as compared to the other feature sets. performance evaluation using independent dataset to assess the generalized predictive ability of the proposed approach, its performance was further tested using an independent test dataset. the independent dataset was collected based on published literature that includes cytochrome p kdr, rdl and ache proteins. specifically, cytochrome p genes were reported by hsu et al.. out of we used as they are available in ncbi. similarly, kdr, rdl and ache genes were collected from ncbi, based on the study of zuo et al., wondji et al. and li and han respectively. sequences of the independent test set are provided in additional file development of prediction server an online prediction server was developed using html and php, where the combination of best feature set and classifier was used. a developed r-code was executed in background upon submission of the sequences in fasta format to the server. the user has to submit the protein sequences having only standard amino acid residues. this server can be used to predict the likelihood of any unknown protein sequence being an insecticide resistant protein with certain probability.  analysis of amino acid compositions the composition of amino acids in four different groups of insecticide resistant proteins is shown in fig.  it is seen that the proportions of leucine are higher, whereas the proportions of cystene and tryptophan are lower in all the four categories.fig. composition of amino acids in all the four categories of insecticide resistant proteins. it is observed that proportions of leucine are higher, whereas proportions of cystene and tryptophan are lower in all the four categories analysis of kernel functions based on a sample dataset consisting of positive and negative sequences that were drawn randomly from the available positive and negative sequences, performance of svm was analyzed. roc curves for all the four kernels as well as for all the five feature sets are shown in fig. 2a and the corresponding auc-roc values are shown in bar plots. from the roc curves it is not clear that which kernel is better, whereas from auc-roc plots it is clear that the values of auc-roc are higher for the rbf kernel, irrespective of the feature set used. though in rbf kernel the auc-roc for acf feature set is highest, it is difficult to choose the best feature set while other three kernels are taken into account. therefore, all the feature sets and the rbf kernel were used for further analysis.fig. a roc curves of svm for different kernels and features, b bar plots of corresponding auc-roc values. it is seen that the auc-roc values are higher for rbf kernel as compared to other kernels cross-validation performance analysis for all the four datasets as well as for all the feature sets, performance metrics averaged over 10-fold as well as sample sets are given in table  moreover, to analyze the trend in accuracies, performance metrics are also plotted in line graphs. it is observed that the sensitivities are less as compared to the specificities. further, higher accuracies are observed for the dataset having resistant proteins with <90% pair-wise sequence identity, whereas lower accuracies are observed for the dataset having resistant proteins with <40% pair-wise sequence identity. though the specificities are observed almost unchanged, sensitivities are observed to be increased with increase in the percentage of pair-wise sequence identity in the positive dataset. besides, it is seen that the most of the performance metrics for dpc and ctd feature sets are higher as compared to the other feature sets. in particular, overall accuracy, mcc and auc-roc are observed to be highest for dpc feature set. since the number of sequences in the positive dataset having sequences with <90% pair-wise sequence identity is larger as compared to the dataset having sequences with <40% pair-wise sequence identity, the former one is used in subsequent analyses.table estimates of different performance metrics for svm with rbf kernel in discriminating resistant from non-resistant proteins, under all the feature sets as well as different percentage of sequence identity in the positive dataset id: maximum percentage of pair-wise sequence identity present in the positive dataset sn sensitivity, sp specificity, ac accuracy, pre precision, mcc matthew’s correlation coefficient, auc-roc area under roc curves fig. performance metrics of svm with rbf kernel for different feature sets and different percentage of pair-wise sequence identity in the positive set. it can be seen that the performance metrics are higher for dpc feature set as compared to other feature sets, irrespective of the percentage of sequence identity in the positive dataset analysis of mechanism-based classification the values of performance metrics, with regard to classification of resistant proteins involved in target-based mechanism and detoxification-based mechanism, measured over loocv as well as sample sets are given in table  performance metrics for all the feature sets are observed ≥90% and are found to be highest in case of dpc feature set. more specifically, overall accuracy for the dpc feature set is observed >97%, with >95% mcc and >97% auc-roc. though the number of features for aac and paac feature sets are almost same, classification accuracies for aac feature set are seen to be higher than that of paac feature set. since the sensitivity and specificity are >90%, it is inferred that hardly one sequence is misclassified in each category.table estimates of performance metrics for classification of detoxification and target-based resistant proteins, under different feature sets sn sensitivity, sp specificity, ac accuracy, pre precision, mcc matthew’s correlation coefficient, auc-roc area under roc curves discriminating target-based resistant proteins from non-resistant proteins with regard to classification of target-based resistant proteins and non-resistant proteins, performance metrics over loocv and sample sets are given in table  the values of performance metrics are observed to be higher for dpc feature set. specifically, accuracies in terms of all the performance metrics are observed ≥90% for dpc feature set, whereas the values of mcc and auc-roc for rest of the feature sets are observed to be <90%.table estimates of performance metrics for discriminating target-based resistant proteins from non-resistant proteins, under different features sn sensitivity, sp specificity, ac accuracy, pre precision, mcc matthew’s correlation coefficient, auc-roc area under roc curves discriminating detoxification-based resistant proteins from non-resistant proteins the classification was also made between detoxification-based resistant proteins and non-resistant proteins, by using svm with rbf kernel. performances metrics were computed over folds of cross validation as well as sample sets are presented in table  it is observed that the accuracies are higher for dpc feature set and lower for aac feature set. in particular, the values of all the performances metrics for both ctd and dpc feature sets are ≥90%. barring sensitivity, the values of performance metrics in discriminating the detoxification-based resistant proteins from non-resistant proteins are higher as compared to that of discriminating target-based resistant proteins from non-resistant proteins.table estimates of different performance metrics for discriminating detoxification-based resistant proteins from non-resistant proteins sn sensitivity, sp specificity, ac accuracy, pre precision, mcc matthew’s correlation coefficient, auc-roc area under roc curves comparative analysis for comparing the proposed approach with blast algorithms, we prepared two different datasets. the first dataset contains resistant proteins and randomly drawn non-resistant proteins, and the second dataset contains resistant proteins and randomly drawn non-resistant proteins. furthermore, performance of the proposed approach was compared based on dpc feature set only as higher accuracies were obtained for this feature set as compared to the other feature sets. in both the datasets, no hits were found for most of the query sequences with e-values and however, hits were found for all the query sequences with e-value therefore, comparison was made based on e-value only, and the accuracies averaged over 10-folds are given in table  it is observed that the overall accuracies of the proposed approach are ~10% higher than that of blastp, psi-blast and delta-blast, in both datasets. though, true positive rates of the blast algorithms are higher, false positive rates are much lower at the same time. among the blast algorithms, delta-blast performed better than both blastp and psi-blast, with both the datasets. barring sensitivity, the proposed approach performed better than blast algorithms in terms of all the performance metrics. it is further seen that the specificities are higher for the first dataset as compared to the second dataset.table performance metrics for the proposed approach, blast, psi-blast and delta-blast, in discriminating the resistant proteins from non-resistant proteins, where the positive dataset consists of <40% and <90% pair-wise sequence identity here, auc-roc values were not computed, as in blast algorithms accuracies are computed based on number of hits sn sensitivity, sp specificity, ac accuracy, pre precision, mcc matthew’s correlation coefficient performance analysis based on independent test dataset both the datasets mentioned in “comparative analysis” section were used to train the model for prediction of the level of each test sequence. furthermore, none of the test sequences were present in the training set. it is observed that out of are correctly predicted while first dataset is used as training set. on the other hand, all the instances are correctly identified as insecticide resistant proteins with second dataset as training set. besides, it is seen that most of the sequences are correctly predicted with > probabilities irrespective of the training datasets. more clearly, test sequences of cytochrome p and sequences of ache are misclassified in the first training dataset.table performance of the proposed approach based on an independent dataset of insecticide resistant proteins fig. heat map of the probabilities with which test sequences are predicted in two different training datasets. all the sequences are correctly predicted as resistant proteins in the second training dataset, whereas are correctly predicted with the first training dataset. it is further seen that most of the test sequences are correctly predicted with high probabilities online prediction server: dirprot a web server dirprot has been developed to discriminate the insecticide resistant proteins from non-resistant proteins. this server has been trained with the svm for prediction of insecticide resistant proteins based on dpc feature set. the web pages showing the execution and results for an example dataset are shown in fig. 5a and b respectively. help pages are also provided to guide the user regarding generation of features, prediction method and input–output. the sequences in fasta format along with the annotations and probabilities with which they are predicted as resistance proteins are shown in the result page. for reproducible research, the trained datasets are also provided in the server. the prediction server is made freely accessible at http://cabgrid.res.in:8080/dirprot for academic users.fig. a server page of dirprot, b result page after execution with an example dataset. the result page is displayed in a tabular form, where the last column is the probabilities with which the each sequences are predicted as insecticide-resistant proteins discussion extensive use of chemical insecticides has been selecting resistant population of insect species to different insecticides, worldwide. around insect species have been reported to resist different insecticides till the end of. insecticidal resistance has been associated with the genetic changes in insects. for instance, a mutation in an insect can alter the behavior, metabolism and physiology by which insect may gain advantage in resisting to different insecticides. most of the earlier studies are dealt with the mutational changes associated with the insecticide resistance. though insecticide resistance is an important researchable issue, there is no computational tool available for prediction of insecticide resistant proteins. therefore, we made an attempt to present the first computational approach for prediction of insecticide resistant proteins. we considered four different categories of insecticide resistant proteins corresponding to four different classes of insecticide resistance genes viz., cytochrome p ache, rdl and kdr. the leucine content was predominantly found in all the four categories of proteins, which has been reported to play an important role in insecticide resistance. for instance, prince et al. reported that leucine-rich repeat receptor-like kinase “brassinosteroid insensitive1-associated kinase1” contributes to the innate immunity to aphids in arabidopsis. the valine to leucine and the leucine to isoleucine mutations were identified in three pesticide-resistant strains of bed bug. further, the composition of tryptophan which has been reported to present in the active site that interacts with trimethyl-ammonium cationic group of ache was found lowest. hassani et al. described that lysine and tryptophan are the most reactive residues that play important role in disrupting the function of neuronal sodium channels by ts gamma, which is the most potent neurotoxin in the venom of the brazilian scorpion tityus serrulatus. for classification of insecticide resistant and non-resistant proteins, initially the sequences were transformed into numeric feature vectors based on different feature generation techniques viz., aac, dpc, paac, acf and ctd. the encoded numeric vectors were then used as input in binary svm classifier. prediction accuracies were found to be higher for rbf kernel as compared to the other three kernels of svm. further, the classification accuracies were found higher for dpc feature set as compared to the other feature sets, which may be due to the fact that in dpc the local ordering of amino acids were taken into account. furthermore, in cross validation analysis, the sensitivity was found to be increased with increase in the percentage of pair-wise sequence identity in the positive dataset. this may be due to the fact that with increase in the pair-wise sequence identity in the positive dataset, it is less-likely that a positive sequence will be misclassified in the negative dataset. the accuracy in discriminating the target-based and detoxification-based resistance proteins from non-resistant proteins was also found to be higher. besides, higher discrimination accuracy was also observed between target-based and detoxification-based resistance proteins. thus, it can be inferred that the composition of di-peptides are not only different between resistant and non-resistant proteins but also among insecticide resistant proteins involved in different insecticide resistance mechanisms. the performance of the proposed approach was compared with blast, psi-blast and delta-blast algorithms. though, prediction was made for three e-values i.e.,, and no hits were found for most of the query sequences for the first two e-values. thus performance metrics were computed based on e-value only, which is also the default e-value in blast algorithms. though the resistant proteins were predicted with higher accuracy, the specificities were found much lower. it was also found that the specificities are higher for the first dataset as compared to the second dataset. one of the possible reasons for this may be that when the pair-wise sequence identity is <40% in the positive class, sequence similarity between the classes will be less. on the other hand, when the pair-wise sequence identity is <90% in the positive class, sequence similarity between the positive and negative classes will be more by which the likelihood of a sequence of the negative class to be predicted in the positive class will be more and vice versa. in terms of overall accuracy, the proposed approach outperformed all the three variations of blast algorithm. among the blast algorithms, delta-blast performed better followed by psi-blast and blast. the performance of the proposed approach was also assessed using an independent test dataset consisting of resistant protein sequences. out these sequences, all were correctly predicted when the pair-wise sequence identity was <90% in the positive dataset of training set, whereas were correctly predicted in for the training dataset having positive sequences with <40% pair-wise sequence identity. nevertheless, the proposed approach achieved higher accuracy for predicting the insecticide resistant proteins.  this paper presents the first computational approach for predicting the insecticide resistant proteins. based on this approach, a web server has also been developed that can be easily used by the scientists and researchers to computationally identify the insecticide resistant proteins. the proposed computational approach is believed to supplement the wet-lab experiments for identifying and targeting the insecticide resistant proteins to develop dynamic and efficient insecticides. additional file additional file 1: it contains the list of insect species and corresponding insecticide resistant gene types that were used in this study. this file also contains the insecticide resistant protein sequences that were used as independent dataset. abbreviations aacamino acid composition acaccuracy acfauto correlation function acheacetylcholinesterases auc-rocarea under roc curves ctdcomposition-trainsition-distribution dpcdipeptide composition fnfalse negative fpfalse positive gabaγ-amino butyric acid gstsglutathione s-transferases kdrknock down resistance mccmatthew’s correlation coefficient oshoptimal separating hyper-plane paacpseudo amino acid composition preprecision rbfradial basis function rdlresistance to dieldrin rocreceiving operating characteristics snsensitivity spspecificity svmsupport vector machine tntrue negative tptrue positive 
the proposed approach not only outperformed blastp, psi-blast and delta-blast algorithms, but also achieved >92% accuracy while assessed using an independent dataset of insecticide resistant proteins. the encoded numeric vectors were then used as input in support vector machine for classification of insecticide resistant and non-resistant proteins. to distinguish the insecticide resistant proteins from non-resistant proteins, no computational tool is available till date. the proposed approach achieved an overall accuracy of >90% in discriminating resistant from non-resistant proteins. further, the two classes of resistant proteins i.e., detoxification-based and target-based were discriminated from non-resistant proteins with >95% accuracy. this paper presents the first computational approach for discriminating the insecticide resistant proteins from non-resistant proteins. five different sets of feature viz., amino acid composition, di-peptide composition, pseudo amino acid composition, composition-transition-distribution and auto-correlation function were used to map the protein sequences into numeric feature vectors. thus, development of such a computational tool will be helpful in predicting the insecticide resistant proteins, which can be targeted for developing appropriate insecticides. higher accuracies were obtained under rbf kernel than that of other kernels. based on the proposed approach, an online prediction server dirprot has also been developed for computational prediction of insecticide resistant proteins, which is accessible at http://cabgrid.res.in:8080/dirprot/. the proposed approach is believed to supplement the efforts needed to develop dynamic insecticides in wet-lab by targeting the insecticide resistant proteins. insecticide resistancesvmdi-peptide compositioncytochrome p450gabahttp://dx.doi.org/http://dx.doi.org/ /501100001503indian council of agricultural researchagril.edn.4-1/2013-a&p dated.2014issue-copyright-statement© the author 2017 resistance to different insecticides is conferred by the proteins encoded from certain class of genes of the insects. insecticide resistance is a major challenge for the control program of insect pests in the fields of crop protection, human and animal health etc.
 insecticide resistance is a major challenge for the control program of insect pests in the fields of crop protection, human and animal health etc. resistance to different insecticides is conferred by the proteins encoded from certain class of genes of the insects. to distinguish the insecticide resistant proteins from non-resistant proteins, no computational tool is available till date. thus, development of such a computational tool will be helpful in predicting the insecticide resistant proteins, which can be targeted for developing appropriate insecticides.  five different sets of feature viz., amino acid composition, di-peptide composition, pseudo amino acid composition, composition-transition-distribution and auto-correlation function were used to map the protein sequences into numeric feature vectors. the encoded numeric vectors were then used as input in support vector machine for classification of insecticide resistant and non-resistant proteins. higher accuracies were obtained under rbf kernel than that of other kernels. further, accuracies were observed to be higher for dpc feature set as compared to others. the proposed approach achieved an overall accuracy of >90% in discriminating resistant from non-resistant proteins. further, the two classes of resistant proteins i.e., detoxification-based and target-based were discriminated from non-resistant proteins with >95% accuracy. besides, >95% accuracy was also observed for discrimination of proteins involved in detoxification- and target-based resistance mechanisms. the proposed approach not only outperformed blastp, psi-blast and delta-blast algorithms, but also achieved >92% accuracy while assessed using an independent dataset of insecticide resistant proteins.  this paper presents the first computational approach for discriminating the insecticide resistant proteins from non-resistant proteins. based on the proposed approach, an online prediction server dirprot has also been developed for computational prediction of insecticide resistant proteins, which is accessible at http://cabgrid.res.in:8080/dirprot/. the proposed approach is believed to supplement the efforts needed to develop dynamic insecticides in wet-lab by targeting the insecticide resistant proteins. electronic supplementary material the online version of this article contains supplementary material, which is available to authorized users. keywords insecticide resistancesvmdi-peptide compositioncytochrome p450gabahttp://dx.doi.org/http://dx.doi.org/ /501100001503indian council of agricultural researchagril.edn.4-1/2013-a&p dated.2014issue-copyright-statement© the author 2017 