An improved sequence based prediction protocol for DNA-binding proteins using SVM and comprehensive feature analysis
© Zou et al.; licensee BioMed Central Ltd. 2013
Received: 29 September 2012
Accepted: 4 March 2013
Published: 9 March 2013
DNA-binding proteins (DNA-BPs) play a pivotal role in both eukaryotic and prokaryotic proteomes. There have been several computational methods proposed in the literature to deal with the DNA-BPs, many informative features and properties were used and proved to have significant impact on this problem. However the ultimate goal of Bioinformatics is to be able to predict the DNA-BPs directly from primary sequence.
In this work, the focus is how to transform these informative features into uniform numeric representation appropriately and improve the prediction accuracy of our SVM-based classifier for DNA-BPs. A systematic representation of some selected features known to perform well is investigated here. Firstly, four kinds of protein properties are obtained and used to describe the protein sequence. Secondly, three different feature transformation methods (OCTD, AC and SAA) are adopted to obtain numeric feature vectors from three main levels: Global, Nonlocal and Local of protein sequence and their performances are exhaustively investigated. At last, the mRMR-IFS feature selection method and ensemble learning approach are utilized to determine the best prediction model. Besides, the optimal features selected by mRMR-IFS are illustrated based on the observed results which may provide useful insights for revealing the mechanisms of protein-DNA interactions. For five-fold cross-validation over the DNAdset and DNAaset, we obtained an overall accuracy of 0.940 and 0.811, MCC of 0.881 and 0.614 respectively.
The good results suggest that it can efficiently develop an entirely sequence-based protocol that transforms and integrates informative features from different scales used by SVM to predict DNA-BPs accurately. Moreover, a novel systematic framework for sequence descriptor-based protein function prediction is proposed here.
DNA binding proteins (DNA-BPs) that interact with DNA are pivotal to the cell function such as DNA replication, transcription, packaging recombination and other fundamental activities associated with DNA. DNA-BPs represent a broad category of proteins, known to be highly diverse in sequence, structure and function. Structurally, they have been divided into eight structural/functional groups, which were further classified into 54 structural families . Functionally, protein-DNA interactions play various roles across the entire genome as previously mentioned . At present, several experimental techniques (such as filter binding assays, genetic analysis, chromatin immunoprecipitation on microarrays, and X-ray crystallography) have been used for identifying DNA-BPs. However, experimental approaches for identifying the DNA-BPs are costly and time consuming. Therefore, a reliable identification of DNA-BPs as well as DNA-binding sites with effective computational approach is an important research topic in the proteomics fields, which can play a vital role in proteome function annotation and discovery of potential therapeutics for genetic diseases and reliable diagnostics.
Computational prediction of proteins that interact with DNA is a difficult task, and state of the art methods have shown only limited success in this arena at present. Previously, there have been several machine-learning methods developed for prediction of DNA-BPs in the literature. Broadly, these methods can be divided into two categories: i) analysis from protein structure [3-6] and ii) prediction from amino acid sequence [7-12]. The accuracy of structure-based prediction methods is usually higher, but it can’t be used in high throughput annotation with the limited number of protein structures. Theoretically, the sequence of a protein contains all the necessary information to predict its function . Until now, many methods for predicting protein function directly from amino acid sequences are useful tools in the study of uncharacterized protein families and in comparative genomics . There are two major problems in the task of computational protein function prediction, which are the choice of the protein representation and the choice of the classification algorithm. To obtain good predictive model, various machine-learning algorithms such as support vector machine (SVM) [8, 10, 15-18], neural network [3, 6, 19], random forest [12, 20], naïve Bayes classifiers [21, 22], nearest neighbor  and ensemble classifiers [24, 25] have been used to build classification models. Among these, the most widely used algorithm for prediction of DNA-BPs is SVM.
In context to the current study, SVM learns the features specific to the DNA-BPs and generates support vectors decisive for possible classification of any given sequence as DNA-BPs and achieved satisfactory results. The most important challenge for SVM-based prediction is to find a suitable way to fully describe the information implied in protein-DNA interactions . There are several different protein features and feature extraction methods that can be used [8, 27-29] and a comprehensive survey of these methods can be found in related research work [30, 31]. However, the underlying principle of protein-DNA interactions is still largely unknown. It is desirable to explore the implications of those already identified features and newly undiscovered properties by machine learning approaches to further advance the prediction accuracy and understand the binding mechanism of DNA-BPs.
Thus, a systematic comparison of all protein features known to perform well is investigated in this article. We propose a novel method for predicting DNA-BPs using the SVM algorithm in conjunction with comprehensive feature analysis based on protein sequence. A recent work about mechanisms of protein folding research  has shown that the property factors of protein can be naturally clustered into two classes. One class is comprised of properties that favor sequentially localized interaction clusters; the other class is in favor of globally distributed interactions. Following the methodology introduced earlier in related protein function prediction work [29-31], we consider a feature vector (xi) to represent proteins which are derived from sequences broadly from three main levels: Global sequence descriptors, Nonlocal sequence descriptors and Local sequence descriptors. Feature vectors extracted from different sequence levels contain information about characteristics of the proteins at different scales which may be helpful in describing the information implied in protein-DNA interactions and improving the final model accuracy.
Four types of datasets are used here: i) DNAdset consists of partial sequences (binding regions or DNA binding domains) ii) DNAaset consists of full-length DNA-binding proteins. It’s reported that models trained on DNA domains or partial sequences are not suitable for predicting DNA binding proteins and vice versa, so separate methods are necessary for predicting DNA-binding domains and DNA-binding proteins . iii) An independent test set called DNAiset used for testing and comparing. iiii) DNArset with non-equal number of positives and negatives used for evaluating our method in real life.
The domain dataset also called DNAdset, consists of 231 DNA-BPs and 231 non-binding proteins with known structures which were obtained from a union of datasets used in previously related studies [8, 9, 16]. After clustering by CD-HIT  and careful inspection, these proteins have less than 40% sequence identity between each pair and without irregular amino acid characters such as “X” and “Z”. Thus the obtained DNAdset consists of 462 proteins, half of which are DNA-BPs and the other half are non-binding proteins. A complete list of all the PDB codes for DNAdset can be found in Additional file 1.
To evaluate effectiveness of our methods by comparing with previously famous studies, we used the benchmark DNAaset from reported papers [10, 17]. The dataset consists of 1153 DNA-binding proteins and 1153 non-binding proteins obtained from Swiss-Prot. No two protein sequences have similarity more than 25% and without irregular amino acid characters such as “X” and “Z”.
In order to evaluate performance of our models on dataset not used for training or testing and compare with other reported methods, we obtained an independent dataset called DNAiset from newly determined DNA-binding protein structures from PDB by keyword searching (released on 2012-01-01 and later) and non-binding proteins used in a reported prediction method . To reduce the redundancy and homology bias, CD-HIT  and PISCES  programs were used to ensure no two protein sequences have similarity more than 30% between DNAiset and two training sets (i.e. DNAdset and DNArset). Finally, the DNAiset has 80 DNA-binding protein chains selected from PDB and 192 non-binding proteins obtained from a newly developed web server named iDNA-Prot . A complete list of all the PDB codes for DNAiset can be found in Additional file 1.
Equal number of positives and negatives is important for developing a powerful predictor for a protein system. It’s also important for evaluating the performance of the prediction model where one can simply calculate the accuracy. All the above datasets in our study have equal number of DNA-binding proteins and non-binding proteins. However, in a real-world situation, DNA-binding proteins are only a fraction of all proteins. It’s one of the (relatively) new problems called imbalanced dataset which has received an increasing attention since the workshop at AAAI 2000 . This raises questions on whether models developed on equal numbers will be effective in real life. Will the method have a significantly poorer performance with more negatives in a test case? Thus, we created a more realistic dataset called DNArset to answer it. This dataset has 231 DNA-BPs used in DNAdset and 1500 non-binding proteins used by Kumar et al. as their “DNArset” .
Support vector machine
Support vector machine (SVM) is a machine learning algorithm based on statistical learning theory presented by Vapnik (1998). It takes a set of feature vectors as the input, along with their output, which is used for training of model. Application of SVM in bioinformatics to various topics has been explored . In this study, publicly available LIBSVM package version 3.11  is used for the implementation of SVM and the RBF is taken as the kernel function, the tunable parameters are optimized based on grid search method to deliver high accuracy. All feature descriptors derived below were normalized in the range of [0, 1] by using formula (value-minimum)/ (maximum-minimum) before training SVM.
To develop a powerful function predictor for a protein system, one of the keys is to formulate the datasets with an effective mathematical expression that can truly reflect their intrinsic correlation with the attribute to be predicted . To realize this, we assess four kinds of features including composition information, physicochemical property, evolutionary information and structural/functional property. The feature vector representations of these features are generated from three different levels, including Global sequence descriptors, Nonlocal descriptors and Local descriptors. Here, the composition information including overall amino acid composition (global descriptors), Dipeptide composition (nonlocal descriptors) and split amino acid composition (local descriptors). The other three kinds of properties are transformed by three different feature transformation methods which are introduced detailedly in Methods section.
Overall amino acid composition (OAAC)
Where p i represents the occurrence frequency of the i-th native amino acid in the protein, n i is the number of the i-th native amino acid in sequence, L is the length of the sequence in protein P.
Dipeptide composition (DPC)
Where D s (i, j) is the number of Dipeptide represented by amino acid type i and j with s skips, f s (i, j) represents the occurrence frequency. N is the length of the sequence.
Then, we concatenate the vector elements of DP 0 , DP 1 and DP 2 together and the mRMR method  is adopted to select the first 400 descriptors from the total of 1200 dimensions as DPC used in this paper.
Split amino acid composition (SAAC)
Split Amino Acid Composition (SAAC) was introduced where the protein sequence is divided into three parts: N terminal, C terminal and a region between them and composition of each part is calculated separately according to equations (1) and (2). Many previous literatures adopted SAAC for protein function prediction and achieved good results [40, 41]. In our SAAC method, the detailed strategy of splitting the protein sequence is illustrated in Split amino acid (SAA) Transformation section.
List of the AAIndex indices used in this paper
Normalized frequency of beta-sheet (Chou-Fasman, 1978b)
Normalized hydrophobicity scales for alpha+beta-proteins (Cid et al., 1992)
Normalized average hydrophobicity scales (Cid et al., 1992)
Number of hydrogen bond donors (Fauchere et al., 1988)
Positive charge (Fauchere et al., 1988)
Helix termination parameter at posision j+1 (Finkelstein et al., 1991)
Alpha-helix indices for alpha/beta-proteins (Geisow-Roberts, 1980)
Beta-strand indices for beta-proteins (Geisow-Roberts, 1980)
Average relative probability of beta-sheet (Kanehisa-Tsong, 1980)
Net charge (Klein et al., 1984)
Side chain interaction parameter (Krigbaum-Rubin, 1971)
Conformational preference for all beta-strands (Lifson-Sander, 1979)
Retention coefficient in HPLC, pH7.4 (Meek, 1980)
Short and medium range non-bonded energy per atom (Oobatake-Ooi, 1977)
Normalized frequency of alpha-helix in all-alpha class (Palau et al., 1981)
Weights for beta-sheet at the window position of 3 (Qian-Sejnowski, 1988)
Side chain orientational preference (Rackovsky-Scheraga, 1977)
Mean polarity (Radzicka-Wolfenden, 1988)
Side chain hydropathy, corrected for solvation (Roseman, 1988)
Optimal matching hydrophobicity (Sweet-Eisenberg, 1983)
Bulkiness (Zimmerman et al., 1968)
Isoelectric point (Zimmerman et al., 1968)
Normalized positional residue frequency at helix termini C4’ (Aurora-Rose, 1998)
Free energy in beta-strand conformation (Munoz-Serrano, 1994)
Hydropathy scale based on self-information values in the two-state model (20% accessibility) (Naderi-Manesh et al., 2001)
Hydropathy scale based on self-information values in the two-state model (36% accessibility) (Naderi-Manesh et al., 2001)
Apparent partition energies calculated from Chothia index (Guy, 1985)
Optimized relative partition energies - method C (Miyazawa-Jernigan, 1999)
To use the evolution information, the position-specific scoring matrix (PSSM)  profiles are adopted, which have been widely used in protein function prediction and other bioinformatics problems with notable improvement of performance [10, 45]. Here, the PSSM profiles are generated by using the PSI-Blast program  to search the non-redundant (NR) database (released on 14 May 2011) through three iterations with 0.001 as the E-value cutoff for multiple sequence alignment. The final PSSM scoring matrix has 20×L elements (excluding dummy residue X), where L is the length of protein.
Secondary structure composition
Secondary structure is an important structural feature of protein that can significantly improve the function prediction performance [46, 47]. In this study, secondary structure calculation is carried out by PSIPRED v3.0 , which is one of the state-of-the-art protein secondary structure prediction methods with an accuracy of up to 80%. PSIPRED predicts secondary structure for each residue in a protein and provides a confidence score for three types of secondary structures: helices, β-sheets and coil regions, thus we get 3×L feature values where L is the length of protein.
Disorder feature score
Over the past decade, there has been a growing acknowledgement that a large proportion of proteins within most proteomes contain disordered regions. Protein with disordered regions can play important functional roles. Its flexibility is advantageous to proteins that recognize multiple target molecules including biomacromolecules like DNA with high specificity and low affinity [49-51]. The IUPred  method is used to score the disorder status of each amino acid which recognizes intrinsically disordered protein regions from amino acid sequences by estimating their total pairwise interresidue interaction energy. The prediction type option for IUPred is set as long and we get L feature values where L is the length of protein.
Feature transformation method
Three different methods are used here to transform pre-selected protein features into uniform length descriptors to capture various types of information implied in proteins. The following part describes in detail the methodology for each of these different transformation methods.
Overall composition-transition-distribution (OCTD)
The original CTD method was first introduced by Dubchak et al.  as a global description of protein sequence for predicting protein folding class. Recently, CTD has been adopted by more and more leading investigators for protein function and structure studies [54, 55] . Composition (C) is the number of amino acids of a particular property divided by the total number of amino acids. Transition (T) characterizes the percent frequency with which amino acids of a particular property is followed by amino acids of a different property. Distribution (D) measures the chain length within which the first, 25, 50, 75 and 100% of the amino acids of a particular property is located respectively. Here, we develop a new variant method of CTD named Overall Composition-Transition-Distribution (OCTD). Initially, we represent the sequence numerically for a particular feature. Then, we normalize the feature values in the range of [0, 1] using formula (value-minimum)/ (maximum-minimum), amino acids were grouped into two classes according to its feature values of threshold 0.5. Finally, the CTD method is adopted to represent the amino acid properties distribution pattern of a specific property along the protein sequence.
Autocross-covariance (ACC) Transformation
Where i is one of the properties, j is the position in the sequence, L is the length of the protein sequence S i,j is the feature value of i at position j,
Thus, the number of AC variables can be calculated as P×LG, where P is the number of feature value, LG is the maximum of lg (lg=1,2,…,LG).
Split amino acid (SAA) Transformation
There have been several ways to calculate protein local features, the proposed method split amino acid composition (SAAC) where composition of N-terminal, middle and C-terminal of protein is computed separately [40, 41]. The PNPRD method which divided the protein sequence into 10 local regions based on positively and negatively charged residues  is some kind of variation from CTD method. There are also more concise method by splitting each protein into 10 local regions of varying length and CTD method was used to exact descriptors from local regions [60, 61].
Here, we adopt the powerful split amino acid (SAA) method to represent local composition of different protein features selected before. Firstly, each sequence is split into three parts: N-terminal, middle and C-terminal. The N-terminal part is further divided into four regions to consider ambiguity in the length and position of signal sequences. Then, the mean of different features corresponding to the six divided local regions are obtained to generate a fixed number of local descriptors.
Where x and y are two random variables, p(x, y) is the joint probabilistic density, p(x) and p(y) are the marginal probabilistic densities respectively. The mathematical description of the algorithm was detailedly presented in the Peng’s previous study . To calculate MI, the joint probabilistic density and the marginal probabilistic densities of the two vectors were used. A parameter t is introduced here to deal with these variables. Suppose mean to be the average value of one feature in all samples, and std to be the standard deviation, the feature of each sample would be classified into one of the three groups according to the boundaries: mean ± (t×std). Here, t was assigned to be 1.
Where f i is the i-th sorted feature in the feature list.
For each feature subset, we use SVM to construct predictor which is evaluated by five-fold cross-validation. As a result, we get a curve named IFS curve, with MCC value as its y-axis and index i of S i as its x-axis. When the overall IFS curve reaches at the peak, meanwhile, the corresponding predictor is chosen as the ultimate prediction model.
Ensemble learning method
The idea of ensemble learning methodology is to build a predictive model by integrating multiple models, treating them as a committee of decision makers. As a growing body of studies indicates that every single learning strategy has its own shortcomings and none of them could consistently perform well over all datasets. To overcome this problem, ensemble methods have been suggested as the promising measures [62, 63]. In general, an ensemble consists of a set of models and a method to combine them. We have twelve different SVM models after the first round of evaluation for exploring the performance of SVM-based modules constructed by different types of features (see Results section for details). Two popular model combination strategies: majority voting and stacking are adopted here.
In majority voting scheme, a classification of an unlabeled instance is performed according to the predicted class that obtains the highest number of votes. That is, we have twelve different classifiers in this work, if a majority of the twelve modules predict a protein as binding, then the prediction result of this protein is taken as binding. When equal number occurred, we found most proteins in this situation are non-binding so the threshold is assigned to six for binding prediction in majority voting scheme.
Stacking is a technique for achieving the highest generalization accuracy [25, 63]. Different from the majority voting scheme, stacking will learning twice. The basic idea is to learn a function that combines the predictions from the individual classifiers. Instead of using the original input attributes, it uses the predicted probabilities for every class label from the base-level classifiers as the input attributes of the target for the second-level machine learning. In our work, we use the decision values from the twelve SVM modules as the input feature vectors for final meta SVM-predictor in stacking scheme.
Indices used to evaluate the prediction method
Definition and formula
(TP + TN)/(TP + TN + FP + FN)
area under the receiving operating characteristic curve
2 · TP/(2TP + FP + FN)
TP/(TP + FN)
TN/(TN + FP)
Individual SVM-modules results
The performance of different kinds of feature descriptors
After investigating individual coding scheme, it’s confirmed that all twelve kinds of descriptors are reasonable for discriminating DNA-BPs. We intend to combine the well performed descriptors above-mentioned which mean a comprehensive presentation of protein functional features related to DNA binding. Thus, two different strategies are proposed here for achieving this goal which are ensemble learning and feature selection.
To achieve the reduction of noise and redundancy to improve the classification accuracy and the combination of more interpretable features that can help identify DNA-BPs, the proposed mRMR-IFS feature selection framework is adopted.
The mRMR program in this study is downloaded from http://penglab.janelia.org/proj/mRMR/. Using the mRMR program, we obtain the ranked mRMR list of the first 1000 features from the original 2040 descriptors for DNAdset and DNAaset separately. Within the list, a smaller index of a feature indicates that it is deemed as a more important feature in the prediction. The mRMR list is retained and will be used in the IFS procedure for feature selection and further analysis.
Ensemble learning results
The performance of feature-selection method and ensemble learning
Performance on independent dataset
In this study, we further evaluated the performance of our optimal SVM models (trained on DNAdset) on an independent dataset called DNAiset, which consists of 80 DNA-BPs and 192 non-binding proteins. The feature-selection-based model correctly predicted 72 and 170 out of 80 DNA-BPs and 192 non-binding proteins respectively, while the ensemble-based model correctly predicted 68 and 166 out of 80 DNA-BPs and 192 non-binding proteins respectively. This demonstrates that our SVM models perform equally well on the independent dataset.
Performance on realistic dataset
In a real-world situation, the number of non-binding proteins is significantly higher than DNA-BPs. Thus, it is important to build and evaluate SVM models on more realistic data rather than equal number of DNA-BPs and non-binding proteins. Hence, we obtained a realistic dataset (DNArset), which has 231 DNA-BPs and 1500 non-binding proteins. Firstly, we developed SVM model using feature selection based method on DNArset and achieved the maximum MCC of 0.720 with an accuracy of 94.16%. Then we also developed ensemble learning models and the ensemble stacking method which achieved the maximum MCC of 0.729 with an accuracy of 94.28%, while the majority voting method has a significant poorer performance with more negative proteins which demonstrates the instability of it. The detailed five fold cross validation results and mRMR-IFS results are shown in Additional File 4. Lastly, we applied models trained by DNArset with feature selection method on DNAiset which can correctly predict 60 out of 80 positives and 172 out of 192 negatives (see Additional file 4). The results further confirmed the prediction effectiveness of our method.
Comparison with existing methods
It is important to compare the performance of our protocol with existing methods in order to evaluate its effectiveness. The performance of feature-selection based models and ensemble learning models for DNA-binding prediction by fivefold cross-validation test is summarized in Table 4. The feature-selection based protocol has an accuracy of 0.940 and 0.789 for DNAdset and DNAaset respectively, accordingly ensemble learning protocol returns a little lower accuracy of 0.907 for DNAdset but a higher accuracy of 0.811 for DNAaset. For DNAaset, the performances of previously reported studies developed from it are , with accuracy of 0.742, MCC of 0.49, Sen of 0.735 and Sp of 0.749 and , with accuracy of 0.755, MCC of 0.51, Sen of 0.820 and Sp of 0.690. As shown in Table 4, the best performance in our protocol with an accuracy of 0.811, MCC of 0.614, Sen of 0.814 and Sp of 0.799 is much better.
Comparison of the predicted results by our method and some web-servers on DNAiset
In this work, we investigate the idea of ensemble of informative features from different levels for predicting DNA-BPs which is motivated by a recently research result that amino acid physical properties can fall into distinct levels . The overall protocol is aimed at representing the four important kinds of properties of protein appropriately by different transformation methods and seeking the optimal feature descriptors for presentation of DNA-BPs. The performances of individual modules indicate the usefulness of features from various levels and their dissimilarity. Based on the obtained different kinds of feature descriptors, we take two strategies for the construction of the final prediction models which are mRMR-IFS feature selection protocol and ensemble learning approach. Encouragingly, we get good performance of Acc of 0.940 for DNAdset with the mRMR-IFS method and Acc of 0.811 for DNAaset with ensemble learning approach, and the performance on independent test set is also good.
Our experiments indicate that it may be helpful to develop a successful machine method to predict the DNA-BPs by exploiting protein sequence comprehensively. However more explorations about amino acid properties are still needed in this direction and further work on interpreting these features and exploring mechanisms of protein-DNA interactions are underway.
Correspondence and requests for reprints to:
Prof. Honglin Li
School of Pharmacy, East China University of Science and Technology 130 Mei Long Road, Shanghai 200237
This work was supported by the Fundamental Research Funds for the Central Universities, the National Natural Science Foundation of China (grants 21173076, 81102375, 81230090, 81222046 and 81230076), the Special Fund for Major State Basic Research Project (grant 2009CB918501), the Shanghai Committee of Science and Technology (grant 11DZ2260600), and the 863 Hi-Tech Program of China (grant 2012AA020308). Honglin Li is also sponsored by Program for New Century Excellent Talents in University (grant NCET-10-0378).
- Luscombe NM, Austin SE, Berman HM, Thornton JM: An overview of the structures of protein-DNA complexes. Genome Biol 2000,1(1):1-37.View Article
- Ren B, Robert F, Wyrick JJ, Aparicio O, Jennings EG, Simon I, Zeitlinger J, Schreiber J, Hannett N, Kanin E: Genome-wide location and function of DNA binding proteins. Science 2000,290(5500):2306-2309.View ArticlePubMed
- Ahmad S, Sarai A: Moment-based prediction of DNA-binding proteins. J Mol Biol 2004,341(1):65-71.View ArticlePubMed
- Zhao H, Yang Y, Zhou Y: Structure-based prediction of DNA-binding proteins by structural alignment and a volume-fraction corrected DFIRE-based energy function. Bioinformatics 2010,26(15):1857-1863.PubMed CentralView ArticlePubMed
- Tjong H, Zhou HX: DISPLAR: an accurate method for predicting DNA-binding sites on protein surfaces. Nucleic Acids Res 2007,35(5):1465-1477.PubMed CentralView ArticlePubMed
- Stawiski EW, Gregoret LM, Mandel-Gutfreund Y: Annotating nucleic acid-binding function based on protein structure. J Mol Biol 2003,326(4):1065-1079.View ArticlePubMed
- Cai YD, Lin SL: Support vector machines for predicting rRNA-, RNA-, and DNA-binding proteins from amino acid sequence. Biochim Biophys Acta 2003,1648(1-2):127-133.View ArticlePubMed
- Fang Y, Guo Y, Feng Y, Li M: Predicting DNA-binding proteins: approached from Chou’s pseudo amino acid composition and other specific sequence features. Amino Acids 2008,34(1):103-109.View ArticlePubMed
- Gao M, Skolnick J: A threading-based method for the prediction of DNA-binding proteins with application to the human genome. PLoS Comput Biol 2009,5(11):e1000567.PubMed CentralView ArticlePubMed
- Kumar M, Gromiha M, Raghava G: Identification of DNA-binding proteins using support vector machines and evolutionary profiles. BMC Bioinforma 2007,8(1):463.View Article
- Shao X, Tian Y, Wu L, Wang Y, Jing L, Deng N: Predicting DNA- and RNA-binding proteins from sequences with kernel methods. J Theor Biol 2009,258(2):289-293.View ArticlePubMed
- Lin WZ, Fang JA, Xiao X, Chou KC: IDNA-prot: identification of DNA binding proteins using random forest with grey model. PLoS One 2011,6(9):e24756.PubMed CentralView ArticlePubMed
- Cai YD, Doig AJ: Prediction of Saccharomyces cerevisiae protein functional class from functional domain composition. Bioinformatics 2004,20(8):1292-1300.View ArticlePubMed
- Brameier M, Haan J, Krings A, MacCallum R: Automatic discovery of cross-family sequence features associated with protein function. BMC Bioinforma 2006,7(1):16.View Article
- Brown J, Akutsu T: Identification of novel DNA repair proteins via primary sequence, secondary structure, and homology. BMC Bioinforma 2009,10(1):25.View Article
- Bhardwaj N, Langlois RE, Zhao G, Lu H: Kernel-based machine learning protocol for predicting DNA-binding proteins. Nucleic Acids Res 2005,33(20):6486-6493.PubMed CentralView ArticlePubMed
- Huang HL, Lin IC, Liou YF, Tsai CT, Hsu KT, Huang WL, Ho SJ, Ho SY: Predicting and analyzing DNA-binding domains using a systematic approach to identifying a set of informative physicochemical and biochemical properties. BMC Bioinforma 2011,12(Suppl 1):S47.View Article
- Xiong Y, Liu J, Wei DQ: An accurate feature-based method for identifying DNA-binding residues on protein surfaces. Proteins 2011,79(2):509-517.View ArticlePubMed
- Ahmad S, Andrabi M, Mizuguchi K, Sarai A: Prediction of mono- and di-nucleotide-specific DNA-binding sites in proteins using neural networks. BMC Struct Biol 2009, 9: 30.PubMed CentralView ArticlePubMed
- Nimrod G, Schushan M, Szilágyi A, Leslie C, Ben-Tal N: iDBPs: a web server for the identification of DNA binding proteins. Bioinformatics 2010,26(5):692-693.PubMed CentralView ArticlePubMed
- Yan C, Terribilini M, Wu F, Jernigan R, Dobbs D, Honavar V: Predicting DNA-binding sites of proteins from amino acid sequence. BMC Bioinforma 2006,7(1):262.View Article
- Govindan G, Nair AS: New Feature Vector for Apoptosis Protein Subcellular Localization Prediction. In Advances in Computing and Communications Communications. Volume 190 . Edited by: Abraham A. Kochi: Springer Berlin Heidelberg; 2011:294-301.
- Qian ZL, Cai YD, Li YX: A novel computational method to predict transcription factor DNA binding preference. Biochem Biophys Res Commun 2006,348(3):1034-1037.View ArticlePubMed
- Nanni L, Lumini A: Combing ontologies and dipeptide composition for predicting DNA-binding proteins. Amino Acids 2008,34(4):635-641.View ArticlePubMed
- Xia JF, Zhao XM, Huang DS: Predicting protein-protein interactions from protein sequences using meta predictor. Amino Acids 2010,39(5):1595-1599.View ArticlePubMed
- Liu ZP, Wu LY, Wang Y, Zhang XS, Chen LN: Bridging protein local structures and protein functions. Amino Acids 2008,35(3):627-650.View ArticlePubMed
- Chou KC: Some remarks on protein attribute prediction and pseudo amino acid composition. J Theor Biol 2011,273(1):236-247.View ArticlePubMed
- Yuan Y, Shi X, Li X, Lu W, Cai Y, Gu L, Liu L, Li M, Kong X, Xing M: Prediction of interactiveness of proteins and nucleic acids based on feature selections. Mol Divers 2010,14(4):627-633.View ArticlePubMed
- Song J, Tan H, Takemoto K, Akutsu T: HSEpred: predict half-sphere exposure from protein sequences. Bioinformatics 2008,24(13):1489-1497.View ArticlePubMed
- Nanni L, Brahnam S, Lumini A: High performance set of PseAAC and sequence based descriptors for protein classification. J Theor Biol 2010,266(1):1-10.View ArticlePubMed
- Zhang Z, Kochhar S, Grigorov MG: Descriptor-based protein remote homology identification. Protein Sci 2005,14(2):431-444.PubMed CentralView ArticlePubMed
- Rackovsky S: Global characteristics of protein sequences and their implications. Proc Natl Acad Sci USA 2010,107(19):8623-8626.PubMed CentralView ArticlePubMed
- Huang Y, Niu B, Gao Y, Fu L, Li W: CD-HIT Suite: a web server for clustering and comparing biological sequences. Bioinformatics 2010,26(5):680-682.PubMed CentralView ArticlePubMed
- Wang G, Dunbrack RL Jr: PISCES: a protein sequence culling server. Bioinformatics 2003,19(12):1589-1591.View ArticlePubMed
- Chawla NV, Japkowicz N, Kotcz A: Editorial: special issue on learning from imbalanced data sets. SIGKDD Explor Newsl 2004,6(1):1-6.View Article
- Chang CC, Lin CJ: LIBSVM: A library for support vector machines. ACM Transact Intell Syst Technol 2011,2(3):27.
- Feng ZP: Prediction of the subcellular location of prokaryotic proteins based on a new representation of the amino acid composition. Biopolymers 2001,58(5):491-499.View ArticlePubMed
- Bhasin M, Raghava GPS: ESLpred: SVM-based method for subcellular localization of eukaryotic proteins using dipeptide composition and PSI-BLAST. Nucleic Acids Res 2004,32(suppl 2):W414-W419.PubMed CentralView ArticlePubMed
- Peng H, Long F, Ding C: Feature selection based on mutual information: criteria of max-dependency, max-relevance, and min-redundancy. IEEE Trans Pattern Anal Mach Intell 2005,27(8):1226-1238.View ArticlePubMed
- Tantoso E, Li KB: AAIndexLoc: predicting subcellular localization of proteins based on a new representation of sequences using amino acid indices. Amino Acids 2008,35(2):345-353.View ArticlePubMed
- Afridi T, Khan A, Lee Y: Mito-GSAAC: mitochondria prediction using genetic ensemble classifier and split amino acid composition. Amino Acids 2012,42(4):1443-1454.View ArticlePubMed
- Han P, Zhang X, Feng Z-P: Predicting disordered regions in proteins using the profiles of amino acid indices. BMC Bioinforma 2009,10(Suppl 1):S42.View Article
- Kawashima S, Pokarowski P, Pokarowska M, Kolinski A, Katayama T, Kanehisa M: AAindex: amino acid index database, progress report 2008. Nucleic Acids Res 2008,36(Database issue):D202-205.PubMed CentralPubMed
- Schaffer AA, Aravind L, Madden TL, Shavirin S, Spouge JL, Wolf YI, Koonin EV, Altschul SF: Improving the accuracy of PSI-BLAST protein database searches with composition-based statistics and other refinements. Nucleic Acids Res 2001,29(14):2994-3005.PubMed CentralView ArticlePubMed
- Chen SA, Ou YY, Lee TY, Gromiha MM: Prediction of transporter targets using efficient RBF networks with PSSM profiles and biochemical properties. Bioinformatics 2011,27(15):2062-2067.View ArticlePubMed
- Song J, Tan H, Wang M, Webb GI, Akutsu T: TANGLE: Two-level support vector regression approach for protein backbone torsion angle prediction from primary sequences. PLoS One 2012,7(2):e30361.PubMed CentralView ArticlePubMed
- Chu WY, Huang YF, Huang CC, Cheng YS, Huang CK, Oyang YJ: ProteDNA: a sequence-based predictor of sequence-specific DNA-binding residues in transcription factors. Nucleic Acids Res 2009,37(suppl 2):W396-W401.PubMed CentralView ArticlePubMed
- Jones DT: Protein secondary structure prediction based on position-specific scoring matrices. J Mol Biol 1999,292(2):195-202.View ArticlePubMed
- Wright PE, Dyson HJ: Intrinsically unstructured proteins: re-assessing the protein structure-function paradigm. J Mol Biol 1999,293(2):321-331.View ArticlePubMed
- Lobley A, Swindells MB, Orengo CA, Jones DT: Inferring function using patterns of native disorder in proteins. PLoS Comput Biol 2007,3(8):e162.PubMed CentralView ArticlePubMed
- Weiss MA, Ellenberger T, Wobbe CR, Lee JP, Harrison SC, Struhl K: Folding transition in the DNA-binding domain of GCN4 on specific binding to DNA. Nature 1990,347(6293):575-578.View ArticlePubMed
- Dosztányi Z, Csizmok V, Tompa P, Simon I: IUPred: web server for the prediction of intrinsically unstructured regions of proteins based on estimated energy content. Bioinformatics 2005,21(16):3433-3434.View ArticlePubMed
- Dubchak I, Muchnik I, Holbrook SR, Kim SH: Prediction of protein folding class using global description of amino acid sequence. Proc Natl Acad Sci USA 1995,92(19):8700-8704.PubMed CentralView ArticlePubMed
- Govindan G, Nair AS: Composition, Transition and Distribution (CTD) - A dynamic feature for predictions based on hierarchical structure of cellular sorting . Hyderabad: India Conference (INDICON); 2011. 2011 Annual IEEE; 16-18 DecView Article
- Cai CZ, Han LY, Ji ZL, Chen X, Chen YZ: SVM-Prot: Web-based support vector machine software for functional classification of a protein from its primary sequence. Nucleic Acids Res 2003,31(13):3692-3697.PubMed CentralView ArticlePubMed
- Wold S, Jonsson J, Sjörström M, Sandberg M, Rännar S: DNA and peptide sequences and chemical processes multivariately modelled by principal component analysis and partial least-squares projections to latent structures. Anal Chim Acta 1993,277(2):239-253.View Article
- Doytchinova IA, Flower DR: VaxiJen: a server for prediction of protective antigens, tumour antigens and subunit vaccines. BMC Bioinforma 2007, 8: 4.View Article
- Guo Y, Yu L, Wen Z, Li M: Using support vector machine combined with auto covariance to predict protein-protein interactions from protein sequences. Nucleic Acids Res 2008,36(9):3025-3030.PubMed CentralView ArticlePubMed
- Lee BJ, Shin MS, Oh YJ, Oh HS, Ryu KH: Identification of protein functions using a machine-learning approach based on sequence-derived properties. Proteome science 2009, 7: 27.PubMed CentralView ArticlePubMed
- Yang L, Xia JF, Gui J: Prediction of protein-protein interactions from protein sequence using local descriptors. Protein Pept Lett 2010,17(9):1085-1090.View ArticlePubMed
- Davies MN, Secker A, Freitas AA, Clark E, Timmis J, Flower DR: Optimizing amino acid groupings for GPCR classification. Bioinformatics 2008,24(18):1980-1986.View ArticlePubMed
- Si J, Zhang Z, Lin B, Schroeder M, Huang B: MetaDBSite: a meta approach to improve protein DNA-binding sites prediction. BMC Syst Biol 2011,5(Suppl 1):S7.PubMed CentralView ArticlePubMed
- Rokach L: Ensemble-based classifiers. Artif Intell Rev 2010,33(1):1-39.View Article
- Sathyapriya R, Vijayabaskar MS, Vishveshwara S: Insights into Protein-DNA Interactions through structure network analysis. PLoS Comput Biol 2008,4(9):e1000170.PubMed CentralView ArticlePubMed
- Szilagyi A, Skolnick J: Efficient prediction of nucleic acid binding function from low-resolution protein structures. J Mol Biol 2006,358(3):922-933.View ArticlePubMed
- Ghosh S, Marintcheva B, Takahashi M, Richardson CC: C-terminal phenylalanine of bacteriophage T7 single-stranded DNA-binding protein is essential for strand displacement synthesis by T7 DNA polymerase at a nick in DNA. J Biol Chem 2009,284(44):30339-30349.PubMed CentralView ArticlePubMed
- Rohs R, West SM, Sosinsky A, Liu P, Mann RS, Honig B: The role of DNA shape in protein-DNA recognition. Nature 2009,461(7268):1248-1253.PubMed CentralView ArticlePubMed
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.