Skip to main content

Prediction of the functional class of metal-binding proteins from sequence derived physicochemical properties by support vector machine approach


Metal-binding proteins play important roles in structural stability, signaling, regulation, transport, immune response, metabolism control, and metal homeostasis. Because of their functional and sequence diversity, it is desirable to explore additional methods for predicting metal-binding proteins irrespective of sequence similarity. This work explores support vector machines (SVM) as such a method. SVM prediction systems were developed by using 53,333 metal-binding and 147,347 non-metal-binding proteins, and evaluated by an independent set of 31,448 metal-binding and 79,051 non-metal-binding proteins. The computed prediction accuracy is 86.3%, 81.6%, 83.5%, 94.0%, 81.2%, 85.4%, 77.6%, 90.4%, 90.9%, 74.9% and 78.1% for calcium-binding, cobalt-binding, copper-binding, iron-binding, magnesium-binding, manganese-binding, nickel-binding, potassium-binding, sodium-binding, zinc-binding, and all metal-binding proteins respectively. The accuracy for the non-member proteins of each class is 88.2%, 99.9%, 98.1%, 91.4%, 87.9%, 94.5%, 99.2%, 99.9%, 99.9%, 98.0%, and 88.0% respectively. Comparable accuracies were obtained by using a different SVM kernel function. Our method predicts 67% of the 87 metal-binding proteins non-homologous to any protein in the Swissprot database and 85.3% of the 333 proteins of known metal-binding domains as metal-binding. These suggest the usefulness of SVM for facilitating the prediction of metal-binding proteins. Our software can be accessed at the SVMProt server


Metal-binding proteins play important roles in structural stability and complex formation [15], gene expression regulation and alteration [1, 69], DNA processing [10], signaling processes and cellular events [11], transport [8, 12, 13], metabolism control [1, 4, 12, 14], metal homeostasis [15, 16], antibody recognition [17], and other events such as cellular respiration, muscle movement, and antioxidant defense [18]. Approximately 1/3 of structurally-determined proteins are metal-bound[19], and large percentages of metals present in human body are bound to proteins [4, 20]. Identification of metal-binding proteins and knowledge of metal-protein interactions is important for elucidating the function and functional mechanism of proteins and biological processes.

Metal-binding proteins have been identified by such experimental approaches as absorbance spectroscopy[21], gel electrophoresis [22], metal-affinity columns and shift assay [23], chromatography [16], mass spectroscopy [22], NMR [6], and combined spectroscopic studies [24]. However, some of these methods generally require a purified or semi-purified target of interest, do not facilitate identification of unknown targets form complex protein mixtures, or require multi-step processes and very specialized equipment, which limit their application ranges [23]. Therefore, there is a need to explore other methods including computational approaches for facilitating the identification of metal-biding proteins to complement these experimental methods.

Several computational methods have been explored for identifying and characterizing metal-binding proteins. In many metalloproteins, the metal ions tightly bind to the proteins and their metal-bound structures could be accurately determined by x-ray crystallography [3, 5, 10, 17]. Thus structural information has been used for predicting metal-binding sites based on the detection of principal liganding residues and metal-ligand complex architectures [25, 26], the use of common local structural parameters [25], combination of sequence and structural profiles [27], analysis of bond strength contributions [28], and the computation of force fields [29, 30]. But for those proteins with loosely or temporarily bound metals, such as enzymes that use metal ions as cofactors, the specific metal binding sites are often poorly characterized or unknown [6]. Therefore, sequence-based computational methods appear to be useful for these types of proteins and those without 3D structures. Apart from sequence similarity methods, the recently explored sequence-based methods include metal-binding sites sequence motifs [31, 32], multiple sequence alignments against known metal-binding proteins[26], and neural networks of sequence segments of amino acids of higher metal-binding propensity[33]. Moreover, combinatorial use of multiple structural, sequence alignments and annotation methods has been found to be highly useful for improving prediction accuracy of metal-binding proteins[26].

Because of the sequence, structural and functional diversity of metal-binding proteins [15, 8, 9, 1115, 17], it is desirable to explore additional methods that predict metal-binding proteins directly from sequence or sequence-derived properties. This work explored a statistical learning method, support vector machines (SVM), as such an approach. SVM has been successfully used for predicting the functional classes of molecule-binding proteins such as RNA-binding proteins [34, 35], DNA-binding proteins [35], lipid-binding proteins [36], and transporters [37] from sequence-derived structural and physicochemical properties and irrespective of sequence similarity. Metal-binding proteins involve a substantially more diverse spectrum of proteins than most of the other classes of proteins. For instance, the zinc-binding proteins of 16,072 sequence entries belong to 765 Pfam [38] domain families, while the EC2.7 enzymes and RNA-binding proteins of similar number of sequence entries (14,171 and 14,208) belong to 548 and 378 Pfam families respectively. The diverse spectrum of proteins poses a more critical test for constructing a SVM prediction system.

Metal-binding proteins are diverse in sequence, structure, and function[15, 8, 9, 1115, 17]. Nonetheless, metal cations generally bind to centers of high hydrophilicity and reduce the enthalpy of a system upon binding [30, 39], and metal ions bind to a shell of polar hydrophilic residues surrounded by a shell of non-polar residues [25]. The binding sites of some metal-ligand complexes have specific structural architectures [25]. To some extent, these metal binding features are similar to those of other molecule-binding features of proteins such as RNA-binding proteins, DNA-binding proteins and transporters that are also diverse in sequence, structure and function whose binding capability are mediated by certain structural and physiochemical characteristics [36, 40, 41]. Therefore, it is expected that SVM is also applicable to the prediction of metal-binding proteins.

In this paper, we developed SVM prediction systems for 10 metal-binding classes and for all metal-binding proteins. These classes are calcium-binding, cobalt-binding, copper-binding, iron-binding, magnesium-binding, manganese-binding, nickel-binding, potassium-binding, sodium-binding and zinc-binding. In addition to the estimate of the prediction accuracy by using an independent set of proteins, the performance of our developed SVM prediction systems was further evaluated by four additional tests to determine the usefulness of SVM for predicting novel metal-binding proteins and the applicability of other kernel functions. One is the evaluation of the prediction accuracies when homologous proteins are considered as one. The second is the prediction of metal-binding proteins non-homologous to any protein in Swissprot database[42]. The third is to study whether the known metal-binding domains can be predicted as metal-binding by our SVM systems. The fourth is to study the performance of SVM with a different kernel function.

