Skip to main content

Prediction of amyloid fibril-forming segments based on a support vector machine

Abstract

Background

Amyloid fibrillar aggregates of proteins or polypeptides are known to be associated with many human diseases. Recent studies suggest that short protein regions trigger this aggregation. Thus, identifying these short peptides is critical for understanding diseases and finding potential therapeutic targets.

Results

We propose a method, named Pafig (P rediction of a myloid fi bril-forming seg ments) based on support vector machines, to identify the hexpeptides associated with amyloid fibrillar aggregates. The features of Pafig were obtained by a two-round selection from AAindex. Using a 10-fold cross validation test on Hexpepset dataset, Pafig performed well with regards to overall accuracy of 81% and Matthews correlation coefficient of 0.63. Pafig was used to predict the potential fibril-forming hexpeptides in all of the 64,000,000 hexpeptides. As a result, approximately 5.08% of hexpeptides showed a high aggregation propensity. In the predicted fibril-forming hexpeptides, the amino acids – alanine, phenylalanine, isoleucine, leucine and valine occurred at the higher frequencies and the amino acids – aspartic acid, glutamic acid, histidine, lysine, arginine and praline, appeared with lower frequencies.

Conclusion

The performance of Pafig indicates that it is a powerful tool for identifying the hexpeptides associated with fibrillar aggregates and will be useful for large-scale analysis of proteomic data.

Background

Understanding protein aggregation has become increasingly important with the discovery of a correlation between amyloid-like fibrils resulting from protein aggregation and diseases such as Alzheimer's disease, Parkinson's disease, transmissible spongiform encephalopathies, and type II diabetes [14].

Not only the intrinsically disordered proteins, but also natively folded proteins can aggregate into amyloid-like fibrils, as has been found in β2-microglobulin, lysozyme, transthyretin, and the prion protein [57]. It is believed that specific continuous regions within the amyloid fibril-forming proteins can act as facilitators or inhibitors of amyloid fibril formation and determine the proteins' aggregation tendency [6, 812]. Therefore, recognizing the specific regions that determine the aggregation propensity of a protein is of fundamental interest, since this will help in understanding the mechanism of amyloid formation and lead to effective treatments for amyloid illnesses [13].

As reviewed recently [13], there are two types of computational approaches used to investigate the aggregation propensity of peptides or proteins and to identify the segments most prone to form fibrils (hot spots). The first approach uses phenomenological models based on the physicochemical properties of the amino acids (e.g. β-propensity, hydrophobicity, aromatic content, and charge) to predict changes in aggregation rates with mutation as well as absolute aggregation rates and hot spots [5, 6, 1419]. The second approach combines atomistic simulations of a protein segment with the microcrystal structure of short fibril-forming peptides to gain insight into aggregation propensity [1, 2022]. This approach may help to elucidate the structural details of ordered aggregates. In addition to the approaches described above, a sequence pattern obtained by saturation mutagenesis analysis [11] has been proposed to identify amyloidogenic stretches in proteins.

In this study, we proposed the Pafig (P rediction of a myloid fi bril-forming seg ments) to identify fibril-forming segments in proteins based on a support vector machine (SVM) [23, 24]. The predictive model of Pafig was also a phenomenological model, which was based on 41 physicochemical properties selected by a two-round selection from 531 physicochemical properties in the Amino acid index database (AAindex) [25, 26]. Because short regions of a protein were responsible for its amyloidogenic behavior [1, 5, 6, 11, 18], Pafig was trained by hexpeptides, which were decomposed by scanning for segments that could form fibrils with a six-residue sliding window. Pafig is simple, fast, and suitable for large-scale calculations such as proteomics analyses.

Methods

Datasets and physicochemical properties

In this study, we constructed a dataset, the Hexpepset dataset, to train the model and test the robustness of Pafig. The six-residue peptides that could form fibrils were defined as positive samples and those that could not form fibrils as negative samples. The Hexpepset dataset consisted of 2452 hexpeptides (1226 positive samples and 1226 negative samples). The positive samples in the Hexpepset dataset were collected by scanning known fibril-forming fragments [1, 5, 6, 8, 14, 18, 20, 27] (Additional file 1) with a six-residue window. Because a large difference between the positive and negative samples can hamper training of the SVM [28, 29], the negative set also contained 1226 hexpeptides selected from two different sources. The first negative sample part contained 876 samples, which was selected by scanning the fragments that had been proved not to form fibrils by experiments[18, 20]. The second negative sample part consisted of 350 samples. These samples were randomly selected from the hexpeptides obtained by scanning the fragments except the experimentally determined amyloidogenic regions of five proteins (Transthyretin, Genbank Accession No. AAB35640.1; Major prion protein precursor, Genbank Accession No. P04156; Apo-AI, Genbank Accession No. P02647; Alpha-synuclein, Genbank Accession No. P37840.1; Beta-2-Microglobulin, Genbank Accession No. 1LDS) [1, 6, 27]. The Hexpepset dataset could be downloaded from website of Pafig [30].

