BioMed Research International

Volume 2014 (2014), Article ID 935719, 8 pages

http://dx.doi.org/10.1155/2014/935719

Research Article

Predicting the Types of J-Proteins Using Clustered Amino Acids

Pengmian Feng,¹ Hao Lin,² Wei Chen,³ and Yongchun Zuo⁴

¹School of Public Health, Hebei United University, Tangshan 063000, China

²Key Laboratory for Neuroinformation of Ministry of Education, Center of Bioinformatics, School of Life Science and Technology, University of Electronic Science and Technology of China, Chengdu 610054, China

³Department of Physics, School of Sciences, Center for Genomics and Computational Biology, Hebei United University, Tangshan 063000, China

⁴The National Research Center for Animal Transgenic Biotechnology, Inner Mongolia University, Hohhot 010021, China

Received 24 January 2014; Revised 4 March 2014; Accepted 13 March 2014; Published 2 April 2014

Academic Editor: Dong Wang

Copyright © 2014 Pengmian Feng et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract.

J-proteins are molecular chaperones and present in a wide variety of organisms from prokaryote to eukaryote. Based on their domain organizations, J-proteins can be classified into 4 types, that is, Type I, Type II, Type III, and Type IV. Different types of J-proteins play distinct roles in influencing cancer properties and cell death. Thus, reliably annotating the types of J-proteins is essential to better understand their molecular functions. In the present work, a support vector machine based method was developed to identify the types of J-proteins using the tripeptide composition of reduced amino acid alphabet. In the jackknife cross-validation, the maximum overall accuracy of 94% was achieved on a stringent benchmark dataset. We also analyzed the amino acid compositions by using analysis of variance and found the distinct distributions of amino acids in each family of the J-proteins. To enhance the value of the practical applications of the proposed model, an online web server was developed and can be freely accessed.

1. Introduction

J-protein, also known as Hsp40 (heat shock protein 40 kD), is a molecular chaperone protein and is found ubiquitously in both prokaryotes and eukaryotes [1, 2]. J-proteins represent a large family of molecular chaperones and have cooperative functions with Hsp70. Most of the J-proteins contain a “J” domain through which they can interact with and stimulate Hsp70. Based on the structural differences, J-proteins can be classified into four types, that is, Type I, Type II, Type III, and Type IV J-proteins. Type I J-proteins contain an N-terminal J-domain that is separated from the rest of the proteins by a linker “G/F” region (glycine/phenylalanine region) [3, 4]. Distal to G/F region is the zinc-binding cysteine-rich sequence named as “Zinc-finger domain” which distinguishes Type I proteins from other types of J-proteins [4], and Zinc-finger domain is followed by the C-terminal domain [1, 2]. Type II proteins possess all the domains in Type I except the zinc-finger domain [3]. Type III J-proteins contain a C-terminal J-domain but lack both G/F and zinc-finger domains [3]. Type IV, also known as the J-like protein [5], is a group of recently identified proteins that lacks histidine, proline, and aspartate signature motifs in their sequences [4].

By binding Hsp70 and Hsp90, J-proteins play important roles in chaperone cycle regulation and control many physiological functions [4], such as assisting the folding of nascent and damaged proteins, translocation of polypeptides across cellular membranes, and degradation of misfolded proteins [6]. Studies carried out in the past decade have also shown the regulatory roles of J-proteins in cell death. In association with Hsp70, J-proteins not only involve in the folding of caspase-activated DNase which is responsible for the apoptosis-induced DNA fragmentation [7] but also protect the macrophages from nitric-oxide-mediated apoptosis [8]. Gotoh and his colleagues have demonstrated the role of J-protein in the inhibition of Bax translocation to the mitochondria to prevent nitric-oxide-induced cell apoptosis [9]. Kurisu et al. found that MDG1/ERdj4, a member of the human J-protein family, can interact with GRP78/BiP and protect against the cell death induced by endoplasmic reticulum stress in human [10]. The regulation of cell death by J-protein was also reported in plant. Liu and Whitham found that the overexpression of J-protein stimulated the hypersensitive response (HR)-like cell death in soybean [11]. Cancer progressions are also reported to be closely related to J-proteins, but different types of J-proteins play distinct roles [12, 13]. Type I J-protein is tumour promoting, while Type II J-protein acts as tumour suppressors [13]. Therefore, reliably annotating the types of J-proteins is of major importance in order to clarify their distinct biological functions in cell death. However, to the best of our knowledge, there is no computational method for predicting the types of J-proteins.