Results and discussion

Overall prediction accuracy

The statistics of the datasets and prediction results of specific metal-binding classes and all metal-binding proteins are given in Table 1. In this Table, TP, FN, TN FP, SE, and SP stand for true positive (correctly predicted metal-binding proteins of specific class), false negative (specific class of metal-binding proteins incorrectly predicted as non-class-members), true negative (correctly predicted non-class-members), false positive (non-class-members incorrectly predicted as specific class of metal-binding proteins), the predicted sensitivity (accuracy for members in each metal-binding class), and the predicted specificity (accuracy for non-members of each metal-binding class). The SE for the class of calcium-binding, cobalt-binding, copper-binding, iron-binding, magnesium-binding, manganese-binding, nickel-binding, potassium-binding, sodium-binding, zinc-binding, and all metal-binding proteins is 86.3%, 81.6%, 83.5%, 94.0%, 81.2%, 85.4%, 77.6%, 90.4%, 90.9%, 74.9% and 78.1% respectively. The corresponding SP is 88.2%, 99.9%, 98.1%, 91.4%, 87.9%, 94.5%, 99.2%, 99.9%, 99.9%, 98.0%, and 88.0% respectively.

Table 1 Statistics of the datasets and prediction accuracy of individual class of metal-binding proteins and that of all metal-binding proteins. The predicted results are given in TP (true positive), FN (false negative), TN (true negative), FP (false positive), sensitivity SE = TP/(TP + FN) (accuracy for class members), specificity SP = TN/(TN + FP) (accuracy for non-members), and overall accuracy Q = (TN + TP)/(TP + FN + TN + FP). The number of members and non-members in the testing and independent evaluation sets is TP + FN or TN + FP respectively

A direct comparison with results from previous metal-binding protein prediction studies may not be most appropriate because of the differences in the protein classes predicted, datasets, protein descriptors, prediction methods and parameters. Nonetheless, a tentative comparison may provide some crude estimate regarding the level of accuracy of our method with respect to those achieved by other studies of metal-binding proteins. The reported SEs of other studies are in the range of 87%~93% for calcium-binding proteins [28] and 90~97% for all metal-binding proteins [25, 27, 29, 33]. Thus the corresponding SEs of 93.8% and 78.9% of our SVM prediction systems are comparable to those of other studies despite of the use of a significantly higher number, and thus more diverse range, of proteins in our studies.

The prediction accuracy of the non-members of each metal-binding class appears to be better than that of the members. The higher prediction accuracy for non-members likely results from the availability of more diverse set of non-members than that of members, which enables SVM to perform a better statistical learning for recognition of non-members. Based on the statistics provided on the webpage of Pfam database [38], there are over 8,000 families of proteins, from which one can generate a diverse set of non-members for each metal-binding class. Because of the differences in the number of members and that of non-members in each class, there is an imbalance between each dataset. SVM based on an imbalanced datasets tends to produce feature vectors that push the hyperplane towards the side with smaller number of data [43], which can lead to a reduced accuracy for the set either with a smaller number of samples or of less diversity. This might partly explain why the prediction accuracy for members is generally lower than that for non-members. It is however inappropriate to simply reduce the size of non-members to artificially match that of members, since this compromises the diversity needed to fully represent all non-members. Computational methods for re-adjusting biased shift of hyperplane are being explored [44]. Application of these methods may help improving SVM prediction accuracy in this and other cases involving unbalanced data.

Prediction of novel metal-binding proteins

One particular application of SVM is the prediction of novel metal-binding proteins that are non-homologous to other proteins [45]. To test this capability, Swiss-Prot database [42] was searched for metal-binding proteins having no single homologous protein in the database based on PSI-BLAST [46] results. A similarity E-value threshold of 0.1 was used for homolog search to ensure maximum exclusion of proteins that have a homolog. Those proteins found in the SVM training sets are then removed. As shown in Table 2, a total of 87 proteins are found from this process, 58 or 66.7% of these proteins were correctly predicted as metal-binding by our SVM classification systems respectively. Therefore, our SVM classification systems appear to show reasonably good capability for predicting novel metal-binding proteins based on the set of proteins tested. 9 of the 29 incorrectly predicted novel metal-binding proteins are nucleic acid binders. These include endonucleases Cfr10I, CviJI, EcoRV, PvuII and BslI, transcription activators chrR and rep2, meiosis-specific protein HOP1, protein suppressor of variegation 3–7. One possible reason for the misclassification of these nucleic acid binding proteins is that spatial conservation rather than sequence conservation plays the dominant role in the formation of active site where metal ions are clustered [47], which is more difficult to predict than the classes of proteins with more apparent sequence signatures.

Table 2 Prediction results of novel metal-binding proteins by SVMProt, where "+" represents proteins correctly predicted as metal-binding proteins, and "-" represents proteins incorrectly predicted as non-metal-binding proteins

Prediction of metal-binding proteins with specific structural characteristics

A number of metal-binding proteins contain metal-binding domains [31, 48, 49] or motifs [31, 32]. Several families of such metal-binding proteins have been documented, and examples of these families are zinc finger family[50], EF hand family[51], and Fer4 family[52]. These families have distinguished structural features responsible for metal-binding. Thus the performance of SVM prediction of metal-binding proteins can be evaluated by examining whether or not proteins containing one of these domains or motifs can be correctly classified as metal-binding.

A search of protein family and sequence databases showed that there are 2462, 780, and 534 metal-binding protein sequences known to contain zinc finger, EF hand, and Fer4 domain respectively. The majority of these sequences are included in the training and testing set of all metal-binding proteins. In the corresponding independent evaluation set, there are 890, 215, and 192 sequences containing zinc finger, EF hand, and Fer4 domain respectively. Most of these protein sequences were correctly classified as metal-binding by SVM. There are only 17, 8, and 6 misclassified sequences in the zinc finger, EF hand, and Fer4 domain families respectively. Thus our results showed the capability of our SVM prediction systems for recognizing these metal-binding proteins.