There are totally 544 physicochemical properties in the amino acid index database version 9.0 (AAindex) [25, 26], which is a collection of published amino acid indices representing different physicochemical and biological properties of amino acids. Each physicochemical property consists of a set of 20 numerical values for amino acids. The property having the value 'NA' in a value set of amino acid index was discarded. Finally, 531 properties were used for the following mining method.

Physicochemical properties selection and Encoding Schemes of Pafig

The input feature of Pafig consisted of 123 elements based on 41 physicochemical properties. The properties selection and encoding schemes of Pafig was illustrated in Figure 1 and calculated as follows:

The 531 physicochemical properties were downloaded from Aaindex [25]. The values of each property were scaled to zero mean and a standard deviation of 1.

According to a physicochemical property p i , a hexpeptide hk could be mapped to feature vector F i with three elements. F i k = [ m i 1 k , m i 2 k , m i 3 k ] MathType@MTEF@5@5@+=feaagaart1ev2aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacPC6xNi=xH8viVGI8Gi=hEeeu0xXdbba9frFj0xb9qqpG0dXdb9aspeI8k8fiI+fsY=rqGqVepae9pg0db9vqaiVgFr0xfr=xfr=xc9adbaqaaeGaciGaaiaabeqaaeqabiWaaaGcbaGaemOray0aa0baaeaacqWGPbqAaeaacqWGRbWAaaGaeyypa0Jaei4waSLaemyBa02aa0baaeaacqWGPbqAcqaIXaqmaeaacqWGRbWAaaGaeiilaWIaemyBa02aa0baaeaacqWGPbqAcqaIYaGmaeaacqWGRbWAaaGaeiilaWIaemyBa02aa0baaeaacqWGPbqAcqaIZaWmaeaacqWGRbWAaaGaeiyxa0faaa@449F@ , m i 1 k = n = 1 n = 6 [ x i k ( n ) ] 2 MathType@MTEF@5@5@+=feaagaart1ev2aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacPC6xNi=xH8viVGI8Gi=hEeeu0xXdbba9frFj0xb9qqpG0dXdb9aspeI8k8fiI+fsY=rqGqVepae9pg0db9vqaiVgFr0xfr=xfr=xc9adbaqaaeGaciGaaiaabeqaaeqabiWaaaGcbaGaemyBa02aa0baaeaacqWGPbqAcqaIXaqmaeaacqWGRbWAaaGaeyypa0ZaaabCaeaacqGGBbWwcqWG4baEdaqhaaqaaiabdMgaPbqaaiabdUgaRbaacqGGOaakcqWGUbGBcqGGPaqkcqGGDbqxdaahaaqabeaacqaIYaGmaaaabaGaemOBa4Maeyypa0JaeGymaedabaGaemOBa4Maeyypa0JaeGOnaydacqGHris5aaaa@45F4@ , m i 2 k = n = 1 n = 6 [ x i k ( n ) ] 2 × δ ( n ) MathType@MTEF@5@5@+=feaagaart1ev2aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacPC6xNi=xH8viVGI8Gi=hEeeu0xXdbba9frFj0xb9qqpG0dXdb9aspeI8k8fiI+fsY=rqGqVepae9pg0db9vqaiVgFr0xfr=xfr=xc9adbaqaaeGaciGaaiaabeqaaeqabiWaaaGcbaGaemyBa02aa0baaeaacqWGPbqAcqaIYaGmaeaacqWGRbWAaaGaeyypa0ZaaabCaeaacqGGBbWwcqWG4baEdaqhaaqaaiabdMgaPbqaaiabdUgaRbaacqGGOaakcqWGUbGBcqGGPaqkcqGGDbqxdaahaaqabeaacqaIYaGmaaaabaGaemOBa4Maeyypa0JaeGymaedabaGaemOBa4Maeyypa0JaeGOnaydacqGHris5aiabgEna0kabes7aKjabcIcaOiabd6gaUjabcMcaPaaa@4CC9@ . If the value of x i k ( n ) MathType@MTEF@5@5@+=feaagaart1ev2aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacPC6xNi=xH8viVGI8Gi=hEeeu0xXdbba9frFj0xb9qqpG0dXdb9aspeI8k8fiI+fsY=rqGqVepae9pg0db9vqaiVgFr0xfr=xfr=xc9adbaqaaeGaciGaaiaabeqaaeqabiWaaaGcbaGaemiEaG3aa0baaSqaaiabdMgaPbqaaiabdUgaRbaakiabcIcaOiabd6gaUjabcMcaPaaa@3356@ was greater than zero, the value of δ(n) is 1; otherwise it was 0. m i 3 k = n = 1 n = 6 [ x i k ( n ) ] 2 × δ ( n ) MathType@MTEF@5@5@+=feaagaart1ev2aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacPC6xNi=xH8viVGI8Gi=hEeeu0xXdbba9frFj0xb9qqpG0dXdb9aspeI8k8fiI+fsY=rqGqVepae9pg0db9vqaiVgFr0xfr=xfr=xc9adbaqaaeGaciGaaiaabeqaaeqabiWaaaGcbaGaemyBa02aa0baaeaacqWGPbqAcqaIZaWmaeaacqWGRbWAaaGaeyypa0ZaaabCaeaacqGGBbWwcqWG4baEdaqhaaqaaiabdMgaPbqaaiabdUgaRbaacqGGOaakcqWGUbGBcqGGPaqkcqGGDbqxdaahaaqabeaacqaIYaGmaaaabaGaemOBa4Maeyypa0JaeGymaedabaGaemOBa4Maeyypa0JaeGOnaydacqGHris5aiabgEna0kabes7aKjabcIcaOiabd6gaUjabcMcaPaaa@4CCB@ . If the value of x i k ( n ) MathType@MTEF@5@5@+=feaagaart1ev2aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacPC6xNi=xH8viVGI8Gi=hEeeu0xXdbba9frFj0xb9qqpG0dXdb9aspeI8k8fiI+fsY=rqGqVepae9pg0db9vqaiVgFr0xfr=xfr=xc9adbaqaaeGaciGaaiaabeqaaeqabiWaaaGcbaGaemiEaG3aa0baaSqaaiabdMgaPbqaaiabdUgaRbaakiabcIcaOiabd6gaUjabcMcaPaaa@3356@ was smaller than zero, the value of δ(n) was 1; otherwise it was 0. Here x i k ( n ) MathType@MTEF@5@5@+=feaagaart1ev2aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacPC6xNi=xH8viVGI8Gi=hEeeu0xXdbba9frFj0xb9qqpG0dXdb9aspeI8k8fiI+fsY=rqGqVepae9pg0db9vqaiVgFr0xfr=xfr=xc9adbaqaaeGaciGaaiaabeqaaeqabiWaaaGcbaGaemiEaG3aa0baaSqaaiabdMgaPbqaaiabdUgaRbaakiabcIcaOiabd6gaUjabcMcaPaaa@3356@ was one of physicochemical property corresponding to the nth amino acid.