Keeping these in mind, in the present work, we proposed a model to predict the four functional types of J-proteins based on reduced amino acid alphabet compositions. According to a recent review [14], the rest of the papers are organized as follows: (i) construct a valid benchmark dataset to train and test the predictor; (ii) formulate the samples with an effective mathematical expression that can truly reflect their intrinsic correlation with the target to be predicted; (iii) select a powerful machine learning method to operate the prediction; (iv) perform cross-validation tests to objectively evaluate the anticipated accuracy of the predictor; (v) provide a web server for the prediction method.

2. Materials and Methods

2.1. Dataset

The sequences of J-protein were taken from the HSPIR database at http://pdslab.biochem.iisc.ernet.in/hspir/, which currently contains 3,901 J-protein sequences [15]. To reduce homologous bias, J-proteins that have ≥40% pairwise sequence identity to each other were removed by using the CD-HIT program [16]. By doing so, we obtained a benchmark dataset containing 1,245 J-proteins that were classified into four types: 63 Type I J-proteins, 53 Type II J-proteins, 1,107 Type III J-proteins, and 22 Type IV J-proteins (Table 1). The benchmark dataset can be freely downloaded from http://lin.uestc.edu.cn/server/iJPred/data.

Table 1: Breakdown of the benchmark dataset used in current study.

Total number
Subfamily
Number

1245
Type I J-protein
63

Type II J-protein
53

Type III J-protein
1107

Type IV J-protein
22

2.2. Reduced Amino Acid Alphabet

Based on the physiochemical properties, the 20 native amino acids can be clustered into a smaller number of representative residues called reduced amino acid alphabet (RAAA) [17–19]. Compared with the traditional amino acid composition, RAAA not only simplifies the complexity of protein system but also improves the ability in finding structurally conserved regions and structural similarity of entire proteins.

Recently, a structural alphabet called protein blocks (PBs) was proposed by de Brevern et al. [20, 21] and has been widely used in computational proteomics as indicated in a review [22]. To aid the design of mutations, Etchebest and his colleagues defined a novel type of RAAA based on PBs [23], where the 20 native amino acids can form five different cluster profiles, that is, CP(13), CP(11), CP(9), CP(8), and CP(5) as shown in Table 2. Ever since it was proposed, RAAA has been widely used for protein family classifications [24–27].

Table 2: Scheme for reduced amino acid alphabet based on protein blocks method.

Cluster profiles
Protein blocks method

CP(13)
G-IV-FYW-A-L-M-E-QRK-P-ND-HS-T-C

CP(11)
G-IV-FYW-A-LM-EQRK-P-ND-HS-T-C

CP(9)
G-IV-FYW-ALM-EQRK-P-ND-HS-TC

CP(8)
G-IV-FYW-ALM-EQRK-P-ND-HSTC

CP(5)
G-IVFYW-ALMEQRK-P-NDHSTC

Hence, in the present study, the J-proteins were encoded using the RAAA as formulated by the discrete feature vector P:

where T is the transposing operator and is the occurrence frequency of the th -peptide RAAA and defined as

where is the number of the th-peptide (, 2, or 3) RAAA in a J-protein with length of . For the different cluster profiles (Table 2) and different values of , the vector dimension () in (1) will be different. The corresponding dimensions of reduced amino acid () composition, reduced dipeptide () composition, and reduced tripeptide () composition were listed in Table 3.

Table 3: Feature vector dimension of -peptide composition with different cluster profiles.

-peptide
Cluster profiles

CP(13)
CP(11)
CP(9)
CP(8)
CP(5)

13
11
9
8
5

169
121
81
64
25

2197
1331
729
512
125

2.3. Support Vector Machine (SVM)

SVM is a powerful and popular method for pattern recognition that has been widely used in the realm of bioinformatics [28–41]. The basic idea of SVM is to transform the data into a high dimensional feature space and then determine the optimal separating hyperplane using a kernel function. To handle a multiclass problem, “one-versus-one (OVO)” and “one-versus-rest (OVR)” methods are generally applied to extend the traditional SVM. For a brief formulation of SVM and how it works, see the papers [28, 29].