Prediction performance for metal-binding domains

Some metal-binding proteins are known to contain multiple domains that include a metal-binding domain plus one or more domains characterized by DNA binding, protein-protein interaction and other motifs [5356]. Our SVM prediction systems were trained by using physicochemical properties derived from the entire protein sequence. There is a need to evaluate how the inclusion of all the other "extra" domains may affect the prediction performance of our SVM systems. For such a purpose, our SVM systems were tested to determine to which extent they can predict known metal-binding domains as metal-binding without having to include representatives of these domains in our training sets. Metal-binding domains are searched from the Pfam database [38] by using key word of ten biological common metals- calcium, cobalt, copper, magnesium, manganese, nickel, potassium, sodium and zinc, followed by manual evaluation of the hits to select those with such annotations as experimental proved metal contained structure, metal chelating site, metal-based functional site, as well as metal channel, transporter, and carrier. A total of 333 distinct metal-binding domains are selected from this process, which include 127 domains in multi-domain metal-binding proteins. It is found that 85.3% and 81.1% of these are predicted as metal-binding. Moreover, 82.5% of the 1368 multi-domain metal-binding proteins in our independent set is correctly predicted. Hence, the inclusion of "extra" domains appears to have a limited effect on the performance of our developed SVM systems, which show certain level of capability for predicting metal-binding domains as well as metal-binding proteins.

SVM prediction performance by using a different kernel function

Apart from the Gaussian kernel function of sequence-derived physicochemical properties used in this work, several other kernel functions have been developed and applied for SVM analysis of proteins and DNAs [5759]. It is of interest to test the usefulness of some of these kernel functions for predicting metal-binding proteins. The string-kernel function has been extensively used and it has shown promising potential for protein and DNA studies [57, 58]. This kernel function is constructed by comparison of sequences of classes of proteins or DNAs and the assignment of individual weights to amino acids or nucleotides to describe physicochemical or other characteristics of the proteins and DNAs. In this work, this kernel function is used to develop three SVM systems for predicting the class of nickel-binding proteins, potassium-binding proteins, and sodium-binding proteins. Spectrum kernel with mismatches[60] is used to generate the string-kernel for each protein. Testing results by using the independent set of proteins for each class show that the SE is 75.1%, 89.5% and 88.7%, and the SP is 99.0%, 98.7% and 97.8% for each of these classes respectively. Thus comparable prediction performance can be achieved by using string-kernel SVM, which suggests the usefulness of this and other kernel functions for SVM prediction of metal-binding proteins.

Comparison of SVM prediction performance with those of other methods

To compare the prediction performance of SVM with those of other methods, our SVM classification system for predicting zinc-binding proteins was used to scan the human genome, and the predicted zinc-binding proteins were compared with those predicted by one or combinations of three other methods[26]. These methods have been used for searching potential zinc-binding proteins in human genome by means of (i) zinc-binding pattern identification via structural comparison with all available X-ray structural data, (ii) multiple sequence alignment based on libraries of zinc-binding domains, and (iii) analysis of sequence annotations[26]. SVM predicted a total of 4,518 zinc-binding proteins compared to that of 3,207 by at least one of the other three methods, 2,770 of which are mutually predicted. The percentages of mutually predicted proteins are significantly higher for those proteins predicted by using combinations of the other three methods. The numbers of proteins predicted by at least two and three of the other methods are 2,430 and 1,684 respectively, 2,256 and 1,615 of which are also predicted by SVM. Therefore, SVM is capable of predicting most of the zinc-binding proteins predicted by the combinations of the other three methods. SVM appears to predict a higher number of zinc-binding proteins than each of the other three methods. Apart from the expected prediction error, the reported problems of the other three methods associated with structurally uncharacterized, non-conserved, unclearly annotated zinc-binding proteins[26] may also contribute to the discrepancy between SVM and the other methods. For example, two SVM predicted proteins that are not predicted by the other three methods, forkhead box protein P1 and TRAF-interacting protein, are annotated as zinc-binding in GO and described to contain ZINC_FINGER_C2H2_1 in PROSITE.

Contribution of feature properties to the classification of metal-binding proteins

In this work, a total of nine feature properties were used to describe physicochemical characteristics of each protein, which have been routinely used for the prediction of other molecular-binding proteins [61]. It has been reported that, not all feature vectors contribute equally to the classification of proteins, some have been found to play relatively more prominent role than others in specific aspects of proteins [62]. It is therefore of interest to examine which feature properties play more prominent role in classification of metal-binding proteins.

In an earlier study, contribution of individual feature property to protein classification is investigated by separately conducting classification using each feature property [36, 40, 41]. The same method was employed here. An analysis on the classification of the group of all metal-binding proteins seems to suggest that, in order of prominence, the hydrophobicity, solvent accessibility, polarity and composition play more prominent role than other feature properties. Hydrophobicity have been shown to be important for metal-protein interactions such that metal binding sites usually appear in clusters with hydrophobic environment. High-affinity metal binding sites in some proteins are located at sequence segments with specific amino acid composition[63], and specific sequence motifs have been used for predicting metal-binding proteins[64, 65]. It was also found that polarity and solvent accessibility of the binding site influences the functional properties of metal-binding proteins[66]. Therefore, our prediction results are consistent with these experimental findings.


SVM appears to be a potentially useful tool for the prediction of metal-binding proteins of different classes. The prediction accuracy may be further enhanced with the further expansion of our knowledge about metal-binding proteins particularly for those small metal-binding classes, more refined representation of the structural and physicochemical properties of proteins, and the improvement of prediction algorithms such as the better treatment of imbalanced dataset. The SVM-derived metal-binding protein classification systems developed in this work can be assessed, free of charge for academic use, at the SVMProt server [67].


Selection of metal-binding and non-metal-binding proteins