The Support vector machine was used to evaluate the classification ability of each property. The property, which obtained the overall accuracy greater or equal than 60%, was selected as the candidate of Pafig.

A standard genetic algorithm [31] was used to select the final physicochemical properties of Pafig with population size of 10, crossover probability of 0.8, mutation probability of 0.01 and predetermined number of 200 generations. The overall accuracy (Q2) was adopted as the fitness.

Figure 1
figure1

Physicochemical property selection and encoding schemes of Pafig.

Support vector machine

Support vector machine (SVM) was used to identify an optimal hyperplane in order to separate two classes of samples [23, 24]. SVM uses kernel functions to map the original data to a feature space of higher dimensions and also locates an optimal separating hyperplane. For SVM implementation, we used LIBSVM [32] with a Radial Basis Function (RBF kernel) and a Polynomial kernel. The parameters were selected with the LIBSVM parameter selection tool (easy.py).

Prediction system assessment

True positives (TP) and true negatives (TN) were identified as the positive and negative samples, respectively. False positives (FP) were negative samples identified as positive. False negatives (FN) were positive samples identified as negative. The prediction performance was tested with sensitivity (TP/(TP+FN)), specificity (TN/(TN+FP)), overall accuracy (Q2), and the Matthews correlation coefficient (MCC). The Q2 and MCC were calculated as follows:

M C C = T P × T N F P × F N ( T N + F N ) × ( T N + F P ) × ( T P + F N ) × ( T P + F P ) Q 2 = T P + T N T P + T N + F P + F N MathType@MTEF@5@5@+=feaagaart1ev2aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacPC6xNi=xI8qiVKYPFjYdHaVhbbf9v8qqaqFr0xc9vqFj0dXdbba91qpepeI8k8fiI+fsY=rqGqVepae9pg0db9vqaiVgFr0xfr=xfr=xc9adbaqaaeGaciGaaiaabeqaaeqabiWaaaGcbaqbaeaabiqaaaqaaiabd2eanjabdoeadjabdoeadjabg2da9KqbaoaalaaabaGaemivaqLaemiuaaLaey41aqRaemivaqLaemOta4KaeyOeI0IaemOrayKaemiuaaLaey41aqRaemOrayKaemOta4eabaWaaOaaaeaacqGGOaakcqWGubavcqWGobGtcqGHRaWkcqWGgbGrcqWGobGtcqGGPaqkcqGHxdaTcqGGOaakcqWGubavcqWGobGtcqGHRaWkcqWGgbGrcqWGqbaucqGGPaqkcqGHxdaTcqGGOaakcqWGubavcqWGqbaucqGHRaWkcqWGgbGrcqWGobGtcqGGPaqkcqGHxdaTcqGGOaakcqWGubavcqWGqbaucqGHRaWkcqWGgbGrcqWGqbaucqGGPaqkaeqaaaaaaOqaaKqzaeGaemyuaeLaeGOmaiJaeyypa0tcfa4aaSaaaeaacqWGubavcqWGqbaucqGHRaWkcqWGubavcqWGobGtaeaacqWGubavcqWGqbaucqGHRaWkcqWGubavcqWGobGtcqGHRaWkcqWGgbGrcqWGqbaucqGHRaWkcqWGgbGrcqWGobGtaaaaaaaa@77E6@

Results

Properties selection