In the current study, the LIBSVM 2.84 package [42] was used as an implementation of SVM, which can be downloaded from http://www.csie.ntu.edu.tw/~cjlin/libsvm/. The OVO method was employed for making predictions using the popular radial basis function (RBF). The regularization parameter and the kernel width parameter were determined via an optimization procedure using a grid search approach using the fivefold cross-validation. In grid research, the search spaces for parameter and range from 2¹⁵ to and from to with the steps of and 2, respectively.

2.4. Performance Evaluation

The performance of the method was measured in terms of sensitivity (Sn), specificity (Sp), Matthew’s correlation coefficient (MCC), and overall accuracy (OA) defined as follows:

where , , , and represent true positive, true negative, false positive, and false negative of family ; is the number of subsets and equals to 4, while is the number of the total J-proteins in benchmark dataset.

3. Results and Discussion

3.1. Cross-Validation

Three cross-validation methods, namely, subsampling (or K-fold cross-validation) test, independent dataset test, and jackknife test, are often used to evaluate the quality of a predictor [43]. Among the three methods, the jackknife test is deemed the least arbitrary and most objective as elucidated in [44] and hence has been widely recognized and increasingly adopted by investigators to examine the quality of various predictors [31, 34, 45–50]. Accordingly, the jackknife test was used to examine the performance of the model proposed in the current study. In the jackknife test, each sequence in the training dataset is in turn singled out as an independent test sample and all the rule parameters are calculated without including the one being identified.

The jackknife results obtained by the proposed model on the benchmark dataset based on the five different cluster profiles of the tripeptide (i.e., ) case were listed in Table 4. As it can be seen from Table 4, the best success rate of 94.06% was achieved when the predictions were based on CP(8) with a dimension of 512. For comparison, the results of the amino acid (i.e., ) and dipeptide (i.e., ) cases were also calculated and listed in Table 5, from which we can see that none of them has higher success rates than the case of .

Table 4: Results obtained in identifying J-protein functional types with tripeptide case .

Subfamily
Metrics
Feature dimension of for each cluster profile

CP(13)
CP(11)
CP(9)
CP(8)
CP(5)

2197
1331
729
512
125

Type I J-protein
Sn
63.49%
74.60%
77.78%
74.60%
60.31%

Sp
99.56%
98.94%
99.11%
98.76%
98.93%

MCC
0.74
0.76
0.79
0.75
0.66

Type II J-protein
Sn
37.73%
45.28%
39.62%
49.06%
24.53%

Sp
100%
99.31%
99.39%
99.05%
99.56%

MCC
0.60
0.57
0.53
0.57
0.41

Type III J-protein
Sn
99.81%
98.82%
99.09%
98.56%
99.19%

Sp
44.44%
58.78%
55.72%
62.02%
40.00%

MCC
0.63
0.68
0.67
0.69
0.56

Type IV J-protein
Sn
0
27.27%
13.64%
31.81%
4.54%

Sp
100.00%
100.00%
100.00%
100.00%
100.00%

MCC
0
0.52
0.37
0.56
0.21

OA
93.57%
94.06%
93.98%
94.06%
92.36%

Table 5: Results obtained in identifying functional types with (a) single amino acid case and (b) dipeptide case .

(a) For the single amino acid case

Subfamily
Metrics
Feature dimension of for each cluster profile

CP(20)
CP(13)
CP(11)
CP(9)
CP(8)
CP(5)

20
13
11
9
8
5

Type I J-protein
Sn
71.42%
65.08%
68.25%
52.38%
50.79%
22.22%

Sp
98.58%
98.66%
98.66%
98.93%
98.48%
98.03%

MCC
0.71
0.67
0.69
0.60
0.56
0.26

Type II J-protein
Sn
33.96%
30.19%
33.96%
16.98%
16.98%
15.09%

Sp
99.82%
99.21%
99.12%
99.47%
99.56%
99.09%

MCC
0.54
0.42
0.45
0.30
0.31
0.23

Type III J-protein
Sn
98.74%
98.28%
98.10%
99.09%
98.73%
98.19%

