Improved residue contact prediction using support vector machines and a large feature set
© Cheng and Baldi; licensee BioMed Central Ltd. 2007
Received: 28 December 2006
Accepted: 02 April 2007
Published: 02 April 2007
Predicting protein residue-residue contacts is an important 2D prediction task. It is useful for ab initio structure prediction and understanding protein folding. In spite of steady progress over the past decade, contact prediction remains still largely unsolved.
Here we develop a new contact map predictor (SVMcon) that uses support vector machines to predict medium- and long-range contacts. SVMcon integrates profiles, secondary structure, relative solvent accessibility, contact potentials, and other useful features. On the same test data set, SVMcon's accuracy is 4% higher than the latest version of the CMAPpro contact map predictor. SVMcon recently participated in the seventh edition of the Critical Assessment of Techniques for Protein Structure Prediction (CASP7) experiment and was evaluated along with seven other contact map predictors. SVMcon was ranked as one of the top predictors, yielding the second best coverage and accuracy for contacts with sequence separation >= 12 on 13 de novo domains.
We describe SVMcon, a new contact map predictor that uses SVMs and a large set of informative features. SVMcon yields good performance on medium- to long-range contact predictions and can be modularly incorporated into a structure prediction pipeline.
Predicting protein inter-residue contacts is an important 2D structure prediction problem . Contact prediction can help improve analogous fold recognition [2, 3] and ab initio 3D structure prediction . Several algorithms for reconstructing 3D structure from contacts have been developed in both the structure prediction and determination (NMR) literature [5–8]. Contact map prediction is also useful for inferring protein folding rates and pathways [9, 10].
Due to its importance, contact prediction has received considerable attention over the last decade. For instance, contact prediction methods have been evaluated in the fifth, sixth, and seventh editions of the Critical Assessment of Techniques for Protein Structure Prediction (CASP) experiment [11–15]. A number of different methods for predicting contacts have been developed. These methods can be classified roughly into two non-exclusive categories: (1) statistical correlated mutations approaches [16–22]; and (2) machine learning approaches [23–34]. The former uses correlated mutations of residues to predict contacts. The latter uses machine learning methods such as neural networks, self organizing map, hidden Markov models, and support vector machines to predict 2D contacts from the primary sequence, as well as other 1D features such as relative solvent accessibility and secondary structure.
In spite of steady progress, contact map prediction remains however a largely unsolved challenge. Here we describe a method that uses support vector machines together with a large set of informative features to improve contact map prediction. On the same data set, SVMcon outperforms the latest version of the CMAPpro contact map predictor [28, 35] and is ranked as one of the top predictors in the blind and independent CASP7 experiment.
Results and Discussion
We first compare SVMcon with the latest version of CMAPpro on the same benchmark dataset. Then we describe the performance of SVMcon along with several other predictors during the CASP7 experiment.
Comparison with CMAPpro on the same Benchmark
SVMcon is trained to predict medium- to long-range contacts (sequence separation >= 6) as in , which are not captured by local secondary structure. We train SVMcon on the same dataset used to train CMAPpro [28, 35] and test both programs on the same test dataset. The training dataset contains 485 proteins and the test dataset contains 48 proteins. The sequence identity between the training and testing datasets is below 25%, a common threshold for ab initio prediction .
We use sensitivity and specificity to evaluate the performance of SVMcon and CMAPpro. Sensitivity is the percentage of native contacts that are predicted to be contacts. Specificity is the percentage of predicted contacts that are present in the native structure. The contact threshold is set at 8 Å between Ca atoms. The sensitivity and specificity of a predictor depend also on the threshold used to separate 'contact' from 'non-contact' predictions. To compare SVMcon and CMAPpro fairly, we choose to evaluate them at their break-even point, where sensitivity is equal to specificity as in . At the break-even point, the sensitivity and specificity of SVMcon is 27.1%, 4% higher than CMAPpro. Thus on the same benchmark dataset, SVMcon yields a sizable improvement.
We also compare the accuracy of SVMcon with the random uniform baseline algorithm consisting of random independent coin flips to decide whether each residue pair is in contact or not. Since the medium-and long-range contacts account for 2.8% of the total number of residue pairs with linear separation >= 6, the probability for the coin to produce a contact is set to 2.8%. As a result, the sensitivity and specificity of the random baseline algorithm is 2.8% at the break-even point. Thus SVMcon yields a nine-fold improvement over the random baseline.
Since the contact prediction accuracy varies significantly with individual proteins and their structure classes , we calculate the contact prediction specificity (or called accuracy) and sensitivity (or called coverage) for each test protein (Table 1). For each one, we select up to L (protein length) predicted contacts ranked by SVM scores because the total number of true contacts is approximately linear to the protein length . The results show that in many cases (e.g. 1QJPA, 1DZOA, 1MAIA, 1LSRA, 1F4PA, 1MSCA, 1IG5A, 1ELRA, 1J75A), the prediction accuracy and coverage are > 30%.
Detailed Contact Prediction Results on 48 Test Proteins for Sequence Separation >= 6, 12, and 24 respectively.
Separation >= 6
Separation >= 12
Separation >= 12
Contact Prediction Results of Proteins in the Six SCOP Structure Classes.
Separation >= 6
Separation >= 12
Separation >= 24
Comparison with seven other Predictors during CASP7
SVMcon participated in the CASP7 experiment in 2006 and was evaluated along with seven other contact map predictors. The CASP7 evaluation procedure focuses on inter-residue contact predictions with linear sequence separation >= 12 and >= 24 respectively . Up to L/5 of the top predicted contacts were assessed, where L is the length of the target protein. These evaluation metrics are also similar to those used in the past editions of the Critical Assessment of Fully Automated Structure Prediction Methods [43–45] and in the EVA contact evaluation server . We use the similar procedure to compute accuracy (specificity) and coverage (sensitivity) for all server predictors.
We compare SVMcon with the other contact map predictors on the 13 out of 15 CASP7 de novo domains whose structures have been released. The contact map predictors participating in CASP7 include SVMcon, BETApro , SAM-T06 , PROFcon , GajdaPairings, Distill [34, 48], Possum , and GPCPRED . Their contact predictions were downloaded from the CASP7 website.
CASP7 Results of Inter-Residue Contact Predictions of Eight Predictors.
Separation >= 12
Separation >= 24
Another caveat is that the evaluation dataset and scheme we used may be slightly different from the official CASP7 evaluation. Thus, here we only try to evaluate the current state of the art of contact predictors instead of ranking them. For the offical contact evaluation scheme and results, readers are advised to check the assessment paper of the CASP7 contact predictions published in the upcoming supplement issue of the journal Proteins.
Overall, these results on the CASP7 dataset show that the accuracy and coverage of protein contact prediction are still low. However, these results are an important step towards reaching the milestone of an accuracy level of about 30%, required for deriving moderately accurate (low resolution) 3D protein structures from scratch [5, 8, 39–42] (Also, Dr. Yang Zhang, personal communication at the CASP7 conference). In particular, these predictors tend to predict different correct contacts. Thus, a consensus combination of contact map predictors may yield more accurate contact map predictions, which in turn could significantly improve 3D structure reconstruction. Since the stakes of contact map prediction are high, a community-wide effort for improving contact map prediction should be worthwhile (Dr. Burkhard Rost's lecture slides at Columbia University).
It is also worth pointing out that the CASP7 de novo dataset is too small to reliably estimate the accuracy of the predictors. So one should not over-interpret these results. Indeed, when we use a larger CASP de novo dataset of 24 domains classified by Dr. Dylan Chivian from Dr. David Baker's group to evaluate the predictors (results not shown), the accuracy of SVMcon and BETApro are very close for both sequence separations >= 12 and 24, and both remain among the top predictors.
We have described a new contact map predictor (SVMcon) that uses support vector machines to integrate a large number of useful information including profiles, secondary structure, solvent accessibility, contact potentials, residue types, segment window information [24, 32], and protein-level information . The method yields a 4% improvement over the state-of-the art contact map predictor CMAPpro. In the blind CASP7 experiment, SVMcon is ranked as one of the top contact predictors. The method represents an effort toward a good 2D structure prediction. It can be used to improve ab initio structure prediction  and analogous fold recognition [2, 3]. The web server, software, and source code are available at the SVMcon website .
In the comparison with CMAPpro [28, 35], we use the same training and testing datasets. The datasets are redundancy reduced. The pairwise sequence identity of any two sequences is less than 25%. The training and testing datasets contain 485 sequences and 48 sequences respectively.
We use PSI-BLAST to search each sequence against the NCBI non-redundant database and generate a multiple sequence alignment. The number of PSI-BLAST iterations is set to 3. The e-value for selecting a sequence is set to 0.001. The e-value for including a sequence into the construction of a profile is set to 10-10. Multiple sequence alignments are converted into profiles, where each position is associated with a vector denoting the probability of each residue type.
We extract only medium- and long- range residue pairs with sequence separation >= 6 as in , which are not captured by local secondary structures. We use a 8 Å distance threshold between Ca atoms to classify residue pairs into two categories: contact (positive, < 8 Å) or non-contact (negative, >= 8 Å). Since the majority of residue pairs are negative examples, to balance the number of positive and negative examples in the training set we randomly sample and retain only 5% of the negative examples while keeping all positive examples. In total, there are 220,994 negative examples and 94,110 positive examples in the training data set. We keep all negative and positive examples in the test data set. The test data set has 10,498 positive examples and 367,299 negative examples.
We extract five categories of features for each pair of residues at positions i and j to evaluate their contact likelihood. In addition to the new features (e.g. pairwise information features), the choice of most features combines ideas from our previous contact map predictors, disulfide bond predictors [33, 50], and beta sheet topology predictors , and from the PROFcon , the best predictor in CASP6.
Local window feature
We extract local features using a 9-residue window centered at each residue in each residue pair. For each position in the window, we use 21 inputs for the probabilities of the 20 amino acids plus gap, computed from multiple sequence alignments, 3 binary inputs for secondary structure (helix: 100, sheet: 010, coil: 001), 2 binary inputs for relative solvent accessibility (exposed: 10, buried: 01) at 25% threshold, one input for the entropy (- ∑ k p k logp k ) as a measure of local conservation. Here p k is the probability of occurrence of the k-th residue (or gap) at the position under consideration. Secondary structure and relative solvent accessibility are predicted using the SSpro and ACCpro programs in the SCRATCH suite [27, 35, 51]. Thus the two local windows produce 2 × 9 × 27 features.
Pairwise information features
For each pair of positions (i, j) in a multiple sequence alignment, we compute the following features. One input corresponds to the mutual information of the profiles of the two positions (∑ kl p kl log (p kl /(p k p l )), where p kl is the empirical probability of residues (or gap) k and l appearing at the two positions i and j simultaneously. Two other pairwise inputs are computed using the cosine () and correlation () measures on the profiles at positions i and j. Thus some information about correlated mutations is used in the inputs. We also use three inputs to represent Levitt's contact potential , Jernigan's pairwise potential , and Braun's pairwise potential  for the residue pairs in the target sequence.
Residue type features
We classify residues into four categories: non-polar (G, A, V, L, I, P, M, F, W), polar (S, T, N, Q, C, Y), acidic (D, E), basic (K, R, H). These four residue types induce 10 different combinations: non-polar/non-nopolar, non-polar/polar, non-polar/acidic, non-polar/basic, polar/polar, polar/acidic, polar/basic, acidic/acidic, acidic/basic, and basic/basic. We use 10 binary inputs to encode the type of a residue pair.
Central segment window feature
Central segment window corresponding to a window centered at position ⌊(i + j)/2)⌋ has been shown to be useful for predicting whether the residues at position i and j are in contact or not [24, 32]. We use a central segment window of size 5. For each position in the window, we use the same 27 features as the local window features. So the total number of features for the central window is 5 × 27. We also compute the amino acid composition (21 features), secondary structure composition (3 features), relative solvent accessibility composition (2 features) in the central segment window. The sequence separation (|i - j + 1|) for residue pair (i, j) are classified into one of 16 length bins (< 6, 6, 7, 8, 9, 10, 11, 12, 13, 14, < 19, < 24, <= 29, <= 39, <= 49, >= 50) using a binary vector of length 16, as in .
Protein information features
We also compute the global amino acid composition (21 features), secondary structure composition (3 features), and relative solvent accessibility composition (2 features) of the target sequence. In addition, we classify sequence lengths into four bins (<= 50, <= 100, <= 150, and > 150) using a binary vector of length 4, as in .
Feature selection is useful to improve the performance of machine learning methods, particularly when there is a large number of features as in this study. However, since there are more than 310,000 training data points, it takes about 12 days to conduct a round of training and testing on a Pentium-IV computer. Thus a thorough feature selection is currently not feasible. So we tried only to remove some features (pairwise profile correlation, pairwise mutual information, residue type, and protein information features) once a time to test how they affect the prediction accuracy. We find that removing these features slightly improve the accuracy by about 0.2%. However, it is not clear if the improvement is due to the random variation or due to the removal of the features. But at least, these features are not essential or being compensated by other similar features. Thus, a more thorough feature selection should be conducted to improve the performance when more computing power is available.
where α i are non-negative weights assigned to the training data point x i during training by minimizing a quadratic objective function and b is the bias. Thus the function f (x) can be viewed as a weighted linear combination of similarities between training data points x i and the target data point x. Only data points with strictly positive weight α in the training dataset affect the final solution. The corresponding data points x i are called the support vectors. For contact map prediction, a new data point x is predicted to be positive or negative by taking the sign of f (x).
We use SVM-light [59–61] to implement SVM classification on our data. We experimented with several common kernels including linear kernels, Gaussian radial basis kernels (RBF), polynomial kernels, and sigmoidal kernels. In our experience, on this data the RBF kernel K (x, y) = (or ) works the best. Using the RBF kernel, f (x) is actually a weighted sum of Gaussians centered on the support vectors. Almost any separating boundary or regression function can be obtained with such a kernel , thus it is important to tune the SVM parameters carefully in order to achieve good generalization performance and avoid overfitting.
We only adjust the width parameter γ of the RBF kernel, leaving all other parameters to their default value. γ is the inverse of the variance (σ2) of the RBF and controls the width of the Gaussian functions centered on the support vectors. The bigger is γ, the more peaked are the Gaussians, and the more complex are the resulting decision boundaries . After experimenting with several values of γ, we selected γ = 0.025.
Work supported by NIH grant LM-07443-01 and NSF grants EIA-0321390 and IIS-0513376 to PB. JC is currently supported by a UCF faculty start-up grant.
- Rost B, Liu J, Przybylski D, Nair R, Wrzeszczynski K, Bigelow H, Ofran Y: Prediction of protein structure through evolution. In Handbook of Chemoinformatics – From Data to Knowledge. Edited by: Gasteiger J, Engel T. New York: Wiley; 2003:1789–1811.Google Scholar
- Olmea O, Rost B, Valencia A: Effective use of sequence correlation and conservation in fold recognition. J Mol Biol 1999, 295: 1221–1239.View ArticleGoogle Scholar
- Cheng J, Baldi P: A Machine Learning Information Retrieval Approach to Protein Fold Recognition. Bioinformatics 2006, 22: 1456–1463.View ArticlePubMedGoogle Scholar
- Bonneau R, Ruczinski I, Tsai J, Baker D: Contact order and ab initio protein structure prediction. Protein Sci 2002, 11: 1937–1944.PubMed CentralView ArticlePubMedGoogle Scholar
- Aszodi A, Gradwell M, Taylor W: Global fold determination from a small number of distance restraints. J Mol Biol 1995, 251: 308–326.View ArticlePubMedGoogle Scholar
- Vendruscolo M, Kussell E, Domany E: Recovery of protein structure from contact maps. Folding and Design 1997, 2: 295–306.View ArticlePubMedGoogle Scholar
- Skolnick J, Kolinski A, Ortiz A: MONSSTER: a method for folding globular proteins with a small number of distance restraints. J Mol Biol 1997, 265: 217–241.View ArticlePubMedGoogle Scholar
- Zhang Y, Skolnick J: Automated structure prediction of weakly homologous proteins on a genomic scal. P.N.A.S 2004, 101: 7594–7599.PubMed CentralView ArticlePubMedGoogle Scholar
- Plaxco K, Simons K, Baker D: Contact order, transition state placement and the refolding rates of single domain proteins. Journal of Molecular Biology 1998, 277: 985–994.View ArticlePubMedGoogle Scholar
- Punta M, Rost B: Protein folding rates estimated from contact predictions. J Mol Biol 2005, 507–512.Google Scholar
- Moult J, Hubbard T, Bryant SH, Fidelis K, Pedersen JT: Critical assessment of methods of protein structure prediction (CASP): round II. Proteins Suppl 1997, 1: 2–6.View ArticleGoogle Scholar
- Moult J, Hubbard T, Bryant SH, Fidelis K, Pedersen JT: Critical assessment of methods of protein structure prediction (CASP): round III. Proteins Suppl 1999, (3):22–29.Google Scholar
- Moult J, Fidelis K, Zemla A, Hubbard T: Critical assessment of methods of protein structure prediction (CASP) – round V. Proteins 2003, 53(Suppl 6):334–339.View ArticlePubMedGoogle Scholar
- Moult J, Fidelis K, Tramontano A, Rost B, Hubbard T: Critical assessment of methods of protein structure prediction (CASP) – round VI. Proteins 2005, 61(S7):3–7.View ArticlePubMedGoogle Scholar
- Grana O, Baker D, MacCallum R, Meiler J, Punta M, Rost B, Tress M, Valencia A: CASP6 assessment of contact prediction. Proteins 2005, 61: 214–224.View ArticlePubMedGoogle Scholar
- Goebel U, Sander C, Schneider R, Valencia A: Correlated mutations and residue contacts in proteins. Proteins 1994, 18: 309–317.View ArticleGoogle Scholar
- Olmea O, Valencia A: Improving contact predictions by the combination of correlated mutations and other sources of sequence information. Fold Des 1997, 2: s25-s32.View ArticlePubMedGoogle Scholar
- Shindyalov I, Kolchanov N, Sander C: Can three-dimensional contacts in protein structure be predicted by analysis of correlated mutation? Protein Eng 1994, 7: 349–358.View ArticlePubMedGoogle Scholar
- Hamilton N, Burrage K, Ragan M, Huber T: Protein contact prediction using patterns of correlation. Proteins 2004, 56: 679–684.View ArticlePubMedGoogle Scholar
- Valencia A, Pazos F: Computational methods for the prediction of protein interactons. Curr Opin Struc Biol 2002, 12: 368–373.View ArticleGoogle Scholar
- Halperin I, Wolfson HJ, Nussinov R: Correlated mutations: Advances and limitations. A Study on fusion proteins and on the Cohesin-Dockerin families. Proteins 2006.Google Scholar
- Kundrotas PJ, Alexov EG: Predicting residue contacts using pragmatic correlated mutations method: reducing the false positives. BMC Bioinformatics 2006, 7: 503.PubMed CentralView ArticlePubMedGoogle Scholar
- Fariselli P, Olmea O, Valencia A, Casadio R: Prediction of contact maps with neural networks and correlated mutations. Protein Engineering 2001, 13: 835–843.View ArticleGoogle Scholar
- Lund O, Frimand K, Gorodkin J, Bohr H, Bohr J, Hansen J, Brunak S: Protein distance constraints predicted by neural networks and probability density functions. Prot Eng 1997, 10(11):1241–1248.View ArticleGoogle Scholar
- Fariselli P, Casadio R: Neural network based predictor of residue contacts in proteins. Protein Engineering 1999, 12: 15–21.View ArticlePubMedGoogle Scholar
- Fariselli P, Olmea O, Valencia A, Casadio R: Progress in predicting inter-residue contacts of proteins with neural networks and correlated mutations. Proteins 2001, (Suppl 5):157–162.Google Scholar
- Pollastri G, Baldi P, Fariselli P, Casadio R: Improved prediction of the number of residue contacts in proteins by recurrent neural networks. Bioinformatics 2001, 17: S234-S242. [Proceedings of the ISMB 2001 Conference]. [Proceedings of the ISMB 2001 Conference].View ArticlePubMedGoogle Scholar
- Pollastri G, Baldi P: Prediction of contact maps by GIOHMMs and recurrent neural networks using lateral propagation from all four cardinal corners. Bioinformatics 2002, 18(Suppl 1):S62-S70. [Proceeding of the ISMB 2002 Conference]. [Proceeding of the ISMB 2002 Conference].View ArticlePubMedGoogle Scholar
- MacCallum R: Striped Sheets and Protein Contact Prediction. Bioinformatics 2004, 20(Supplement 1):i224-i231. [Proceedings of the ISMB 2004 Conference]. [Proceedings of the ISMB 2004 Conference].View ArticlePubMedGoogle Scholar
- Shao Y, Bystroff C: Predicting inter-residue contacts using templates and pathways. Proteins 2003, 53(Supplement 6):497–502.View ArticlePubMedGoogle Scholar
- Zhao Y, Karypis G: Prediction of Contact Maps Using Support Vector Machines. Proc of the IEEE Symposium on Bioinformatics and BioEngineering 2003, 26–36.Google Scholar
- Punta M, Rost B: PROFcon: novel prediction of long-range contacts. Bioinformatics 2005, 21: 2960–2968.View ArticlePubMedGoogle Scholar
- Cheng J, Saigo H, Baldi P: Large-Scale Prediction of Disulphide Bridges Using Kernel Methods, Two-Dimensional Recursive Neural Networks, and Weighted Graph Matching. Proteins: Structure, Function, Bioinformatics 2006, 62(3):617–629.View ArticleGoogle Scholar
- Vullo A, Walsh I, Pollastri G: A two-stage approach for improved prediction of residue contact maps. BMC Bioinformatics 2006, 7: 180.PubMed CentralView ArticlePubMedGoogle Scholar
- Cheng J, Randall A, Sweredoski M, Baldi P: SCRATCH: a protein structure and structural feature prediction server. Nucleic Acids Research 2005, (33 web server):w72–76.Google Scholar
- Rost B, Eyrich V: EVA: large-scale analysis of secondary structure prediction. Proteins 2001, 45(S5):192–199.View ArticleGoogle Scholar
- Cheng J, Baldi P: Three-Stage Prediction of Protein Beta-Sheets by Neural Networks, Alignments, and Graph Algorithms. Bioinformatics 2005, 21(suppl 1):i75-i84.View ArticlePubMedGoogle Scholar
- Murzin A, Brenner S, Hubbard T, Chothia C: SCOP: a structural classification of proteins database for the investigation of sequences and structures. Journal of Molecular Biology 1995, 247: 536–540.PubMedGoogle Scholar
- Skolnick J, Kolinski A, Ortiz A: MONSTER: A method for folding globular Proteins with a small number of distance restraints. J Mol Biol 1997, 265: 217–241.View ArticlePubMedGoogle Scholar
- Ortiz A, Kolinski A, Rotkiewicz P, Ilkowski B, Skolnick J: Ab initio folding of proteins using restraints derived from evolutionary information. Proteins Suppl 1999, 3: 177–185.View ArticleGoogle Scholar
- Ortiz A, Kolinski A, Skolnick J: Fold assembly of small proteins using Monte Carlo simulations driven by restraints derived from multiple sequence alignments. J Mol Bio 1998, 227: 419–448.View ArticleGoogle Scholar
- Zhang Y, Kolinski A, Skolnick J: TOUCHSTONE II: a new approach to ab initio protein structure prediction. Biophysical Journal 2003, 85: 1145–1164.PubMed CentralView ArticlePubMedGoogle Scholar
- Fischer D, Barret C, Bryson K, Elofsson A, Godzik A, Jones D, Karplus K, Kelley L, MacCallum R, Pawowski K, Rost B, Rychlewski L, Sternberg M: CAFASP-1: Critical assessment of fully automated structure prediction methods. Proteins 1999, (Suppl 3):209–217.Google Scholar
- Lesk A, Conte LL, Hubbard T: Assessment of novel fold targets in CASP4: predictions of three-dimensional structures, secondary structures, and interresidue contacts. Proteins 2001, 45(S5):98–118.View ArticleGoogle Scholar
- Fischer D, Elofsson A, Rychlewski L, Pazos F, Valencia A, Godzik A, Rost B, Ortiz A, Dunbrack R: CAFASP-2: the second critical assessment of fully automated structure prediction methods. Proteins 2001, 45(S5):171–183.View ArticleGoogle Scholar
- Grana O, Eyrich V, Pazos F, Rost B, Valencia A: EVAcon: a protein contact prediction evaluaton. Nucleic Acid Res 2005, 33: W347-W351.PubMed CentralView ArticlePubMedGoogle Scholar
- Karplus K, Barrett C, Hughey R: Hidden Markov models for detecting remote protein homologies. Bioinformatics 1998, 14(10):846–56.View ArticlePubMedGoogle Scholar
- Bau D, Martin A, Mooney C, Vullo A, Walsh I, Pollastri G: Distill: a suite of web servers for the prediction of one-, two- and three-dimensional structural features of proteins. BMC Bioinformatics 2006, 7: 402.PubMed CentralView ArticlePubMedGoogle Scholar
- Baldi P, Cheng J, Vullo A: Large-scale prediction of disulphide bond connectivity. In Advances in Neural Information Processing Systems (NIPS04 Conference). Volume 17. Edited by: Saul L LB Y Weiss. Cambridge, MA: MIT press; 2005:97–104.Google Scholar
- Pollastri G, Przybylski D, Rost B, Baldi P: Improving the prediction of protein secondary structure in three and eight classes using recurrent neural networks and profiles. Proteins 2002, 47: 228–235.View ArticlePubMedGoogle Scholar
- Huang E, Subbiah S, Tsai J, Levitt M: Using a Hydrophobic Contact Potential to Evaluate Native and Near-Native Folds Generated by Molecular Dynamics Simulations. J Mol Biol 1996, 257: 716–725.View ArticlePubMedGoogle Scholar
- Miyazawa S, Jernigan R: An empirical energy potential with a reference state for protein fold and sequence recognition. Proteins 1999, 36: 357–369.View ArticlePubMedGoogle Scholar
- Zhu H, Braun W: Sequence specificity, statistical potentials, and three-dimensional structure prediction with self-correcting. Protein Sci 1999, 8: 326–342.PubMed CentralView ArticlePubMedGoogle Scholar
- Vapnik V: Statistical Learning Theory. New York, NY: Wiley; 1998.Google Scholar
- Vapnik V: The Nature of Statistical Learning Theory. Berlin, Germany: Springer-Verlag; 1995.View ArticleGoogle Scholar
- Drucker H, Burges C, Kaufman L, Smola A, Vapnik V: Support Vector Regression Machines. In Advances in Neural Information Processing Systems. Volume 9. Edited by: Mozer MC TP MI Jordan. Cambridge, MA: MIT Press; 1997:155–161.Google Scholar
- Schölkopf B, Smola A: Learning with Kernels: Support Vector Machines, Regularization, Optimization, and Beyond. Cambridge, MA: MIT Press; 2002.Google Scholar
- Joachims T: Making large-scale SVM Learning Practical. Advances in Kernel Methods – Support Vector Learning. Edited by: Schölkopf B, Burges C, Smola A. MIT Press; 1999.Google Scholar
- Joachims T: Learning to Classify Text Using Support Vector Machines. Dessertation. Springer. 2002.View ArticleGoogle Scholar
- Vert J, Tsuda K, Scholkopf B: A Primer on Kernel Methods. In Kernel Methods in Computational Biology. Edited by: Scholkopf B JV K Tsuda. Cambridge, MA: MIT Press; 2004:55–72.Google Scholar
- Kraulis P: MOLSCRIPT: A program to produce both detailed and schematic plots of protein structure. Journal of Applied Crystallography 1991, 24: 946–950.View ArticleGoogle Scholar
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.