The properties of Pafig were selected by a two-round selection on Hexpepset dataset from the AAindex [25, 26]. Firstly, a preliminary properties selection on all of the physicochemical properties (531 properties) was performed using the following steps. (1) Every hexpeptide in the Hexpepset dataset was encoded by each physicochemical property. (2) The overall prediction accuracy (Q2) corresponding to each property was calculated by LIBSVM with a 10-Cross validation on Hexpepset dataset. (3) All performance of each physicochemical property was ranked according to their overall prediction accuracy (Q2). (4) The overall prediction accuracy (Q2 = 60%) was used as the cut off to select properties. As a result, 155 physicochemical properties were selected from 531 properties in AAindex and used as the candidates for further selection. Secondly, an efficient bi-objective genetic algorithm with population size of 10, crossover probability of 0.8, mutation probability of 0.01 and predetermined number of 200 generations was utilized to mine informative physicochemical properties and combine the selected properties to build a prediction model. The overall prediction accuracy (Q2) was adopted as the fitness. All results of genetic algorithm were obtained by LIBSVM with a 10-Cross validation on Hexpepset dataset and were shown in Figure 2.

Figure 2
figure2

Average overall accuracy among the individuals and the best overall accuracy of different generations. The results were obtained by LIBSVM with 10-cross validation on Hexpepset dataset. Q2: overall accuracy.

As shown in Figure 2, the SVM classifier obtained the best overall accuracy Q2 = 81% and MCC = 0.63 in an individual of the 186th generation, where the number of selected physicochemical properties was 41. The standard deviations of Q2 and MCC among the 10-fold cross-validation of the best individual were 1% and 0.01, respectively. Table 1 shows the selected physicochemical properties by Genetic Algorithm and the overall accuracy obtained by the corresponding property. The best property of all the selected properties was VINM940104, which obtained the best overall accuracy Q2 = 66% and corresponds to the "Normalized flexibility parameters (B-values) for each residue surrounded by two rigid neighbours" [33].

Table 1 The AAindex identities and the overall accuracy (Q2) of the 41 properties selected by Genetic Algorithm. The overall prediction accuracy (Q2) was calculated using the corresponding property by LIBSVM with RBF kernel of 10-Cross validation on Hexpepset dataset.

Prediction performance of Pafig

We investigated the performances using 41, 155, or 531 physicochemical properties. The 10-fold cross validation tests on the Hexpepset dataset were carried out by LIBSVM with RBF kernel. As shown in Table 2, the classifier employing the 41 properties obtained the best performance, with the overall prediction accuracy (Q2) of 81% and the Matthews correlation coefficient (MCC) of 0.63, which was better than using 155 and 531 physicochemical properties. This improvement is mainly due to the reduction in noise and outliers present in the 378 and 531 physicochemical properties, which influence the performance of the Support Vector Machine [34]. In addition, we also evaluated the effect of different SVM kernel and found that the performance with the Radial Basis Function (RBF) kernel was better than the Polynomial kernel, which only obtained an overall accuracy (Q2) of 77% and a Matthews correlation coefficient (MCC) of 0.54 using the selected 41 physicochemical properties with best parameter of d = 3.

Table 2 Prediction performance on the Hexpepset dataset with different properties. All of the results were obtained by LIBSVM with RBF kernel.

The receiver operating characteristics (ROC) score was usually used as the primary measure of the machine learning method performance and provided an overview of the possible cut-off levels in the test performance [35]. The ROC curves for the random classifier and classifiers with property of VINM940104, properties with overall prediction accuracy (Q2) ≥ 60% and properties selected by Genetic Algorithm were shown in Figure 3. The result revealed that the classifier with the 41 predictive properties selected by genetic algorithm was better than the property of VINM940104 and properties with overall prediction accuracy (Q2) ≥ 60%.

Figure 3
figure3

ROC curves of random classifier, classifier with property of VINM940104, properties with overall prediction accuracy (Q2) ≥ 60% and properties selected by Genetic Algorithm and some published classifiers (Alg1 [27] Alg2 [5] and Alg3 [6]).

In addition, we also compared the performance of Pafig with other methods [5, 6, 27] using the same dataset. The Alg1 [27] is based on the intrinsic aggregation propensities to identify the regions of the protein sequence that are most important for promoting amyloid formation. The Alg2 [5] and Alg3 [6] were used the observed packing density and the relative experimental aggregation propensities of the 20 natural amino acids to detect the amyloid fibril-forming segments. As shown in the Figure 3, we could clearly found that performance of Pafig was better than these popular methods. Moreover, these results also indicate that Pafig is a powerful tool for predicting the amyloid fibril-forming segments in the protein sequence.

Reliability index for Pafig predictions

When machine learning approaches are selected to classify the samples, it is important to know the reliability of the prediction result [3639]. In this study, a reliability index (RI) was assigned to a predicted hexpeptide based on the output of LIBSVM. Provided that an output of LIBSVM for a hexpeptide is O, the RI value is thus computed as RI = INTEGER (20 × abs(O-0.5)). The value of RI could provide information about the certainty of the classification and could be used as an indicator of prediction certainty for a particular hexpeptide. Figure 4 showed the expected prediction accuracies and the fraction of hexpeptides with a given RI value. For example, about 42% of the hexpeptides had an RI ≥ 7, and of these, 89% were correctly predicted. The result was obtained by LIBSVM with 10-fold cross validations on the Hexpepset dataset.