All metal-binding proteins used in this study are collected from a comprehensive search of Swiss-Prot database [42]. A total of 33295 metal-binding protein sequences were obtained. Most of these proteins can be classified into one of the 10 metal-binding classes, and the number of proteins is 5426, 1467, 1645, 9462, 9688, 4214, 705, 1240, 1567 and 16072 in calcium-binding, cobalt-binding, copper-binding, iron-binding, magnesium-binding, manganese-binding, nickel-binding, potassium-binding, sodium-binding and zinc-binding class respectively. Some proteins were found to belong to more than one class. The distribution of all these proteins in different kingdoms and in top 10 host species is given in Table 3, and that of the four largest classes of metal-binding proteins is given in Table 4. From these two Tables one finds that these proteins are from diverse range of species and all species appear to be fairly adequately represented.

Table 3 Distribution of metal-binding proteins in different kingdoms and in top 10 host species of each kingdom. Not all protein sequences studied in this work are included because the host species information of some protein sequences is not yet available in the protein sequence database
Table 4 Distribution of different classes of metal-binding proteins (calcium-binding, magnesium-binding, iron-binding, and zinc-binding) in different kingdoms and in top 10 host species. Not all protein sequences studied in this work are included because the host species information of some protein sequences is not yet available in the protein sequence database

All distinct members in each class were used to construct a positive dataset for the corresponding SVM prediction system. A negative dataset, representing non-class members, are selected by a well-established procedure [45, 68, 69] such that all proteins are grouped into domain families [38] and the representative proteins of those families un-related to the specific metal-binding class are used as negative samples. Members in the other metal-binding classes were included in the negative dataset if they are not a member of the class being studied. These datasets are divided into separate training, testing and independent evaluation sets by the following procedure: First, proteins were clustered into groups based on their distance in the structural and physicochemical feature-space by using the hierarchical clustering method. An upper-limit of the largest separation of 20 was used to for each cluster. One representative protein was randomly selected from each group to form a training set that is sufficiently diverse and broadly distributed in the feature space. One or up to 50% of the remaining proteins in each group were randomly selected to form the testing set. The selected proteins from each group were further checked to ensure that they are distinguished from proteins from other groups. The remaining proteins were used as the independent evaluation set, which are also of reasonable level of diversity. Moreover, an analysis of the "similarity" proteins in each cluster shows that the majority of the proteins in a cluster are non-homologous. Thus, the testing and independent evaluation sets are expected to have certain level of usefulness for performing their task of fine-tuning the parameter of a SVM classification system and for evaluating its prediction performance. The statistics of the members and non-members of the datasets of each metal-binding class is given in Table 1.

Derivation of structural and physicochemical properties from protein sequence

Construction of the feature vector for each protein is based on the formula used in the prediction of RNA-binding proteins [69], protein-protein interaction[70], protein fold recognition [62], and protein functional family prediction [68]. Given the sequence of a protein, its amino acid composition and the properties of every constituent amino acid are computed and then used to generate this vector. The computed amino acid properties include hydrophobicity, normalized Van der Waals volume, polarity, polarizability, charge, surface tension, secondary structure and solvent accessibility [68].

For each of these properties, amino acids are divided into three groups such that those in a particular group are regarded to have approximately the same property. For instance, amino acids can be divided into hydrophobic (CVLIMFW), neutral (GASTPHY), and polar (RKEDQN) groups. Three descriptors, composition (C), transition (T), and distribution (D), are introduced to describe global composition of each of these properties. C is the number of amino acids of a particular property (such as hydrophobicity) divided by the total number of amino acids in a protein sequence. T characterizes the percent frequency with which amino acids of a particular property is followed by amino acids of a different property. 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.

A hypothetical protein sequence AEAAAEAEEAAAAAEAEEEAAEEAEEEAAE, as shown in Figure 1, has 16 alanines (n1 = 16) and 14 glutamic acids (n2 = 14). The composition for these two amino acids are n1 × 100.00/(n1 + n2) = 53.33 and n2 × 100.00/(n1 + n2) = 46.67 respectively. There are 15 transitions from A to E or from E to A in this sequence and the percent frequency of these transitions is (15/29) × 100.00 = 51.72. The first, 25%, 50%, 75% and 100% of As are located within the first 1, 5, 12, 20, and 29 residues respectively. The D descriptor for As is thus 1/30 × 100.00 = 3.33, 5/30 × 100.00 = 16.67, 12/30 × 100.00 = 40.0, 20/30 × 100.00 = 66.67, 29/30 × 100.00 = 96.67. Likewise, the D descriptor for Es is 6.67, 26.67, 60.0, 76.67, 100.0. Overall, the amino acid composition descriptors for this sequence are C = (53.33, 46.67), T = (51.72), and D = (3.33, 16.67, 40.0, 66.67, 96.67, 6.67, 26.67, 60.0, 76.67, 100.0) respectively. Descriptors for other properties can be computed by a similar procedure.

Figure 1
figure 1

The sequence of a hypothetic protein for illustration of derivation of the feature vector of a protein. Sequence index indicates the position of an amino acid in the sequence. The index for each type of amino acids in the sequence (A or E) indicates the position of the first, second, third, ... of that type of amino acid (The position of the first, second, third, ..., A is at 1, 3, 4, ...). A/E transition indicates the position of AE or EA pairs in the sequence.

Overall, there are 21 elements representing these three descriptors: 3 for C, 3 for T and 15 for D. The feature vector of a protein is constructed by combining the 21 elements of all of these properties and the 20 elements of amino acid composition in sequential order.

There is some level of overlap in the descriptors for hydrophobicity, polarity, and surface tension. Thus the dimensionality of the feature vectors may be reduced by principle component analysis (PCA). Our own study suggests that the use of PCA reduced feature vectors only moderately improves the accuracy for some of the families. It is thus unclear to which extent this overlap affects the accuracy of SVM classification. It is noted that reasonably accurate results have been obtained using these overlapping descriptors in various protein classification studies [62, 68, 7072].

Support Vector Machines method