Sp
48.12%
42.86%
45.52%
32.31%
31.54%
17.46%

MCC
59.71%
0.53
0.54
0.48
0.45
0.26

Type IV J-protein
Sn
4.54%
0
0
0
0
0

Sp
99.91%
100.00%
100.00%
100.00%
100.00%
100.00%

MCC
0.15
0.53
0
0
0
0

OA
92.93%
91.97%
92.13%
91.48%
91.08%
89.08%

(b) For the dipeptide case

Subfamily
Metrics
Feature dimension of for each cluster profile

CP(20)
CP(13)
CP(11)
CP(9)
CP(8)
CP(5)

400
169
121
81
64
25

Type I J-protein
Sn
74.42%
60.31%
73.02%
60.32%
58.73%
49.20%

Sp
97.58%
98.59%
98.76%
97.71%
98.32%
97.79%

MCC
0.75
0.63
0.73
0.58
0.60
0.5

Type II J-protein
Sn
39.76%
45.23%
39.62%
39.62%
35.84%
28.30%

Sp
94.31%
99.29%
99.48%
99.03%
98.60%
97.99%

MCC
0.57
0.57
0.54
0.49
0.42
0.31

Type III J-protein
Sn
98.88%
98.10%
98.82%
97.74%
98.01%
97.31%

Sp
46.37%
50.74%
51.14%
50.79%
48.80%
40.34%

MCC
60.08%
0.59
0.62
0.57
0.56
0.46

Type IV J-protein
Sn
13.16%
27.27%
0
22.73%
25.00%
9.09%

Sp
99.91%
99.91%
100.00%
100.00%
100.00%
99.91%

MCC
0.13
0.48
0
0.47
0.47
0.24

OA
91.47%
92.93%
93.25%
91.97%
92.04%
91.16%

In our previous study [27], the six HSP families were successfully classified by using the dipeptide of RAAA. But for the classification of the J-protein subfamilies in the present work, the best predictive result was obtained by using the tripeptide of RAAA. Hsps belong to the same family share more sequence identity than that of different families [5]; hence we need more suitable parameters to encode the protein sequences as used in the current study.

3.2. Comparison with Other Methods

Since there is no published work to predict the types of J-proteins, we could not provide the comparison analysis with existing results to confirm that our presented model is superior to other methods. However, for the purpose of comparison, we compared the results of the present model with that of Random Forest and Naïve Bayes using the same optimal features (the reduced tripeptide compositions based on CP(8)). The results of jackknife test on the benchmark dataset for Random Forest and Naïve Bayes are listed in Table 6. It is shown that the accuracy of SVM is higher than that of Random Forest and Naïve Bayes.

Table 6: Comparative result of SVM with other methods for J-protein types classification.

Subfamily
SVM
Random Forest
Naïve Bayes

Sn
SP
MCC
Sn
SP
MCC
Sn
SP
MCC

Type I J-protein
74.60%
98.76%
0.75
14.29%
99.55%
0.29
74.60%
92.17%
0.47

Type II J-protein
49.06%
99.05%
0.57
13.33%
99.82%
0.31
54.72%
94.67%
0.39

Type III J-protein
98.56%
62.02%
0.69
99.73%
12.70%
0.31
88.62%
65.83%
0.43

Type IV J-protein
31.81%
100.00%
0.56
4.55%
100.00%
0.21
13.64%
100.00%
0.37

OA
94.06%
89.96%
85.14%

3.3. Amino Acids Composition Analysis

To provide an overall view, the frequencies of the 20 naive amino acids were compared among the four types of J-proteins using the analysis of variance (ANOVA), and the average amino acid frequency of one type of J-protein with that of another type was further explored and compared using the Fisher’s least significant difference (LSD) test. The result is given in Figure 1, where the green boxes indicate that the frequency differences among different types of J-proteins are not significant, while blue and red boxes indicate that the frequency differences are significant (; LSD test) among different types of J-proteins (see Figure 1 for more details).