Figure 4
figure4

Average prediction accuracy calculated cumulatively with RI above a given value. For example, 42% of all hexpeptides had an RI ≥ 7, and of these hexpeptides, 89% were correctly predicted. The results were based on the Hexpepset dataset with 10-fold cross validations.

Identification of fibril forming peptides in all hexpeptides

Pafig was used to predict the potential fibril-forming hexpeptides in all of the 64,000,000 hexpeptides. The fraction of possible fibril-forming hexpeptides with different RI cutoffs was shown in Figure 5. We found that 5.08% of fibril-forming hexpeptides had an RI ≥ 7, which was the minority of all hexpeptides [11]. As shown in Figure 5, the hydrophobic amino acids (alanine, phenylalanine, isoleucine, leucine and valine) occurred at the higher frequencies in the predicted fibril formation hexpeptides. These results matched the hypothesis [6, 27, 40] that hydrophobic residues usually induced aggregation. However, the amino acids with positive or negative charges, such as aspartic acid, glutamic acid, histidine, lysine and arginine appeared in the predicted fibril formation hexpeptides with lower frequencies. In addition, proline was found with the lowest frequency in the predicted fibril-forming hexpeptides, as expected, as proline was a β-sheet breaker and most fibrils had crossing β-structures [11, 41].

Figure 5
figure5

The frequency of amino acid in the predicted fibril-forming hexpeptides with different RI cutoffs. The left column is the different RI cutoffs. PM(%) in the right column is the percentage of predicted fibril-forming hexpeptides in all hexpeptides. The other columns represent the frequency of different amino acids in the predicted fibril-forming hexpeptides with different RI cutoffs. Red represents higher frequency, green represents lower frequency and black represents equal levels of frequency.

Identification of the frequency of fibril-forming hexpeptides in proteins of the UniProt database

The UniProt Knowledgebase (Release 13.5) analysed here was taken from the European Bioinformatics Institute (EBI), since it was the central access point for extensively curated protein information, including function, classification and cross-references. Pafig was used to predict the possible fibril forming hexpeptides of every protein in the database. As shown in Table 3, archaea, bacteria and plants had a consistently higher frequency of fibril-forming hexpeptides. However, human had a relatively lower aggregation propensity. These results also suggested that evolution was as a factor to oppose protein aggregation by minimizing the amount of strongly aggregating sequence stretches [42]. Moreover, we classified the UniProt Knowledgebase based on the gene ontology [43]. As shown in Table 4, the frequency of fibril forming hexpeptides was very similar among the gene ontology Biological process, Cellular component and Molecular function classes. However, the sub-classes of GO:0009055(electron carrier activity) and GO:0005215(transporter activity) in the molecular function class, which consisted of many membrane proteins, had higher frequency of fibril-forming hexpeptides. This result indicated that the aggregation propensity of membrane proteins was somewhat higher than other proteins.

Table 3 Fibril-forming Propensity of proteins in the Uniprot Knowledgebase classed based on taxonomy. The Predicted fibril-forming hexpeptides were identified by obtaining RI ≥ 7.
Table 4 Fibril-forming Propensity of proteins in the Uniprot Knowledgebase classed based on Gene Ontology. The Predicted fibril-forming hexpeptides were identified by obtaining RI ≥ 7.

Discussion

Amyloid fibrillar aggregates of proteins or polypeptides are potentially lethal and related to many diseases, such as Alzheimer's disease, Parkinson's disease, transmissible spongiform encephalopathies, and type II diabetes [14]. Here, we have described a method for predicting the amyloid fibril-forming propensities of hexpeptides. Moreover, we used this approach to identify the possible fibril- forming peptides in all hexpeptides and the frequency of amyloid fibril-forming hexpeptides in proteins of the UniProt database. All of these results will help us in understanding the cause of the amyloid fibrillar aggregates of peptides.

In this study, we did not start from some properties which had been proved to relate with the amyloid fibrillar aggregates of peptides, such as the Packing Density, Hydrophobicity, Charge β-sheet propensity and so on [5, 6, 14, 27]. However, we firstly evaluated every physicochemical property in the AAindex and found some physicochemical properties related to amyloid fibrillar aggregates of proteins or polypeptides. Secondly, the genetic algorithm was used to select a physicochemical property set. The results indicated that the selected properties could complement one another to yield a powerful and efficient predictor. In addition, a hexpeptide with every physicochemical property was encoded by three elements rather than the average value of the corresponding property. The performance of employing this encode scheme was better than using one element (data not shown). The detailed encode schemes was shown in the section of Methods.

The features of Pafig did not directly contain the structural features of the proteins. Thus, it is possible that some of the structure information was ignored by Pafig, which also existed as the protein aggregation and fibril forming factors [1, 20]. However, the lack of structural information was likely overcome by the inclusion of different physicochemical properties in the Pafig [39]. Moreover, the sample size of the training dataset of Pafig compared with the number of all hexpeptides was very small, which would affect the performance of Pafig. Therefore, the future work is to collect more data by combining biological knowledge and related sources and add some structure feature into Pafig.