The algorithms of SVM and its applications to proteins are extensively described in the literature [68, 69, 73]. Thus only a brief description is given here. A linear SVM constructs a hyperplane that separates two different classes of feature vectors with a maximum margin. One class represents metal-binding proteins and the other non-metal-binding proteins. This hyperplane is constructed by finding a vector w and a parameter b that minimizes ||w||2 which satisfies the following conditions: w · x i + b ≥ 1, for y i = +1(positive class) and w · x i + b ≤ -1, for y i = -1 (negative class). Here x i is a feature vector, y i is the group index, w is a vector normal to the hyperplane, |b|/||w|| is the perpendicular distance from the hyperplane to the origin, and ||w||2 is the Euclidean norm of w.

A nonlinear SVM projects feature vectors into a high dimensional feature space by using a kernel function such as a Gaussian kernel function

. The linear SVM procedure is then applied to the feature vectors in this feature space. After the determination of w and b, a given vector x can be classified by using sign [(w · x) + b], a positive or negative value indicates that the vector x belongs to the positive or negative class respectively.

The performance of SVM has been measured by the positive, negative and overall prediction accuracies P p = TP/(TP + FN), P n = TN/(TN + FP) and P = (TP + TN)/N, which correspond to the accuracies for proteins of a metal-binding class, non-members of the class, and all members and non-members of the class respectively. Here TP, TN, FP, and FN are the number of true positive (true member), true negative (true non-member), false positive (false member), and false negative (false non-member) respectively, and N is the total number of proteins studied.


  1. 1.

    Wintz H, Fox T, Wu YY, Feng V, Chen W, Chang HS, Zhu T, Vulpe C: Expression profiles of Arabidopsis thaliana in mineral deficiencies reveal novel transporters involved in metal homeostasis. J Biol Chem 2003, 278(48):47644–47653. 10.1074/jbc.M309338200

    CAS  Article  PubMed  Google Scholar 

  2. 2.

    Cox EH, McLendon GL: Zinc-dependent protein folding. Curr Opin Chem Biol 2000, 4(2):162–165. 10.1016/S1367-5931(99)00070-8

    CAS  Article  PubMed  Google Scholar 

  3. 3.

    Michel SL, Berg JM: Building a metal binding domain, one half at a time. Chem Biol 2002, 9(6):667–668. 10.1016/S1074-5521(02)00160-6

    CAS  Article  PubMed  Google Scholar 

  4. 4.

    de la Calle Guntinas MB, Bordin G, Rodriguez AR: Identification, characterization and determination of metal-binding proteins by liquid chromatography. A review. Anal Bioanal Chem 2002, 374(3):369–378. 10.1007/s00216-002-1508-3

    CAS  Article  PubMed  Google Scholar 

  5. 5.

    Yang W, Lee HW, Hellinga H, Yang JJ: Structural analysis, identification, and design of calcium-binding sites in proteins. Proteins 2002, 47(3):344–356. 10.1002/prot.10093

    CAS  Article  PubMed  Google Scholar 

  6. 6.

    Jensen MR, Petersen G, Lauritzen C, Pedersen J, Led JJ: Metal binding sites in proteins: identification and characterization by paramagnetic NMR relaxation. Biochemistry 2005, 44(33):11014–11023. 10.1021/bi0508136

    CAS  Article  PubMed  Google Scholar 

  7. 7.

    Wu H, Yang Y, Jiang SJ, Chen LL, Gao HX, Fu QS, Li F, Ma BG, Zhang HY: DCCP and DICP: construction and analyses of databases for copper- and iron-chelating proteins. Genomics Proteomics Bioinformatics 2005, 3(1):52–57.

    CAS  PubMed  Google Scholar 

  8. 8.

    Hantke K: Iron and metal regulation in bacteria. Curr Opin Microbiol 2001, 4(2):172–177. 10.1016/S1369-5274(00)00184-3

    CAS  Article  PubMed  Google Scholar 

  9. 9.

    Bouton CM, Pevsner J: Effects of lead on gene expression. Neurotoxicology 2000, 21(6):1045–1055.

    CAS  PubMed  Google Scholar 

  10. 10.

    Feng M, Patel D, Dervan JJ, Ceska T, Suck D, Haq I, Sayers JR: Roles of divalent metal ions in flap endonuclease-substrate interactions. Nat Struct Mol Biol 2004, 11(5):450–456. 10.1038/nsmb754

    CAS  Article  PubMed  Google Scholar 

  11. 11.

    Carafoli E: Calcium signaling: a tale for all seasons. Proc Natl Acad Sci U S A 2002, 99(3):1115–1122. 10.1073/pnas.032427999

    PubMed Central  CAS  Article  PubMed  Google Scholar 

  12. 12.

    Harris ED: Cellular copper transport and metabolism. Annu Rev Nutr 2000, 20: 291–310. 10.1146/annurev.nutr.20.1.291

    CAS  Article  PubMed  Google Scholar 

  13. 13.

    O'Halloran TV, Culotta VC: Metallochaperones, an intracellular shuttle service for metal ions. J Biol Chem 2000, 275(33):25057–25060. 10.1074/jbc.R000006200

    Article  PubMed  Google Scholar 

  14. 14.

    Vallee BL, Auld DS: Active-site zinc ligands and activated H2O of zinc enzymes. Proc Natl Acad Sci U S A 1990, 87(1):220–224. 10.1073/pnas.87.1.220

    PubMed Central  CAS  Article  PubMed  Google Scholar 

  15. 15.

    Cobbett C, Goldsbrough P: Phytochelatins and metallothioneins: roles in heavy metal detoxification and homeostasis. Annu Rev Plant Biol 2002, 53: 159–182. 10.1146/annurev.arplant.53.100301.135154

    CAS  Article  PubMed  Google Scholar 

  16. 16.

    Papoyan A, Kochian LV: Identification of Thlaspi caerulescens genes that may be involved in heavy metal hyperaccumulation and tolerance. Characterization of a novel heavy metal transporting ATPase. Plant Physiol 2004, 136(3):3814–3823. 10.1104/pp.104.044503

    PubMed Central  CAS  Article  PubMed  Google Scholar 

  17. 17.

    Zhou T, Hamer DH, Hendrickson WA, Sattentau QJ, Kwong PD: Interfacial metal and antibody recognition. Proc Natl Acad Sci U S A 2005, 102(41):14575–14580. 10.1073/pnas.0507267102

    PubMed Central  CAS  Article  PubMed  Google Scholar 

  18. 18.

    Lieu PT, Heiskala M, Peterson PA, Yang Y: The roles of iron in health and disease. Mol Aspects Med 2001, 22(1–2):1–87. 10.1016/S0098-2997(00)00006-6

    CAS  Article  PubMed  Google Scholar 

  19. 19.

    Barondeau DP, Getzoff ED: Structural insights into protein-metal ion partnerships. Curr Opin Struct Biol 2004, 14(6):765–774. 10.1016/

    CAS  Article  PubMed  Google Scholar 

  20. 20.

    Sandier A, Amiel C, Sebille B, Rouchaud JC, Fedoroff M, Soltes L: Chromatographic method involving inductively coupled plasma atomic emission spectrometric detection for the study of metal-protein complexes. J Chromatogr A 1997, 776(1):93–100. 10.1016/S0021-9673(97)00450-0

    CAS  Article  PubMed  Google Scholar 

  21. 21.

    Reed GH, Poyner RR: Mn2+ as a probe of divalent metal ion binding and function in enzymes and other proteins. Met Ions Biol Syst 2000, 37: 183–207.

    CAS  PubMed  Google Scholar 

  22. 22.

    Binet MR, Ma R, McLeod CW, Poole RK: Detection and characterization of zinc- and cadmium-binding proteins in Escherichia coli by gel electrophoresis and laser ablation-inductively coupled plasma-mass spectrometry. Anal Biochem 2003, 318(1):30–38. 10.1016/S0003-2697(03)00190-8

    CAS  Article  PubMed  Google Scholar 

  23. 23.

    Herald VL, Heazlewood JL, Day DA, Millar AH: Proteomic identification of divalent metal cation binding proteins in plant mitochondria. FEBS Lett 2003, 537(1–3):96–100. 10.1016/S0014-5793(03)00101-7

    CAS  Article  PubMed  Google Scholar 

  24. 24.

    Schnepf R, Haehnel W, Wieghardt K, Hildebrandt P: Spectroscopic identification of different types of copper centers generated in synthetic four-helix bundle proteins. J Am Chem Soc 2004, 126(44):14389–14399. 10.1021/ja0484294

    CAS  Article  PubMed  Google Scholar 

  25. 25.

    Gregory DS, Martin AC, Cheetham JC, Rees AR: The prediction and characterization of metal binding sites in proteins. Protein Eng 1993, 6(1):29–35. 10.1093/protein/6.1.29

    CAS  Article  PubMed  Google Scholar 

  26. 26.

    Andreini C, Banci L, Bertini I, Rosato A: Counting the zinc-proteins encoded in the human genome. J Proteome Res 2006, 5(1):196–201. 10.1021/pr050361j

    CAS  Article  PubMed  Google Scholar 

  27. 27.

    Sodhi JS, Bryson K, McGuffin LJ, Ward JJ, Wernisch L, Jones DT: Predicting metal-binding site residues in low-resolution structural models. J Mol Biol 2004, 342(1):307–320. 10.1016/j.jmb.2004.07.019

    CAS  Article  PubMed  Google Scholar 

  28. 28.

    Nayal M, Di Cera E: Predicting Ca(2+)-binding sites in proteins. Proc Natl Acad Sci U S A 1994, 91(2):817–821. 10.1073/pnas.91.2.817

    PubMed Central  CAS  Article  PubMed  Google Scholar 

  29. 29.

    Schymkowitz JW, Rousseau F, Martins IC, Ferkinghoff-Borg J, Stricher F, Serrano L: Prediction of water and metal binding sites and their affinities by using the Fold-X force field. Proc Natl Acad Sci U S A 2005, 102(29):10147–10152. 10.1073/pnas.0501980102

    PubMed Central  CAS  Article  PubMed  Google Scholar 

  30. 30.

    Khalili M, Saunders JA, Liwo A, Oldziej S, Scheraga HA: A united residue force-field for calcium-protein interactions. Protein Sci 2004, 13(10):2725–2735. 10.1110/ps.04878904

    PubMed Central  CAS  Article  PubMed  Google Scholar 

  31. 31.

    Ettema TJ, Huynen MA, de Vos WM, van der Oost J: TRASH: a novel metal-binding domain predicted to be involved in heavy-metal sensing, trafficking and resistance. Trends Biochem Sci 2003, 28(4):170–173. 10.1016/S0968-0004(03)00037-9

    CAS  Article  PubMed  Google Scholar 

  32. 32.

    Rigden DJ, Galperin MY: The DxDxDG motif for calcium binding: multiple structural contexts and implications for evolution. J Mol Biol 2004, 343(4):971–984. 10.1016/j.jmb.2004.08.077

    CAS  Article  PubMed  Google Scholar 

  33. 33.

    Lin CT, Lin KL, Yang CH, Chung IF, Huang CD, Yang YS: Protein metal binding residue prediction based on neural networks. Int J Neural Syst 2005, 15(1–2):71–84. 10.1142/S0129065705000116

    Article  PubMed  Google Scholar 

  34. 34.

    Birch PJ, Dekker LV, James IF, Southan A, Cronk D: Strategies to identify ion channel modulators: current and novel approaches to target neuropathic pain. Drug Discov Today 2004, 9(9):410–418. 10.1016/S1359-6446(04)03043-0

    CAS  Article  PubMed  Google Scholar 

  35. 35.

    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.

    CAS  Article  PubMed  Google Scholar 

  36. 36.

    Lin HH, Han LY, Zhang HL, Zheng CJ, Xie B, Chen YZ: Prediction of the functional class of lipid-binding proteins from sequence derived properties irrespective of sequence similarity. J Lipid Res 2006, 47(4):824–831. 10.1194/jlr.M500530-JLR200

    CAS  Article  PubMed  Google Scholar 

  37. 37.

    Lin HH, Han LY, Cai CZ, Ji ZL, Chen YZ: Prediction of transporter family from protein sequence by support vector machine approach. Proteins 2006, 62(1):218–231. 10.1002/prot.20605

    CAS  Article  PubMed  Google Scholar 

  38. 38.

    Bateman A, Birney E, Cerruti L, Durbin R, Etwiller L, Eddy SR, Griffiths-Jones S, Howe KL, Marshall M, Sonnhammer EL: The Pfam protein families database. Nucleic Acids Res 2002, 30(1):276–280. 10.1093/nar/30.1.276

    PubMed Central  CAS  Article  PubMed  Google Scholar 

  39. 39.

    Frausto da Silva JJR, Williams RJP: The biological chemistry of the elements: The inorganic chemistry of life. New York: Oxford University Press; 1991.

    Google Scholar 

  40. 40.

    Fierro-Monti I, Mathews MB: Proteins binding to duplexed RNA: one motif, multiple functions. Trends Biochem Sci 2000, 25(5):241–246. 10.1016/S0968-0004(00)01580-2

    CAS  Article  PubMed  Google Scholar 

  41. 41.

    Perez-Canadillas JM, Varani G: Recent advances in RNA-protein recognition. Curr Opin Struct Biol 2001, 11(1):53–58. 10.1016/S0959-440X(00)00164-0

    CAS  Article  PubMed  Google Scholar 

  42. 42.

    Bairoch A, Apweiler R: The SWISS-PROT protein sequence database and its supplement TrEMBL in 2000. Nucleic Acids Res 2000, 28(1):45–48. 10.1093/nar/28.1.45

    PubMed Central  CAS  Article  PubMed  Google Scholar 

  43. 43.

    Veropoulos K, Campbell C, Cristianini N: Controlling the sensitivity of Support Vector machines. In Proceedings of the International Joint Conference on Artificial Intelligence (UCAI99). Edited by: Dean T. Sweden: Morgan Kaufmann; 1999:55–60.

    Google Scholar 

  44. 44.

    Brown MP, Grundy WN, Lin D, Cristianini N, Sugnet CW, Furey TS, Ares M Jr, Haussler D: Knowledge-based analysis of microarray gene expression data by using support vector machines. Proc Natl Acad Sci U S A 2000, 97(1):262–267. 10.1073/pnas.97.1.262

    PubMed Central  CAS  Article  PubMed  Google Scholar 

  45. 45.

    Han LY, Cai CZ, Ji ZL, Cao ZW, Cui J, Chen YZ: Predicting functional family of novel enzymes irrespective of sequence similarity: a statistical learning approach. Nucleic Acids Res 2004, 32(21):6437–6444. 10.1093/nar/gkh984

    PubMed Central  CAS  Article  PubMed  Google Scholar 

  46. 46.

    Altschul SF, Madden TL, Schaffer AA, Zhang J, Zhang Z, Miller W, Lipman DJ: Gapped BLAST and PSI-BLAST: a new generation of protein database search programs. Nucleic Acids Res 1997, 25(17):3389–3402. 10.1093/nar/25.17.3389

    PubMed Central  CAS  Article  PubMed  Google Scholar 

  47. 47.

    Skirgaila R, Grazulis S, Bozic D, Huber R, Siksnys V: Structure-based redesign of the catalytic/metal binding site of Cfr10I restriction endonuclease reveals importance of spatial rather than sequence conservation of active centre residues. J Mol Biol 1998, 279(2):473–481. 10.1006/jmbi.1998.1803

    CAS  Article  PubMed  Google Scholar 

  48. 48.

    Thickman KR, Davis A, Berg JM: Site selection in tandem arrays of metal-binding domains. Inorg Chem 2004, 43(25):7897–7901. 10.1021/ic048850t

    CAS  Article  PubMed  Google Scholar 

  49. 49.

    Ledwidge R, Patel B, Dong A, Fiedler D, Falkowski M, Zelikova J, Summers AO, Pai EF, Miller SM: NmerA, the metal binding domain of mercuric ion reductase, removes Hg2+ from proteins, delivers it to the catalytic core, and protects cells under glutathione-depleted conditions. Biochemistry 2005, 44(34):11402–11416. 10.1021/bi050519d

    CAS  Article  PubMed  Google Scholar 

  50. 50.

    Evans RM, Hollenberg SM: Zinc fingers: gilt by association. Cell 1988, 52(1):1–3. 10.1016/0092-8674(88)90522-3

    CAS  Article  PubMed  Google Scholar 

  51. 51.

    Grabarek Z: Structural basis for diversity of the EF-hand calcium-binding proteins. J Mol Biol 2006, 359(3):509–525. 10.1016/j.jmb.2006.03.066

    CAS  Article  PubMed  Google Scholar 

  52. 52.

    Sweeney WV, Rabinowitz JC: Proteins containing 4Fe-4S clusters: an overview. Annu Rev Biochem 1980, 49: 139–161. 10.1146/

    CAS  Article  PubMed  Google Scholar 

  53. 53.

    Laity JH, Lee BM, Wright PE: Zinc finger proteins: new insights into structural and functional diversity. Curr Opin Struct Biol 2001, 11(1):39–46. 10.1016/S0959-440X(00)00167-6

    CAS  Article  PubMed  Google Scholar 

  54. 54.

    Barrera FN, Poveda JA, Gonzalez-Ros JM, Neira JL: Binding of the C-terminal sterile alpha motif (SAM) domain of human p73 to lipid membranes. J Biol Chem 2003, 278(47):46878–46885. 10.1074/jbc.M307846200

    CAS  Article  PubMed  Google Scholar 

  55. 55.

    Chang S, ran Ma T, Miranda RD, Balestra ME, Mahley RW, Huang Y: Lipid- and receptor-binding regions of apolipoprotein E4 fragments act in concert to cause mitochondrial dysfunction and neurotoxicity. Proc Natl Acad Sci U S A 2005, 102(51):18694–18699. 10.1073/pnas.0508254102

    PubMed Central  CAS  Article  PubMed  Google Scholar 

  56. 56.

    Chen MH, Ben-Efraim I, Mitrousis G, Walker-Kopp N, Sims PJ, Cingolani G: Phospholipid scramblase 1 contains a nonclassical nuclear localization signal with unique binding site in importin alpha. J Biol Chem 2005, 280(11):10599–10606. 10.1074/jbc.M413194200

    CAS  Article  PubMed  Google Scholar 

  57. 57.

    Vishwanathan SVN, Smola AJ: Fast Kernels for String and Tree Matching. In Proceedings of Neural Information Processing Systems 2002 2002.

    Google Scholar 

  58. 58.

    Ratsch G, Sonnenburg S, Scholkopf B: RASE: recognition of alternatively spliced exons in C.elegans. Bioinformatics 2005, 21(Suppl 1):i369-i377. 10.1093/bioinformatics/bti1053

    Article  PubMed  Google Scholar 

  59. 59.

    Kuang R, Ie E, Wang K, Wang K, Siddiqi M, Freund Y, Leslie C: Profile-based string kernels for remote homology detection and motif extraction. J Bioinform Comput Biol 2005, 3(3):527–550. 10.1142/S021972000500120X

    CAS  Article  PubMed  Google Scholar 

  60. 60.

    Leslie C, Kuang R, Eskin E: Inexact matching string kernels for protein classification. In Kernel Methods in Computational Biology. Cambridge: MIT Press; 2003:95–112.

    Google Scholar 

  61. 61.

    Han LY, Cai CZ, Lo SL, Chung MC, Chen YZ: Prediction of RNA-binding proteins from primary sequence by a support vector machine approach. Rna 2004, 10(3):355–368. 10.1261/rna.5890304

    PubMed Central  Article  PubMed  Google Scholar 

  62. 62.

    Ding CH, Dubchak I: Multi-class protein fold recognition using support vector machines and neural networks. Bioinformatics 2001, 17(4):349–358. 10.1093/bioinformatics/17.4.349

    CAS  Article  PubMed  Google Scholar 

  63. 63.

    Hunt JA, Ahmed M, Fierke CA: Metal binding specificity in carbonic anhydrase is influenced by conserved hydrophobic core residues. Biochemistry 1999, 38(28):9054–9062. 10.1021/bi9900166

    CAS  Article  PubMed  Google Scholar 

  64. 64.

    Rapisarda VA, Chehin RN, De Las Rivas J, Rodriguez-Montelongo L, Farias RN, Massa EM: Evidence for Cu(I)-thiolate ligation and prediction of a putative copper-binding site in the Escherichia coli NADH dehydrogenase-2. Arch Biochem Biophys 2002, 405(1):87–94. 10.1016/S0003-9861(02)00277-1

    CAS  Article  PubMed  Google Scholar 

  65. 65.

    Abbott JJ, Pei J, Ford JL, Qi Y, Grishin VN, Pitcher LA, Phillips MA, Grishin NV: Structure prediction and active site analysis of the metal binding determinants in gamma -glutamylcysteine synthetase. J Biol Chem 2001, 276(45):42099–42107. 10.1074/jbc.M104672200

    CAS  Article  PubMed  Google Scholar 

  66. 66.

    Maglio O, Nastri F, Calhoun JR, Lahr S, Wade H, Pavone V, DeGrado WF, Lombardi A: Artificial di-iron proteins: solution characterization of four helix bundles containing two distinct types of inter-helical loops. J Biol Inorg Chem 2005, 10(5):539–549. 10.1007/s00775-005-0002-8

    CAS  Article  PubMed  Google Scholar 

  67. 67.

    SVMProt Server[]

  68. 68.

    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. 10.1093/nar/gkg600

    PubMed Central  CAS  Article  PubMed  Google Scholar 

  69. 69.

    Cai CZ, Han LY, Ji ZL, Chen YZ: Enzyme family classification by support vector machines. Proteins 2004, 55(1):66–76. 10.1002/prot.20045

    CAS  Article  PubMed  Google Scholar 

  70. 70.

    Bock JR, Gough DA: Predicting protein – protein interactions from primary structure. Bioinformatics 2001, 17(5):455–460. 10.1093/bioinformatics/17.5.455

    CAS  Article  PubMed  Google Scholar 

  71. 71.

    Cai YD, Liu XJ, Xu XB, Chou KC: Support Vector Machines for predicting HIV protease cleavage sites in protein. J Comput Chem 2002, 23(2):267–274. 10.1002/jcc.10017

    CAS  Article  PubMed  Google Scholar 

  72. 72.

    Cai YD, Liu XJ, Xu XB, Chou KC: Prediction of protein structural classes by support vector machines. Comput Chem 2002, 26(3):293–296. 10.1016/S0097-8485(01)00113-9

    CAS  Article  PubMed  Google Scholar 

  73. 73.

    Burges CJC: A tutorial on support vector machine for pattern recognition. Data Min Knowl Disc 1998, 2: 121–167. 10.1023/A:1009715923555

    Article  Google Scholar 