Figure 1: Statistical results to show the divergent distributions of the 20 amino acids among the four (I, II, III, and IV) types of J-proteins. The green boxes indicate that the frequency differences among different types of J-proteins are not significant. The blue boxes indicate that the amino acid is significantly enriched (; LSD test) in one type of J-proteins compared with its counterpart. Taking W as an example, the blue box with the coordinate (W, I–IV) indicates that W is enriched in Type I J-proteins compared with Type IV J-proteins. The red boxes indicate that the amino acid is lacking in one type of J-proteins but significantly enriched (; LSD-test) in its counterpart. Also taking W as the example, the two red boxes with the coordinates (W, I–III) and (W, II-III) indicate that W is lacking in both Type I and Type II J-proteins compared with Type III J-proteins, respectively.

We found that, except Asn (N), the frequencies of all the other 19 amino acids are significantly different among the four types of J-proteins. Compared with other three types, Type I J-proteins are enriched in Cys (C), Gly (G), and Thr (T), Type II J-proteins are enriched in Phe (F), Type III J-proteins are enriched in Ala (A) and Leu (L), while Type IV-J proteins are enriched in Met (M), Gln (Q), Glu (E), and Pro (P) but lack Asp (D) and His (H). The lack of D and H residues in Type IV-J proteins leads to their inability to stimulate ATP hydrolysis [5]. Moreover, according to the binomial distribution [51], we also found the overpresented tripeptides in each family and listed them in Supporting Table S1 in Supplementary Material available online at http://dx.doi.org/10.1155/2014/935719, where the over-presented tripeptides with their confidence levels are provided. These results indicate that the distinct distributions of amino acids in the four types of J-proteins may account for their distinct functions in biological processes.

3.4. Web Server Guide

To enhance the value of the practical applications of the proposed model and for the convenience of the vast majority of experimental scientists, an online predictor was developed. The step-by-step guide on how to use it is provided as follows.

(1)Open the web server at http://lin.uestc.edu.cn/server/Jpred and you will see the top page as shown in Figure 2. Click on the Read Me button to see a brief introduction about the predictor and the caveat when using it, and click on the Data button to download the benchmark datasets used to train and test the predictor. The relevant papers that document the algorithm of the predictor can be found by clicking on the Citation button.(2)Either type or copy/paste the query J-protein sequence into the input box at the center of Figure 2. The input protein sequence should be in the FASTA format that can be seen by clicking on the Example button right above the input box.(3)Click on the Submit button to see the predicted result. For example, if you use the four query J-protein sequences in the Example window as the input, after clicking the Submit button, you will obtain the results: the outcome for the 1st query sample is “Type I J-protein;” the outcome for the 2nd query sample is “Type II J-protein;” the outcome for the 3rd query sample is “Type III J-protein;” the outcome for the 4th query sample is “Type IV J-protein.”

Figure 2: A semiscreenshot to show the top page of the web server. It is available at http://lin.uestc.edu.cn/server/Jpred.

4. Conclusion

Cell death is a common phenomenon in developmental processes or in normal physiological conditions and is induced by an array of extra- or intracellular stimuli [7]. However, organisms are equipped with their own physiological defense to cope with environmental stress in order to prevent or induce cell death depending upon the severity of the stress [7]. In mammalian cells, the stress response involves the induction of Hsps, such as Hsp70 and Hsp90. By interacting with J-proteins, these Hsps play pivotal roles in cell death regulations. Since J-proteins act as intermediates, the analysis of J-proteins functions is urgent in order to clarify the regulatory roles of Hsps in cell death.

Based on combination of whole-genome analyses and biochemical evidences, a large number of J-proteins have been identified [6]. However, the exact roles for many of the J-proteins are far from being understood [2, 52]. In order to understand its biological functions, it is highly desirable to know which family a given J-protein belongs to.

By encoding the sequences using the reduced amino acid alphabet information, a predictor was developed to identify the four different families of J-proteins in the present work. To enhance the value of the practical applications of the proposed model and for the convenience of the experimental scientists, an online web server was provided and can be freely accessed at http://lin.uestc.edu.cn/server/Jpred. We hope that the present model will be helpful for scientists who focus on J-proteins and will provide novel insights into the research of cell death.

Conflict of Interests

There is no conflict of interests with any financial organization regarding this paper.

Acknowledgments