Conclusion

In this study, we used 41 physicochemical properties to identify the specific regions (six consecutive residues) associated with amyloid fibrillar aggregates. As this method is computationally efficient and accurate, it can be used to analyze large systems, such as entire proteomics data. Moreover, we have found that the amino acids – alanine, phenylalanine, isoleucine, leucine and valine occurred at the higher frequencies and the amino acids – aspartic acid, glutamic acid, histidine, lysine, arginine and proline appeared with the lower frequencies in the predicted fibril-forming hexpeptides. The amyloid fibrillar aggregation propensity of membrane proteins was somewhat higher than other proteins.

Availability and requirements

Project name: Pafig [30]

Operating systems: Windows

Programming language: Perl

License: GNU General Public License. This license allows the source code to be redistributed and/or modified under the terms of the GNU General Public License as published by the Free Software Foundation. The source code for the application is available at no charge.

Any restrictions to use by non-academics: None

References

  1. 1.

    Zhang Z, Chen H, Lai L: Identification of amyloid fibril-forming segments based on structure and residue-based statistical potential. Bioinformatics 2007, 23(17):2218–2225. 10.1093/bioinformatics/btm325

    CAS  Article  PubMed  Google Scholar 

  2. 2.

    Ross CA, Poirier MA: Protein aggregation and neurodegenerative disease. Nat Med 2004, 10(Suppl):S10–17. 10.1038/nm1066

    Article  PubMed  Google Scholar 

  3. 3.

    Dobson CM: The structural basis of protein folding and its links with human disease. Philos Trans R Soc Lond B Biol Sci 2001, 356(1406):133–145. 10.1098/rstb.2000.0758

    PubMed Central  CAS  Article  PubMed  Google Scholar 

  4. 4.

    Rochet JC, Lansbury PT Jr: Amyloid fibrillogenesis: themes and variations. Curr Opin Struct Biol 2000, 10(1):60–68. 10.1016/S0959-440X(99)00049-4

    CAS  Article  PubMed  Google Scholar 

  5. 5.

    Galzitskaya OV, Garbuzynskiy SO, Lobanov MY: Prediction of amyloidogenic and disordered regions in protein chains. PLoS Comput Biol 2006, 2(12):e177. 10.1371/journal.pcbi.0020177

    PubMed Central  Article  PubMed  Google Scholar 

  6. 6.

    Sanchez de Groot N, Pallares I, Aviles FX, Vendrell J, Ventura S: Prediction of "hot spots" of aggregation in disease-linked polypeptides. BMC Struct Biol 2005, 5: 18. 10.1186/1472-6807-5-18

    Article  PubMed  Google Scholar 

  7. 7.

    Kelly JW: The alternative conformations of amyloidogenic proteins and their multi-step assembly pathways. Curr Opin Struct Biol 1998, 8(1):101–106. 10.1016/S0959-440X(98)80016-X

    CAS  Article  PubMed  Google Scholar 

  8. 8.

    Nelson R, Sawaya MR, Balbirnie M, Madsen AO, Riekel C, Grothe R, Eisenberg D: Structure of the cross-beta spine of amyloid-like fibrils. Nature 2005, 435(7043):773–778. 10.1038/nature03680

    PubMed Central  CAS  Article  PubMed  Google Scholar 

  9. 9.

    Makin OS, Atkins E, Sikorski P, Johansson J, Serpell LC: Molecular basis for amyloid fibril formation and stability. Proc Natl Acad Sci USA 2005, 102(2):315–320. 10.1073/pnas.0406847102

    PubMed Central  CAS  Article  PubMed  Google Scholar 

  10. 10.

    Ventura S, Zurdo J, Narayanan S, Parreno M, Mangues R, Reif B, Chiti F, Giannoni E, Dobson CM, Aviles FX, et al.: Short amino acid stretches can mediate amyloid formation in globular proteins: the Src homology 3 (SH3) case. Proc Natl Acad Sci USA 2004, 101(19):7258–7263. 10.1073/pnas.0308249101

    PubMed Central  CAS  Article  PubMed  Google Scholar 

  11. 11.

    Lopez de la Paz M, Serrano L: Sequence determinants of amyloid fibril formation. Proc Natl Acad Sci USA 2004, 101(1):87–92. 10.1073/pnas.2634884100

    PubMed Central  CAS  Article  PubMed  Google Scholar 

  12. 12.

    Ivanova MI, Sawaya MR, Gingery M, Attinger A, Eisenberg D: An amyloid-forming segment of beta2-microglobulin suggests a molecular model for the fibril. Proc Natl Acad Sci USA 2004, 101(29):10584–10589. 10.1073/pnas.0403756101

    PubMed Central  CAS  Article  PubMed  Google Scholar 

  13. 13.

    Caflisch A: Computational models for the prediction of polypeptide aggregation propensity. Curr Opin Chem Biol 2006, 10(5):437–444. 10.1016/j.cbpa.2006.07.009

    CAS  Article  PubMed  Google Scholar 

  14. 14.

    Tartaglia GG, Cavalli A, Pellarin R, Caflisch A: Prediction of aggregation rate and aggregation-prone segments in polypeptide sequences. Protein Sci 2005, 14(10):2723–2734. 10.1110/ps.051471205

    PubMed Central  CAS  Article  PubMed  Google Scholar 

  15. 15.

    Tartaglia GG, Cavalli A, Pellarin R, Caflisch A: The role of aromaticity, exposed surface, and dipole moment in determining protein aggregation rates. Protein Sci 2004, 13(7):1939–1941. 10.1110/ps.04663504

    PubMed Central  CAS  Article  PubMed  Google Scholar 

  16. 16.

    DuBay KF, Pawar AP, Chiti F, Zurdo J, Dobson CM, Vendruscolo M: Prediction of the absolute aggregation rates of amyloidogenic polypeptide chains. J Mol Biol 2004, 341(5):1317–1326. 10.1016/j.jmb.2004.06.043

    CAS  Article  PubMed  Google Scholar 

  17. 17.

    Chiti F, Stefani M, Taddei N, Ramponi G, Dobson CM: Rationalization of the effects of mutations on peptide and protein aggregation rates. Nature 2003, 424(6950):805–808. 10.1038/nature01891

    CAS  Article  PubMed  Google Scholar 

  18. 18.

    Fernandez-Escamilla AM, Rousseau F, Schymkowitz J, Serrano L: Prediction of sequence-dependent and mutational effects on the aggregation of peptides and proteins. Nat Biotechnol 2004, 22(10):1302–1306. 10.1038/nbt1012

    CAS  Article  PubMed  Google Scholar 

  19. 19.

    Idicula-Thomas S, Balaji PV: Understanding the relationship between the primary structure of proteins and their amyloidogenic propensity: clues from inclusion body formation. Protein Eng Des Sel 2005, 18(4):175–180. 10.1093/protein/gzi022

    CAS  Article  PubMed  Google Scholar 

  20. 20.

    Thompson MJ, Sievers SA, Karanicolas J, Ivanova MI, Baker D, Eisenberg D: The 3D profile method for identifying fibril-forming segments of proteins. Proc Natl Acad Sci USA 2006, 103(11):4074–4078. 10.1073/pnas.0511295103

    PubMed Central  CAS  Article  PubMed  Google Scholar 

  21. 21.

    Yoon S, Welsh WJ: Detecting hidden sequence propensity for amyloid fibril formation. Protein Sci 2004, 13(8):2149–2160. 10.1110/ps.04790604

    PubMed Central  CAS  Article  PubMed  Google Scholar 

  22. 22.

    Lopez De La Paz M, Goldie K, Zurdo J, Lacroix E, Dobson CM, Hoenger A, Serrano L: De novo designed peptide-based amyloid fibrils. Proc Natl Acad Sci USA 2002, 99(25):16052–16057. 10.1073/pnas.252340199

    PubMed Central  Article  PubMed  Google Scholar 

  23. 23.

    Vapnik VN: Statistical Learning Theory. New York: Wiley; 1998.

    Google Scholar 

  24. 24.

    Vapnik VN: The Nature of Statistical Learning Theory. 1st edition. New York: Springer; 1995.

    Google Scholar 

  25. 25.

    Kawashima S, Kanehisa M: AAindex: amino acid index database. Nucleic Acids Res 2000, 28(1):374. 10.1093/nar/28.1.374

    PubMed Central  CAS  Article  PubMed  Google Scholar 

  26. 26.

    Kawashima S, Ogata H, Kanehisa M: AAindex: Amino Acid Index Database. Nucleic Acids Res 1999, 27(1):368–369. 10.1093/nar/27.1.368

    PubMed Central  CAS  Article  PubMed  Google Scholar 

  27. 27.

    Pawar AP, Dubay KF, Zurdo J, Chiti F, Vendruscolo M, Dobson CM: Prediction of "aggregation-prone" and "aggregation-susceptible" regions in proteins associated with neurodegenerative diseases. J Mol Biol 2005, 350(2):379–392. 10.1016/j.jmb.2005.04.016

    CAS  Article  PubMed  Google Scholar 

  28. 28.

    Idicula-Thomas S, Kulkarni AJ, Kulkarni BD, Jayaraman VK, Balaji PV: A support vector machine-based method for predicting the propensity of a protein to be soluble or to form inclusion body on overexpression in Escherichia coli. Bioinformatics 2006, 22(3):278–284. 10.1093/bioinformatics/bti810

    CAS  Article  PubMed  Google Scholar 

  29. 29.

    Lin Y, Lee Y, Wahba G: Support Vector Machines for Classification in Nonstandard Situations. Machine Learning 2002, 46(1–3):191–202. 10.1023/A:1012406528296

    Article  Google Scholar 

  30. 30.

    Pafig[http://www.mobioinfor.cn/pafig]

  31. 31.

    Goldberg DE: Genetic Algorithms in Search, Optimization and Machine Learning. Boston: Addison-Wesley; 1989.

    Google Scholar 

  32. 32.

    LIBSVM[http://www.csie.ntu.edu.tw/~cjlin/]

  33. 33.

    Vihinen M, Torkkila E, Riikonen P: Accuracy of protein flexibility predictions. Proteins 1994, 19(2):141–149. 10.1002/prot.340190207

    CAS  Article  PubMed  Google Scholar 

  34. 34.

    Xia H, Hu B: Feature selection using fuzzy support vector machines. Fuzzy Optim Decis Making 2006, (5):187–192. 10.1007/s10700-006-7336-8

  35. 35.

    Jung E, Kim J, Kim M, Jung DH, Rhee H, Shin JM, Choi K, Kang SK, Kim MK, Yun CH, et al.: Artificial neural network models for prediction of intestinal permeability of oligopeptides. BMC Bioinformatics 2007, 8: 245. 10.1186/1471-2105-8-245

    PubMed Central  Article  PubMed  Google Scholar 

  36. 36.

    Baldi P, Brunak S, Chauvin Y, Andersen CA, Nielsen H: Assessing the accuracy of prediction algorithms for classification: an overview. Bioinformatics 2000, 16(5):412–424. 10.1093/bioinformatics/16.5.412

    CAS  Article  PubMed  Google Scholar 

  37. 37.

    Capriotti E, Calabrese R, Casadio R: Predicting the insurgence of human genetic diseases associated to single point protein mutations with support vector machines and evolutionary information. Bioinformatics 2006, 22(22):2729–2734. 10.1093/bioinformatics/btl423

    CAS  Article  PubMed  Google Scholar 

  38. 38.

    Hua S, Sun Z: Support vector machine approach for protein subcellular localization prediction. Bioinformatics 2001, 17(8):721–728. 10.1093/bioinformatics/17.8.721

    CAS  Article  PubMed  Google Scholar 

  39. 39.

    Tian J, Wu N, Guo X, Guo J, Zhang J, Fan Y: Predicting the phenotypic effects of non-synonymous single nucleotide polymorphisms based on support vector machines. BMC Bioinformatics 2007, 8: 450. 10.1186/1471-2105-8-450

    PubMed Central  Article  PubMed  Google Scholar 

  40. 40.

    Fink AL: Protein aggregation: folding aggregates, inclusion bodies and amyloid. Fold Des 1998, 3(1):R9–23. 10.1016/S1359-0278(98)00002-9

    CAS  Article  PubMed  Google Scholar 

  41. 41.

    Bemporad F, Calloni G, Campioni S, Plakoutsi G, Taddei N, Chiti F: Sequence and structural determinants of amyloid fibril formation. Acc Chem Res 2006, 39(9):620–627. 10.1021/ar050067x

    CAS  Article  PubMed  Google Scholar 

  42. 42.

    Rousseau F, Serrano L, Schymkowitz JW: How evolutionary pressure against protein aggregation shaped chaperone specificity. J Mol Biol 2006, 355(5):1037–1047. 10.1016/j.jmb.2005.11.035

    CAS  Article  PubMed  Google Scholar 

  43. 43.

    Ashburner M, Ball CA, Blake JA, Botstein D, Butler H, Cherry JM, Davis AP, Dolinski K, Dwight SS, Eppig JT, et al.: Gene ontology: tool for the unification of biology. The Gene Ontology Consortium. Nature genetics 2000, 25(1):25–29. 10.1038/75556

    PubMed Central  CAS  Article  PubMed  Google Scholar 

Download references

Acknowledgements

The authors thank Professor Chih-Jen Lin for help using LIBSVM. We would also like to thank Mr Xing Qi for help us revision the paper. This work was supported by the National Natural Science Foundation of China (Grant no. 30770038).

This article has been published as part of BMC Bioinformatics Volume 10 Supplement 1, 2009: Proceedings of The Seventh Asia Pacific Bioinformatics Conference (APBC) 2009. The full contents of the supplement are available online at http://www.biomedcentral.com/1471-2105/10?issue=S1

Author information

Affiliations

Authors

Corresponding author

Correspondence to Ningfeng Wu.

Additional information

Competing interests

The authors declare that they have no competing interests.

Authors' contributions

Jian Tian wrote the code of Pafig. Ningfeng Wu and Yunliu Fan supervised the work. Jian Tian, Ningfeng Wu, and Jun Guo were involved in the preparation of the manuscript. Jian Tian, Ningfeng Wu, Jun Guo and Yunliu Fan read and approved the manuscript.

Electronic supplementary material

Table S1. The fibril-forming peptides

Additional file 1: . The file can be viewed by the software word. (DOC 305 KB)

Rights and permissions

Reprints and Permissions

About this article

Cite this article

Tian, J., Wu, N., Guo, J. et al. Prediction of amyloid fibril-forming segments based on a support vector machine. BMC Bioinformatics 10, S45 (2009). https://doi.org/10.1186/1471-2105-10-S1-S45

Download citation

Keywords

  • Support Vector Machine
  • Fibril
  • Physicochemical Property
  • Reliability Index
  • Receiver Operating Characteristic