Download references


This work was supported in part by grants from Ministry of Science and Technology China (2003CB715900, 2004CB720103), National Natural Science Foundation of China (30500107, 30670953), and Science and technology commission of Shanghai municipality (04DZ19850, 06PJ14972).

This article has been published as part of BMC Bioinformatics Volume 7, Supplement 5, 2006: APBioNet – Fifth International Conference on Bioinformatics (InCoB2006). The full contents of the supplement are available online at

Author information



Corresponding authors

Correspondence to ZW Cao or YZ Chen.

Additional information

Authors' contributions

HH wrote the programs and implemented the study. LY designed the web service. HL, CJ and BX performed the data analysis. ZWC assisted in the study and revised the paper. YZ conceived of the study and wrote the manuscript. All authors read and approved the final manuscript.

Rights and permissions

Open Access This article is published under license to BioMed Central Ltd. This is an Open Access article is distributed under the terms of the Creative Commons Attribution License ( ), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Reprints and Permissions

About this article

Cite this article

Lin, H., Han, L., Zhang, H. et al. Prediction of the functional class of metal-binding proteins from sequence derived physicochemical properties by support vector machine approach. BMC Bioinformatics 7, S13 (2006).

Download citation

  • Published:

  • DOI:


  • Support Vector Machine
  • Prediction Accuracy
  • Linear Support Vector Machine
  • Swissprot Database
  • Support Vector Machine Prediction