The authors wish to express their gratitude to executive editor and three anonymous reviewers whose constructive comments were very helpful in strengthening the presentation of this paper. This work was supported by the National Nature Scientific Foundation of China (nos. 61100092 and 61202256), Nature Scientific Foundation of Hebei Province (no. C2013209105), Foundation of Science and Technology Department of Hebei Province (no. 132777133), and the Fundamental Research Funds for the Central Universities (ZYGX2013J102).

References

A. J. Caplan, D. M. Cyr, and M. G. Douglas, “Eukaryotic homologues of Escherichia coli dnaJ: a diverse protein family that functions with HSP70 stress proteins,” Molecular and Cellular Biology, vol. 4, no. 6, pp. 555–563, 1993.

X. B. Qiu, Y. M. Shao, S. Miao, and L. Wang, “The diversity of the DnaJ/Hsp40 family, the crucial partners for Hsp70 chaperones,” Cellular and Molecular Life Sciences, vol. 63, no. 22, pp. 2560–2570, 2006.

M. E. Cheetham and A. J. Caplan, “Structure, function and evolution of DnaJ: conservation and adaptation of chaperone function,” Cell Stress Chaperones, vol. 3, no. 1, pp. 28–36, 1998.

V. B. Rajan and P. D'Silva, “Arabidopsis thaliana J-class heat shock proteins: cellular stress sensors,” Functional & Integrative Genomics, vol. 9, no. 4, pp. 433–446, 2009.

P. Walsh, D. Bursać, Y. C. Law, D. Cyr, and T. Lithgow, “The J-protein family: modulating protein assembly, disassembly and translocation,” EMBO Reports, vol. 5, no. 6, pp. 567–571, 2004.

E. A. Craig, P. Huang, R. Aron, and A. Andrew, “The diverse roles of J-proteins, the obligate Hsp70 co-chaperone,” Reviews of Physiology, Biochemistry and Pharmacology, vol. 156, pp. 1–21, 2006.

A. S. Sreedhar and P. Csermely, “Heat shock proteins in the regulation of apoptosis: new strategies in tumor therapy—a comprehensive review,” Pharmacology & Therapeutics, vol. 101, no. 3, pp. 227–257, 2004.

T. Gotoh, K. Terada, and M. Mori, “Hsp70-DnaJ chaperone pairs prevent nitric oxide-mediated apoptosis in RAW 264.7 macrophages,” Cell Death & Differentiation, vol. 8, no. 4, pp. 357–366, 2001.

T. Gotoh, K. Terada, S. Oyadomari, and M. Mori, “hsp70-DnaJ chaperone pair prevents nitric oxide- and CHOP-induced apoptosis by inhibiting translocation of Bax to mitochondria,” Cell Death & Differentiation, vol. 11, no. 4, pp. 390–402, 2004.

J. Kurisu, A. Honma, H. Miyajima, S. Kondo, M. Okumura, and K. Imaizumi, “MDG1/ERdj4, an ER-resident DnaJ family member, suppresses cell death induced by ER stress,” Genes to Cells, vol. 8, no. 2, pp. 189–192, 2003.

J. Z. Liu and S. A. Whitham, “Overexpression of a soybean nuclear localized type-III DnaJ domain -containing HSP40 reveals its roles in cell death and disease resistance,” Plant Journal, vol. 74, no. 1, pp. 110–121, 2013.

A. Mitra, L. A. Shevde, and R. S. Samant, “Multi-faceted role of HSP40 in cancer,” Clinical & Experimental Metastasis, vol. 26, no. 6, pp. 559–567, 2009.

J. N. Sterrenberg, G. L. Blatch, and A. L. Edkins, “Human DNAJ in cancer and stem cells,” Cancer Letters, vol. 312, no. 2, pp. 129–142, 2011.

K.-C. Chou, “Some remarks on protein attribute prediction and pseudo amino acid composition,” Journal of Theoretical Biology, vol. 273, no. 1, pp. 236–247, 2011.

R. K. Ratheesh, S. N. Nagarajan, P. A. Arunraj et al., “HSPIR: a manually annotated heat shock protein information resource,” Bioinformatics, vol. 28, no. 21, pp. 2853–2855, 2012.

W. Li and A. Godzik, “Cd-hit: a fast program for clustering and comparing large sets of protein or nucleotide sequences,” Bioinformatics, vol. 22, no. 13, pp. 1658–1659, 2006.

P. D. Thomas and K. A. Dill, “An iterative method for extracting energy-like quantities from protein structures,” Proceedings of the National Academy of Sciences of the United States of America, vol. 93, no. 21, pp. 11628–11633, 1996.

L. A. Mirny and E. I. Shakhnovich, “Universally conserved positions in protein folds: reading evolutionary signals about stability, folding kinetics and function,” Journal of Molecular Biology, vol. 291, no. 1, pp. 177–196, 1999.

A. D. Solis and S. Rackovsky, “Optimized representations and maximal information in proteins,” Proteins, vol. 38, no. 2, pp. 49–164, 2000.

A. G. de Brevern, “New assessment of a structural alphabet,” In Silico Biology, vol. 5, no. 3, pp. 283–289, 2005.

A. G. de Brevern, C. Etchebest, and S. Hazout, “Bayesian probabilistic approach for predicting backbone structures in terms of protein blocks,” Proteins, vol. 41, no. 3, pp. 271–287, 2000.

A. P. Joseph, G. Agarwal, S. Mahajan et al., “A short survey on protein blocks,” Biophysical Reviews, vol. 2, no. 3, pp. 137–147, 2010.

C. Etchebest, C. Benros, A. Bornot, A.-C. Camproux, and A. G. de Brevern, “A reduced amino acid alphabet for understanding and designing protein adaptation to mutation,” European Biophysics Journal, vol. 36, no. 8, pp. 1059–1069, 2007.

Y. C. Zuo and Q. Z. Li, “Using reduced amino acid composition to predict defensin family and subfamily: Integrating similarity measure and structural alphabet,” Peptides, vol. 30, no. 10, pp. 1788–1793, 2009.

W. Chen, P. M. Feng, and H. Lin, “Prediction of ketoacyl synthase family using reduced amino acid alphabets,” Journal of Industrial Microbiology and Biotechnology, vol. 39, no. 4, pp. 579–584, 2011.

Y. L. Chen, Q. Z. Li, and L. Q. Zhang, “Using increment of diversity to predict mitochondrial proteins of malaria parasite: integrating pseudo-amino acid composition and structural alphabet,” Amino Acids, vol. 42, no. 4, pp. 1309–1316, 2010.

P. M. Feng, W. Chen, H. Lin, and K. C. Chou, “iHSP-PseRAAAC: identifying the heat shock protein families using pseudo reduced amino acid alphabet composition,” Analytical Biochemistry, vol. 442, no. 1, pp. 118–125, 2013.

K. C. Chou and Y. D. Cai, “Using functional domain composition and support vector machines for prediction of protein subcellular location,” The Journal of Biological Chemistry, vol. 277, no. 48, pp. 45765–45769, 2002.

Y. D. Cai, G. P. Zhou, and K. C. Chou, “Support vector machines for predicting membrane protein types by using functional domain composition,” Biophysical Journal, vol. 84, no. 5, pp. 3257–3263, 2003.

W. Chen and H. Lin, “Prediction of midbody, centrosome and kinetochore proteins based on gene ontology information,” Biochemical and Biophysical Research Communications, vol. 401, no. 3, pp. 382–384, 2010.

M. Hayat and A. Khan, “MemHyb: predicting membrane protein types by hybridizing SAAC and PSSM,” Journal of Theoretical Biology, vol. 292, pp. 93–102, 2012.

H. Lin and H. Ding, “Predicting ion channels and their types by the dipeptide mode of pseudo amino acid composition,” Journal of Theoretical Biology, vol. 269, no. 1, pp. 64–69, 2011.

X. Xiao, P. Wang, and K.-C. Chou, “iNR-physchem: a sequence-based predictor for identifying nuclear receptors and their subfamilies via physical-chemical property matrix,” PLoS ONE, vol. 7, no. 2, Article ID e30869, 2012.

W. Chen, P. M. Feng, H. Lin, and K. C. Chou, “iRSpot-PseDNC: identify recombination spots with pseudo dinucleotide composition,” Nucleic Acids Research, vol. 41, no. 6, p. e68, 2013.

P. M. Feng, H. Lin, and W. Chen, “Identification of antioxidants from sequence information using Naive Bayes,” Computational and Mathematical Methods in Medicine, vol. 2013, Article ID 567529, 5 pages, 2013.

B. Liu, X. Wang, Q. Zou, Q. Dong, and Q. Chen, “Protein remote homology detection by combining chou's pseudo amino acid composition and profile-based protein representation,” Molecular Informatics, vol. 32, no. 9-10, pp. 775–782, 2013.

B. Liu, X. Wang, Q. Chen, Q. Dong, and X. Lan, “Using amino acid physicochemical distance transformation for fast protein remote homology detection,” PLoS ONE, vol. 7, no. 9, Article ID e46633, 2012.

P. M. Feng, H. Ding, W. Chen, and H. Lin, “Naive Bayes classifier with feature selection to identify phage virion proteins,” Computational and Mathematical Methods in Medicine, vol. 2013, Article ID 530696, 6 pages, 2013.

B. Liu, X. Wang, L. Lin, B. Tang, Q. Dong, and X. Wang, “Prediction of protein binding sites in protein structures using hidden Markov support vector machine,” BMC Bioinformatics, vol. 10, article 381, 2009.

B. Liu, X. Wang, L. Lin, Q. Dong, and X. Wang, “Exploiting three kinds of interface propensities to identify protein binding sites,” Computational Biology and Chemistry, vol. 33, no. 4, pp. 303–311, 2009.

B. Liu, X. Wang, L. Lin, Q. Dong, and X. Wang, “A discriminative method for protein remote homology detection and fold recognition combining Top-n-grams and latent semantic analysis,” BMC Bioinformatics, vol. 9, article 510, 2008.

C. C. Chang and C. J. Lin, “LIBSVM: a library for support vector machines,” 2001, http://www.csie.ntu.edu.tw/~cjlin/libsvm

K. C. Chou and C. T. Zhang, “Prediction of protein structural classes,” Critical Reviews in Biochemistry and Molecular Biology, vol. 30, no. 4, pp. 275–349, 1995.

K. C. Chou and H. B. Shen, “Cell-PLoc: a package of Web servers for predicting subcellular localization of proteins in various organisms,” Nature Protocols, vol. 3, no. 2, pp. 153–162, 2008.

M. Esmaeili, H. Mohabatkar, and S. Mohsenzadeh, “Using the concept of Chou's pseudo amino acid composition for risk type prediction of human papillomaviruses,” Journal of Theoretical Biology, vol. 263, no. 2, pp. 203–209, 2010.

S. Mei, “Predicting plant protein subcellular multi-localization by Chou's PseAAC formulation based multi-label homolog knowledge transfer learning,” Journal of Theoretical Biology, vol. 310, pp. 80–87, 2012.

K. C. Chou, Z. C. Wu, and X. Xiao, “iLoc-Euk: a multi-label classifier for predicting the subcellular localization of singleplex and multiplex eukaryotic proteins,” PLoS ONE, vol. 6, no. 3, Article ID e18258, 2011.

H. Mohabatkar, “Prediction of cyclin proteins using Chou's pseudo amino acid composition,” Protein and Peptide Letters, vol. 17, no. 10, pp. 1207–1214, 2010.

K. C. Chou, Z. C. Wu, and X. Xiao, “ILoc-Hum: using the accumulation-label scale to predict subcellular locations of human proteins with both single and multiple sites,” Molecular BioSystems, vol. 8, no. 2, pp. 629–641, 2012.

H. Mohabatkar, M. Mohammad Beigi, and A. Esmaeili, “Prediction of GABAA receptor proteins using the concept of Chou's pseudo-amino acid composition and support vector machine,” Journal of Theoretical Biology, vol. 281, no. 1, pp. 18–23, 2011.

C. Ding, L. F. Yuan, S. H. Guo, H. Lin, and W. Chen, “Identification of mycobacterial membrane proteins and their types using over-represented tripeptide compositions,” Journal of Protemics, vol. 77, pp. 321–328, 2012.

K. Ohtsuka and M. Hata, “Mammalian HSP40/DNAJ homologs: cloning of novel cDNAs and a proposal for their classification and nomenclature,” Cell Stress and Chaperones, vol. 5, no. 2, pp. 98–112, 2000.