Nothing Special   »   [go: up one dir, main page]

Skip to main content

Accurate classification of membrane protein types based on sequence and evolutionary information using deep learning

Abstract

Background

Membrane proteins play an important role in the life activities of organisms. Knowing membrane protein types provides clues for understanding the structure and function of proteins. Though various computational methods for predicting membrane protein types have been developed, the results still do not meet the expectations of researchers.

Results

We propose two deep learning models to process sequence information and evolutionary information, respectively. Both models obtained better results than traditional machine learning models. Furthermore, to improve the performance of the sequence information model, we also provide a new vector representation method to replace the one-hot encoding, whose overall success rate improved by 3.81% and 6.55% on two datasets. Finally, a more effective model is obtained by fusing the above two models, whose overall success rate reached 95.68% and 92.98% on two datasets.

Conclusion

The final experimental results show that our method is more effective than existing methods for predicting membrane protein types, which can help laboratory researchers to identify the type of novel membrane proteins.

Background

Protein is an important component of all cells and tissues in the human body and is the material basis of life [1–5]. Membrane proteins represent one important protein type that is rich in function, they participate in many important reactions of the cell, including transporting substances into and out of cells as a carrier, acting as a specific receptor for hromones, and carrying out cellular recognition functions, as well as being responsible for signal transduction and cell-cell interactions [6]. Among the genomes that have been completely sequenced, the number of membrane proteins accounts for 30% [7]. Moreover, membrane proteins are of particular importance in drug therapy, as they act as the targets for many drugs. Currently, the effects more than 50% of drugs on the market are exerted by the actions of membrane proteins [8]. Since the function of membrane protein is closely related to its type, the prediction of membrane protein types can contribute to research in the field of bioinformatics.

According to their functions, membrane proteins can be classified into three classes: integral, peripheral and lipid-anchored [9]. Based on the direct interaction relationship between membrane proteins and lipid bilayers, the three classes can be further extended into eight basic types: (1) type I membrane proteins, (2) type II membrane proteins, (3) type III membrane proteins (4) type IV membrane proteins, (5) multipass transmembrane proteins, (6) lipid chain-anchored membrane proteins, (7) GPI-anchored membrane proteins, and (8) peripheral membrane proteins. Among them, Types I, II, III and IV are single-pass transmembrane proteins, and detailed descriptions of their differences are given in [9].

To help laboratory researchers discover the type of novel membrane protein, various computation methods are proposed for membrane protein type recognition. Many of these approaches incorporate machine learning algorithms and statistical analysis techniques, such as k-nearest neighbor (KNN) [10], the naive Bayesian model (NBM) [11], support vector machines (SVM) [12–14], random forests (RF) [15], probabilistic neural network (PNN) [16] and hidden Markov models [7].

The most popular feature extraction methods for predicting membrane protein types are based on sequence information. The pseudo-amino acid composition (PseAAC) [17] method incorporates information about sequence order. Local amino acid composition (LAAC), local dipeptide composition (LDC), global descriptor (GD), Lempel-Ziv complexity (LZC), autocorrelation descriptor (AD), sequence-order descriptor (SD) and Hilbert-Huang transform (HHT) are proposed based on amino acid classification and physicochemical properties [18], and these methods are partially applied in our previous research [19]. The peptide composition method, such as dipeptide composition (DipC) [20] and tripeptide composition (TipC) [21] are also powerful sequence information-based feature extraction methods.

The evolutionary information, which mainly refers to the position-specific scoring matrix (PSSM), is also widely used to predict membrane protein types [20, 22]. To extract rich evolutionary information from PSSM, many researchers have proposed their own methods, such as the reduced position-specific score matrix (RPSSM) [23], evolutionary difference position-specific score matrix (EDPSSM) [24], tri-gram position-specific score matrix (TriPSSM) [25], k-separated-bigrams position-specific score matrix (KPSSM) [26], correlation position-specific scoring matrix (CoPSSM) [27], and pseudo position-specific score matrix (PsePSSM) [9]. These methods are all obtained through complex feature engineering.

In this paper, we improve the performance in membrane protein type prediction. The main contributions of the paper are summarized as follows: First, we present two deep neural network (DNN) models to process the sequence information and evolutionary information separately. These models both achieve better performance than traditional machine learning models. Second, by using the DNN model with convolutional and recurrent layers, we can remove the burden of feature engineering and the reliance on domain experts. Third, we provide a new vector representation based on autoencoder and physicochemical property indexes. Numerous experimental results prove that the new vector has more powerful representation ability than one-hot encoding. Finally, a more effective model is constructed by fusing the above two models with the ensemble learning.

The DNN model has shown its superiority in the field of bioinformatics [28–37]. However, to the best of our knowledge, this paper is the first to propose DNN models for the prediction of membrane protein types. The two DNN models proposed in this paper can process the sequence information and evolutionary information, separately. When processing sequence information, we use a 1D convolutional layer [38] and bidirectional long short term-memory (Bi-LSTM) layers [39], and when processing evolutionary information, a 2D convolutional layer and CapsNet layer [40] are adopted. Among these layers, compared to convolutional network, CapsNet is a new network that can extract local position information.

Methods

Datasets

In this work, we use two datasets to evaluate the performance of our method, hereafter referred to as Dataset 1 and Dataset 2. Dataset 1 consists of two benchmark sets, the training set and the testing set, which are both from paper [41] and used in previous studies [9, 10, 12, 13, 20]. Dataset 2 is constructed by Chen [12]; it also consists of two benchmark sets. Table 1 shows the concrete distribution of Dataset 1 and Dataset 2.

Table 1 The distribution of samples for Dataset 1 and Dataset 2

Architecture of the proposed DNN

To improve the performance of the prediction of membrane protein types, we build the DNN model with a keras framework (http://www.keras.io). Considering the respective characteristics of sequence information and evolutionary information, we build two different models to process them, which are hereafter referred to as the sequence information model and evolutionary information model, respectively.

Sequence information model

In this model, the input sequences are converted into numerical vectors of fixed length by truncating the long sequence and padding the short sequence with ‘0’. The length of the sequence ranges from 50 to 5000 in our dataset. Although the length of the sequence varies greatly, sequences with a length of less than 1500 account for almost 98% of the total dataset. Therefore, taking the utilization of information and the computational complexity into account, we decided to set the fixed length to 1500. Each of the basic 20 amino acids is converted into a number of 0-19. Ultimately, to prevent different encoding from affecting the performance of model, one-hot encoding is used before inputting the vector into the model.

The architecture of the sequence information model, as shown in Fig. 1, consists of one 1D convolutional layer (filters: 256, kernel__size: 15, strides: 10, padding: ‘same’, activation: ‘relu’), two Bi-LSTM layers (units: 128, dropout: 0.5 and rest default setting) and one fully connected layer (softmax). The 1D convolutional layer can reduce the complexity of the Bi-LSTM layer that reduces the shape of the output matrix from 1500 ×256 tensor to 150 ×256. Moreover, the 256 filters in the layer can also extract rich information from the sequence. Next, two Bi-LSTM layers are applied, which can identify features separated by large gaps. In previous studies, a normal LSTM layer is mainly used to process sequence information [29, 30]. Compared to normal LSTM, the Bi-LSTM layer we used can take advantage of both historical input and future information, which performs better for our task. Each Bi-LSTM unit comprises three gates, which are shown in Fig. 1. Their calculation process is as follows:

$$ f_{t}=\sigma (W_{f}\cdot [h_{t-1},x_{t}]+b_{f})\qquad\text{forget\ gate} $$
(1)
Fig. 1
figure 1

The sequence information model uses 1D Conv and Bi-LSTM layers. First, the membrane protein sequences are encoded into a matrix of size 1500 ×20. Then, these matrixes are fed to the 1D Conv layer (filters: 256, kernel size: 15, strides: 10, padding: same, activation: relu). Next, two Bi-LSTM layers (units: 128, dropout: 0.5 and the remainder as default settings) are used to identify features separated by large gaps. Last, the output from the Bi-LSTM layer is passed through a final dense layer that uses a softmax function to obtain the probability value belonging to each membrane protein type

$$ i_{t}=\sigma (W_{i}\cdot [h_{t-1},x_{t}]+b_{i})\qquad\text{input\ gate} $$
(2)
$$ \tilde{C}_{t}=tanh(W_{C}\cdot [h_{t-1},x_{t}]+b_{C})\qquad\text{candidate\ cell\ states} $$
(3)
$$ C_{t}=f_{t}*C_{t-1}+i_{t}*\tilde{C}_{t}\qquad\text{updated\ cell\ states} $$
(4)
$$ o_{t}=\sigma (W_{o}\cdot [h_{t-1},x_{t}]+b_{o})\qquad\text{output\ gate} $$
(5)
$$ h_{t}=o_{t}*tanh(C_{t}))\qquad\text{hidden\ cell\ states} $$
(6)

where Wf, Wi, WC and Wo denote the weight matrix of the forget gate, input gate, candidate cell states and output gate, respectively, and bf, bi, bC, and bo denote the bias of the forget gate, input gate, candidate cell states and output gate, respectively. The forget gate decides which information is discarded from the cell state, where ht−1 represents the output of the previous cell, and xt represents the input of the current cell. σ represents the sigmod function. The input gate decides how much new information is added to the cell states, and it uses two steps to accomplish this task: First, a sigmoid layer determines which information needs to be updated; a tanh layer generates a vector \(\tilde {C}_{t}\), which is the alternative content to update. Then, we combine two parts to update the state of the cell. The output gate will determine what value to output. First, we run a sigmoid layer to determine which part of the cell state will be output. Next, we process the cell state through tanh (obtaining a value between -1 and 1) and multiply it by the output of the sigmoid gate. Eventually, we will only output the part of our output that we determined.

Then, the output from Bi-LSTM layer is passed through a final dense layer that used a softmax function to obtain the probability value belonging to each membrane protein type. We compile our keras model with the ‘Adam’ optimizer (lr: 0.001, beta__1: 0.9, beta__2: 0.999, decay=0). Finally, we find that adding the BatchNormalization layers (default setting) after three mainly layers in our model helps to accelerate convergence.

Evolutionary information model

In this paper, evolutionary information mainly refers to the PSSM, which was first proposed by Jones for predicting the secondary structure of proteins [42]. The PSSM can discover protein sequences that have evolved relationships with search sequences. It is expressed as follows:

$$ PSSM= \left[\begin{array}{cccc} M_{1\rightarrow 1} &M_{1\rightarrow 2} &\cdots &M_{1\rightarrow 20} \\ M_{2\rightarrow 1} &M_{2\rightarrow 2} &\cdots &M_{2\rightarrow 20} \\ \vdots &\vdots &\ddots &\vdots \\ M_{L\rightarrow 1} &M_{L\rightarrow 2} &\cdots &M_{L\rightarrow 20} \end{array}\right] $$
(7)

where L is the length of the sequence, and Mi→j denotes the score of the amino acid in the i-th position of the amino acid being mutated to j-th position of the amino acid during the evolution process. Because membrane protein sequences vary in size, their PSSM shapes are also different. Similar to the sequence information model, we also truncate the long sequence and pad the short sequence with ‘0’ to convert the PSSM matrix into a 1500 ×20 matrix.

The architecture of evolutionary information model is shown in Fig. 2. It consists of two 2D convolutional layers, two average pooling layers, one PrimaryCaps and one fully connected layer. The first 2D convolutional layer (filters: 256, kernel__size: (5,5), strides: 1, padding: ‘valid’, activation: ‘relu’) is used to extract rich and effective features from the PSSM. Next, the first 2D average pooling layer (pool__size: (2,2) and the remainder as default settings) is used to retain the main features while reducing the parameters and calculations of the next layer to prevent overfitting. Although paper [40] argued that the pooling layer may throw away information about the precise position of the entity within the region, we find that the pooling layer and CapsNet work together to obtain the best results in our task. One possible reason is that the pooling layer can improve the robustness of extracted features. Then, the second 2D convolutional layers (filters: 128, kernel__size: (5,5), strides: 1, padding: ‘valid’, activation: ‘relu’) and the second 2D average pooling layer (pool__size: (2,2) and the remainder as default settings) are reused to obtain more advanced features. Moreover, each convolutional layer has a dropout technique with the rate 0.5 to prevent the model from overfitting. The first four layers are designed to increase the representation power of CapsNet.

Fig. 2
figure 2

The evolutionary information model uses 2D Conv and CapsNet layers. First, the PSSM matrix is processed as the matrix with the same shape. These matrixes are fed to the 2D Conv layer (filters: 256, kernel size: (5,5), strides: 1, padding: valid, activation: relu) and average pooling layer (pool size: (2,2) and the remainder as default settings) to extract the rich and effective features. Next, the 2D Conv layer(filters: 128, kernel__size: (5,5), strides: 1, padding: ’valid’, activation: ’relu’) and average pooling layer(pool__size: (2,2)) are reused to extract more advanced features. The PrimaryCaps layer is the convolutional capsule layer, which has size 1 ×1 convolution kernels and 4 channels of 16D capsules. The ProteinCaps layer has eight 16D capsules to represent one of membrane protein types. Finally, the L2-norm of each capsule vector is calculated to indicate the probability of each type

The PrimaryCaps layer contains 4 primary capsules that accept the basic features detected by the first four layers and generate a combination of features. The 4 main capsules of this layer are essentially similar to the convolutional layer. Among them, the 1 ×1 convolution kernel used in our study not only achieves cross-channel interaction and information integration, but also reduces the number of convolution kernel channels. Each capsule applies sixteen 1 ×1×64 convolution kernels to the 372 ×2×128 input tensor, thus generating a 372 ×2×64 output tensor. With a total of 4 capsules, the output is a 372 ×2×4×16 tensor. Since the length of a capsule represents the probability that the entity presented, the CapsNet layer needs a new activation function, which is called the squashing function. It is a novel nonlinear activation function that accepts an input vector and then compresses its length to [0,1] as output, which is calculated as follows:

$$ V_{j}=\frac{\left \| s_{j} \right \|}{1+\left \| s_{j} \right \|^{2}}\frac{s_{j}}{\left \| s_{j} \right \|} $$
(8)

where Vj and sj are the vector output and input of the CapsNet, respectively.

The last layer contains eight membrane protein type capsules, one for each type of membrane protein. The calculation between the PrimaryCaps layer and ProteinCaps layer is shown in Fig. 3. Among them, ui represents the i-th capsule in PrimaryCaps. Wij is the weight matrix that represents the spatial relationship and other important relationships between low-level capsules and high-level capsules. There are eight capsules (Vj,j∈[1,2,⋯,8]) in ProteinCaps, each of which receives the inputs from all capsule outputs in PrimaryCaps. Vj is calculated by the weighted sum of \(\tilde {u}_{ji}\) and then is passed through the squashing function. Here, the weights cij are obtained by an iterative dynamic routing process that can be found in [40].

Fig. 3
figure 3

Computation between the PrimaryCaps and ProteinCaps layers.There are 2976 16D capsules (each ui is an 16D vector) in PrimaryCaps. Each \(\hat {u_{j|i}}\) is produced by multiplying ui by a weight matrix Wi,j (16 ×16). Capsule Vj (16D vector and j∈ [1, 8]) in ProteinCaps is produced by a weighted sum over all \(\hat {u_{j|i}}\) and the squashing nonlinear activation function. The parameter ci,j is determined by the iterative dynamic routing process

The output of the ProteinCaps layer has eight 16D vectors. During training, for each training sample, the loss value is calculated as follows:

$$ L_{c}\,=\,T_{c}\text{max}(0,m^{+}\!-\left \| V_{c} \right \|))^{2}+\lambda (1-T_{c})\text{max}(0,\left \| V_{c} \right \|)-m^{-})^{2} $$
(9)

where the value of Tc is decided by the correct label. If the correct label corresponds to the membrane protein type in ProteinCaps, then Tc=1; otherwise, Tc=0. For the value of the other hyperparameters m+, m− and λ, we used the default values 0.9,0.1 and 0.5, respectively, which are suggested in [40].Then, eight loss values are added to obtain the final loss.

Finally, we use the ‘Adam’ optimizer method (lr: 0.0001, beta__1: 0.9, beta__2: 0.999, decay=0). The suggested settings, three routing iterations and the margin loss function, are used in our work.

Fusion information model

Many studies have proved that the fusing information method has better results than the single information method [12, 20, 43–45]. Therefore, we also present a fusion information model based on ensemble learning. Among these methods, the ensemble learning method that we adopted is a stacking algorithm, and its details can be found in our previous work [19]. Figure 4 shows the complete flow chart of the stacking algorithm, where the meta-classifier that we used is multiple logistic regression (MLR).

Fig. 4
figure 4

The flow of the stacking algorithm. Model_1 refers to the sequence information model, and Model_2 refers to the evolutionary information model. The class probability outputs of two models are combined to serve as a new training set to train the meta-classifier (MLR). When testing unknown membrane protein samples, we take the result of the meta-classifier as the final output

Research has shown that using the class probability output of the base classifier as the input of the meta classifier performs better [46]. Then, we use P(x) as the output of the base classifier, it can be represented as follows:

$${\begin{aligned} P(x)&=\left(P^{\mathrm{T}}_{1}, P^{\mathrm{T}}_{2}\right) \\ &= \left(\underbrace{P^{1}_{1}, P^{1}_{2}, P^{1}_{3}, P^{1}_{4}, P^{1}_{5}, P^{1}_{6}, P^{1}_{7}, P^{1}_{8}}_{{Sequence\_informtaion\_model}}, \underbrace{P^{2}_{1}, P^{2}_{2}, P^{2}_{3}, P^{2}_{4}, P^{2}_{5}, P^{2}_{6}, P^{2}_{7}, P^{2}_{8}}_{Evolutionary\_informtaion\_model}\right) \end{aligned}} $$

Model training

In most experiments, the DNN models are trained using identical training strategies. The dataset is divided into three parts, the training set, validation set and testing set, which play different roles in our work. The training set is used to train the model, and the validation set is used to adjust parameters and select the best model. The testing set is used to evaluate the performance of the model at the end. In our work, 20% of samples are separated from the original training set as a validation set. Then, the distributions of three sets in two datasets are shown in Table 2. All DNN models are implemented using keras 2.1.2. Model training and testing are performed on a personal computer equipped with an Nvidia GTX 1060 GPU.

Table 2 The distribution of three data in two datasets

Model evaluation

We use sensitivity (Se), specificity (Sp), accuracy (ACC) and overall success rate (OSR) to evaluate the classification performance. They are defined as follows:

$$ Se_{i}=TP_{i}/(TP_{i}+FN_{i}) $$
(10)
$$ Sp_{i}=TN_{i}/(TN_{i}+FP_{i}) $$
(11)
$$ ACC_{i}=\frac{TP_{i}+TN_{i}}{TP_{i}+TN_{i}+FN_{i}+FP_{i}} $$
(12)
$$ OSR=\frac{\sum_{i=1}^{k}TP_{i}}{\sum_{i=1}^{k}(TP_{i}+FN_{i})} $$
(13)

where TP, TN, FP and FN represent true positive, true negative, false positive and false negative, respectively.

The Matthews correlation coefficient (MCC) is also used in our work. It is generally considered to be a more balanced indicator, even if the sample content of the two categories varies widely [47]. The MCC is essentially a correlation coefficient between the actual classification and the predicted classification. Its value range is [-1, 1]. A value of 1 indicates a perfect prediction for the subject. A value of 0 indicates that the predicted result is not as good as the result of random prediction, and -1 means that the predicted classification is completely inconsistent with the actual classification. The MCC is defined as follows:

$$ {\begin{aligned} MCC_{i}=\frac{TP_{i}\times TN_{i}-FP_{i} \times FN_{i}}{\sqrt{(TP_{i}+FN_{i})(TP_{i}+FP_{i})(TN_{i}+FP_{i})(TN_{i}+FN_{i})}} \end{aligned}} $$
(14)

The construction of vector representation

One-hot encoding ignores the possible relationship between amino acids. For example, alanine and cysteine have a certain relationship in their physicochemical properties, yet ‘0’ and ‘1’ do not provide much informtaion about this relationship. Moreover, we usually need more data to train because the amino acids are stored as sparse matrixes. Using a vector representation method can effectively solve the above problems. Continuous bag of words (CBOW) and Skim-Gram [48] are two popular vector representation methods that are widely used in the field of natural language processing (NLP). Limited by the size of protein samples and the characteristics of the protein sequence, these methods are no longer applicable. In this paper, we provide a new vector representation method based on autoencoder and the physicochemical properties of amino acids. This method can achieve a more compact representation of input symbols and yield semantically similar symbols close to each other in vector space.

The AAindex is a database of numerical indexes representing various physicochemical properties of amino acids and pairs of amino acids [49]. The database now has a total of 557 indexes, which all come from the published paper. In this research, the physicochemical properties indexes with ‘NA’ are screened out, and the remaining 537 indexes are adopted. Using a physicochemical properties index that represents the amino acids is a good idea to replace the one-hot encoding. However, this method would consume a large amount of storage space and greatly increase the computational complexity. To process the above problems, autoencoder is used to create a more effective representation.

Autoencoder is a neural network that can capture the most important features of data. It consists of the input layer, multiple hidden layers and output layer. The input and output of the autoencoder are consistent, and the goal is to reconstruct itself using some high-level feature recombination. If the number of nodes of the hidden layer is smaller than the input layer and the output layer, it represents the same low-density information and is a centralized representation of the input data obtained from the learning. In our work, we take 537 physicochemical properties indexes of each amino acid as input, reconstruct this informatin with the autoencoder, and then obtain a more effective vector representation from the intermediate layer. The architecture of the autoencoder that we used is shown in Fig. 5.

Fig. 5
figure 5

The structure of autoencoder. It consists of an input layer, five hidden layers and an output layer. The input layer has 537 units; the first two layers and the last two layers of the five hidden layers have 128, 64, 64, and 128 units, respectively, the number of nodes in the middle layer of the hidden layers is an additional hyperparameter obtained by training; and the output layer has 537 units

Results and discussion

The performance of vector representation

To explore the advantage of the new vector representation method, we compared the performance between the one-hot encoding and our method in the sequence information model. We note that our process of parameter adjustment is all on the validation set and not reparameterized to apply on the testing set.

Figure 6 shows the overall success rate of vector representations using different dimensions in Dataset 1 and Dataset 2. To verify whether the vector representation that we proposed is better than one-hot encoding, we also indicated the overall success rate of the one-hot encoding method in the figures. It can be seen from the figure that our method is better in most cases. Taking both the prediction performance and computational complexity into account, we finally chose the vector representation with 10-dimension to replace the one-hot encoding.

Fig. 6
figure 6

Comparison of the overall success rate between our method and one-hot encoding on Dataset 1 and Dataset 2. The x-axis represents the number of dimensions reduced by the autoencoder. The y-axis represents the overall success rate. a and b show the results of the validation set and testing set on Dataset 1; c and d show the results on Dataset 2

Furthermore, to verify whether the vector representation is helpful to improve the performance of the model, we also randomly generate a 10-dimensional vector in the interval [0,1] for comparison experiments; the results are shown in Table 3. From the table, we find that the vector representation that we proposed is more effective than other methods both on Dataset 1 and Dataset 2. In addition, we find that the random method is even worse than the one-hot encoding. Although the random method has a lower dimension, it cannot express the relationship between amino acids well, which makes the result worse. The final experimental results illustrate the reliability of the new vector representation method that we proposed.

Table 3 Comparison of the overall success rate between different vector representation methods on Dataset 1 and Dataset 2

The effective of sequence information model

Table 4 gives the overall success rate of the comparison between the sequence information model and traditional models on Dataset 1 and Dataset 2, where four popular features used for comparative experiments are AAC, DipC, TipC and PseAAC. Among them, AAC, DipC and Tipc are amino acid composition-based methods; PseAAC is the method based on amino acid composition and physicochemical properties. These features are all extracted from the membrane protein sequence. The machine learning algorithm used to compare random forests. Research has illustrated its effectiveness in the prediction of membrane protein types [50]. The optimal parameters in random forests are obtained by using the grid search method and 5-fold cross-validation on the training set. From the table, we find that our DNN model has a better overall success rate than all the machine learning methods. Compared with the DipC, which has the best results among the four sequence features, our model improved by 5.43% and 8.30% on Dataset 1 and Dataset 2, respectively.

Table 4 Comparison of the overall success rate between four traditional models and sequence information model on Dataset 1 and Dataset 2

Tables 5 and 6 give the values of Se, Sp, ACC and MCC for eight membrane protein types using four traditional models and our DNN model on Dataset 1 and Dataset 2, respectively. Although our model is not as good as the traditional model in some indicators, our model achieves better results in most cases. Furthermore, we find that it is difficult for the method combining feature engineering and the machine learning algorithm to predict the small types, especially in Type 4 membrane protein samples where the value of Sn and Mcc is zero in most traditional models. However, our model can predict the small type samples well. The reason may be the powerful feature extraction capabilities of the convolutional layer and Bi-LSTM layer that we used. The experimental results prove that our proposed sequence information model has better performance than traditional models.

Table 5 Comparison of the Se, Sp, ACC and MCC between four traditional models and sequence information model on Dataset 1
Table 6 Comparison of the Se, Sp, ACC and MCC between four traditional models and sequence information model on Dataset 2

The effective of evolutionary information model

Then, we use the same experimental method to compare the performance between the evolutionary information model and seven popular PSSM-based tradition methods proposed by other researchers, including EDPSSM, KPSSM, CPSSM, TriPSSM, RPSSM, CoPSSM and PsePSSM. They are all powerful features extracted from PSSM. Table 7 gives the overall success rate between seven traditional models and the evolutionary information model on Dataset1 and Dataset 2. From the table, we find the best method for the traditional model is TriPSSM, which utilizes PSSM linear probabilities to compute features. Although TriPSSM can obtain the outstanding grades on membrane protein types prediction, we find that it has 8000-dimensional features; the high dimensionality of features in turn increases the computing complexity. Our model is still effective since it can use GPU acceleration to yield better results in less time.

Table 7 Comparison of the overall success rate between seven traditional models and evolutionary information model on Dataset 1 and Dataset 2, respectively

Tables 8 and 9 report the results, measured in Se, Sp, ACC and MCC over all methods on Dataset 1 and Dataset 2. From the tables, we can draw a conclusion that our model is more effective than the traditional model in evolution information extraction. The convolutional layer has a powerful ability to extract features from the evolutionary information, and the CapsNet layer can extract local position information. These abilities explain why our DNN model could achieve better results.

Table 8 Comparison of the Se, Sp, ACC and MCC between seven traditional models and the evolutionary information model on Dataset 1
Table 9 Comparison of the Se, Sp, ACC and MCC between seven traditional models and the evolutionary information model on Dataset 2

The effective of stacking ensemble method

To investigate on how the fusion model impacts the performance, we report the overall success rate of the validation set and testing set using two single models and the fusion model (Table 10). From the table, we find that the fusion model achieves better results both on datasets, demonstrating its superiority in prediction.

Table 10 Comparison of the overall success rate between two single information model and fusion model on Dataset 1 and Dataset 2

Figures 7 and 8 give the comparison results, measured in Se, Sp, ACC and MCC over three models on Dataset 1 and Dataset 2,respectively. We find that the fusion model achieves the best results in most cases. However, the performance of small-scale forecasts has not yet met our expectations. Other effective strategies for solving imbalance problems may contribute to improving our method.

Fig. 7
figure 7

Comparison of two single models with the fusion model on Dataset 1: a Sensitivity (Se); b Specificity (Sp); c Accuracy (ACC) and d Mathews correlation coeffcient (MCC)

Fig. 8
figure 8

Comparison of two single models with the fusion model on Dataset 2: a Sensitivity (Se); b Specificity (Sp); c Accuracy (ACC) and d Mathews correlation coeffcient (MCC)

Furthermore, the prediction of membrane protein types requires not only high classification accuracy but also certain operational efficiency. Table 11 lists the running time for different models when predicting 100 membrane protein samples. We find that the process of fusion takes approximately 1.93e-10s, and the running time of the fusion model is almost the sum of two single models. Though the time complexity of the fusion model is higher, it is also within the tolerance range.

Table 11 Comparison of the running time for different models when predicting 100 membrane protein samples

Comparison with the existing methods

To demonstrate the performance of our DNN model in practical use, we compared our DNN model with eight state-of-the-art machine learning methods in the prediction of membrane protein types. These models all proceed through complex feature engineering, including feature extraction, dimension reduction and so on. Table 12 represents the overall success rate between our DNN model with eight models on Dataset 1 and Dataset 2. Taking the accuracy of prediction and the degree of domain experts dependence into account, we find the DNN model not only achieves the best performance both on Dataset 1 and Dataset 2, but also removes the burden of feature engineering, which illustrates the effectiveness of our model.

Table 12 Comparison of the overall success rate with eight state-of-the-art methods

Conclusion

In this paper, we propose two DNN models to process the sequence information and evolutionary information separately. Among these models, 1D convolutional layers and Bi-LSTM layers are used for the sequence information model, while the evolutionary information model uses 2D convolutional layers and CapsNet layers. In the comparative experiments, we compare these models with the traditional model. Numerous experimental results demonstrate that our two proposed DNN models not only obtain more competitive results in membrane protein type prediction but also can remove the burden of feature engineering and the reliance on domain experts. In addition, to obtain a better prediction performance, the stacking ensemble method is used to fuse the above two models.

Furthermore, we provide a new vector representation method based on autoencoder and the physicochemical properties of amino acids to improve the performance of the sequence information model. Compared with the one-hot encoding, the new vector representation method not only represents the relationship between amino acids well but also can effectively improve the prediction performance.

Finally, our method can be applied not only to membrane protein type prediction but also to other fields of bioinformatics [51–63]. However, there is still room for further investigation. For example, the problem of an imbalanced dataset had a negative effect on the accuracy of small-sized types. Some strategies for processing with the imbalanced dataset may improve the performance of our DNN model.

Availability of data and materials

The related source codes and datasets are available at https://github.com/DragonKnightss/MembraneProteinTypePrediction.

Abbreviations

ACC:

accuracy

AD:

autocorrelation descriptor

Bi-LSTM:

bidirectional long short term memory

CBOW:

Continuous bag of words

CoPSSM:

correlation position-specific scoring matrix

DipC:

dipeptide composition

DNN:

deep neural network

EDPSSM:

evolutionary difference position-specific score matrix

FN:

false negative

FP:

false positive

GD:

global descriptor

HHT:

Hilbert-Huang transform

KNN:

k-nearest neighbor

KPSSM:

k-separated-bigrams position-specific score matrix

LAAC:

Local amino acid composition

LDC:

local dipeptide composition

LZC:

Lempel-Ziv complexity

MCC:

Matthews correlation coefficient

MLR:

multiple logistic regression. NBM: naive bayesian model

NLP:

natural language processing

OSR:

overall success rate

PNN:

probabilistic neural network

PseAAC:

pseudo-amino acid composition

PsePSSM:

pseudo position-specific score matrix

PSSM:

position-specific scoring matrix

RF:

random forests

RPSSM:

reduced position-specific score matrix

SD:

sequence-order descriptor

Se:

sensitivity

Sp:

specificity

SVM:

support vector machines

TipC:

tripeptide composition

TN:

true negative

TP:

true positive

TriPSSM:

tri-gram position-specific score matrix

References

  1. Deng SP, Huang DS. SFAPS: An R package for structure/function analysis of protein sequences based on informational spectrum method. Methods. 2014; 69(3):207–12.

    Article  CAS  PubMed  Google Scholar 

  2. Huang DS, Yu HJ. Normalized Feature Vectors: A Novel Alignment-Free Sequence Comparison Method Based on the Numbers of Adjacent Amino Acids. IEEE/ACM Trans Comput Biol Bioinforma. 2013; 10(2):457–67.

    Article  Google Scholar 

  3. Xia JF, Zhao XM, Song J, et al.APIS: accurate prediction of hot spots in protein interfaces by combining protrusion index with solvent accessibility. BMC Bioinformatics. 2010; 11(1):174–0.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  4. Wang B, Chen P, Huang DS, et al.Predicting protein interaction sites from residue spatial sequence profile and evolution rate. Febs Lett. 2006; 580(2):0–384.

    Article  CAS  Google Scholar 

  5. Zhao XM, Cheung YM, Huang DS. A novel approach to extracting features from motif content and protein composition for protein sequence classification. Neural Netw. 2005; 18(8):1019–28.

    Article  PubMed  Google Scholar 

  6. Almén MS, Nordström KJV, Fredriksson R, et al.Mapping the human membrane proteome: a majority of the human membrane proteins can be classified according to function and evolutionary origin. BMC Biol. 2009; 7(1):50.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  7. Krogh A, Larsson B, Heijne GV, et al.Predicting transmembrane protein topology with a hidden markov model: application to complete genomes. J Mol Biol. 2001; 305(3):0–580.

    Article  CAS  Google Scholar 

  8. Overington JP, Al-Lazikani B, Hopkins AL. How many drug targets are there?Nat Rev Drug Discov. 2006; 5(12):993–6.

    Article  CAS  PubMed  Google Scholar 

  9. Chou KC, Shen HB. MemType-2L: A Web server for predicting membrane proteins and their types by incorporating evolution information through Pse-PSSM. Biochem Biophys Res Commun. 2007; 360(2):0–345.

    Article  CAS  Google Scholar 

  10. Wang L, Yuan Z, Chen X, et al.The prediction of membrane protein types with NPE. IEICE Electron Express. 2010; 7(6):397–402.

    Article  Google Scholar 

  11. Ali F, Hayat M. Classification of membrane protein types using Voting Feature Interval in combination with Chou’s Pseudo Amino Acid Composition. J Theor Biol. 2015; 384:78–83.

    Article  CAS  PubMed  Google Scholar 

  12. Chen YK, Li KB. Predicting membrane protein types by incorporating protein topology, domains, signal peptides, and physicochemical properties into the general form of Chou’s pseudo amino acid composition. J Theor Biol. 2013; 318:1–12.

    Article  CAS  PubMed  Google Scholar 

  13. Nanni L, Lumini A. An ensemble of support vector machines for predicting the membrane protein type directly from the amino acid sequence. Amino Acids. 2008; 35(3):573–80.

    Article  CAS  PubMed  Google Scholar 

  14. Wan S, Mak MW, Kung SY. Mem-ADSVM: A two-layer multi-label predictor for identifying multi-functional types of membrane proteins. J Theor Biol. 2016; 398(7):32–42.

    Article  CAS  PubMed  Google Scholar 

  15. Hayat M, Khan A. Mem-PHybrid: Hybrid features-based prediction system for classifying membrane protein types. Anal Biochem. 2012; 424(1):35–44.

    Article  CAS  PubMed  Google Scholar 

  16. Hayat M, Khan A. Predicting membrane protein types by fusing composite protein sequence features into pseudo amino acid composition. J Theor Biol. 2011; 271(1):10–7.

    Article  CAS  PubMed  Google Scholar 

  17. Chou KC. Prediction of protein cellular attributes using pseudo-amino acid composition. Proteins-Struct Funct Bioinforma. 2010; 43(3):246–55.

    Article  Google Scholar 

  18. Han GS, Yu ZG, Anh V. A two-stage SVM method to predict membrane protein types by incorporating amino acid classifications and physicochemical properties into a general form of Chou’s PseAAC. J Theor Biol. 2014; 344:31–9.

    Article  CAS  PubMed  Google Scholar 

  19. Guo L, Wang S. Prediction for Membrane Protein Types Based on Effective Fusion Representation and MIC-GA Feature Selection. IEEE Access. 2018; 6:75669–81.

    Article  Google Scholar 

  20. Wang T, Xia T, Hu XM. Geometry preserving projections algorithm for predicting membrane protein types. J Theor Biol. 2010; 262(2):208–13.

    Article  CAS  PubMed  Google Scholar 

  21. Anishetty S, Pennathur G, Anishetty R. Tripeptide analysis of protein structures. BMC Struct Biol. 2002; 2(1):9.

    Article  PubMed  PubMed Central  Google Scholar 

  22. Wang T, Yang J, Shen HB, et al.Predicting membrane protein types by the LLDA algorithm. Protein Pept Lett. 2008; 15(9):915–21.

    Article  CAS  PubMed  Google Scholar 

  23. Ding S, Li Y, Shi Z, et al.A protein structural classes prediction method based on predicted secondary structure and PSI-BLAST profile. Biochimie. 2014; 97(2):60–5.

    Article  CAS  PubMed  Google Scholar 

  24. Zhang L, Zhao X, Kong L. Predict protein structural class for low-similarity sequences by evolutionary difference information into the general form of Chou’s pseudo amino acid composition. J Theor Biol. 2014; 355:105–10.

    Article  CAS  PubMed  Google Scholar 

  25. Sharma A, Lyons J, Dehzangi A, et al.A feature extraction technique using bi-gram probabilities of position specific scoring matrix for protein fold recognition. J Theor Biol. 2013; 320:41–6.

    Article  CAS  PubMed  Google Scholar 

  26. Harsh S, Gaurav R, et al.Protein Fold Recognition Using Genetic Algorithm Optimized Voting Scheme and Profile Bigram. J Softw. 2016; 11(8):756–67.

    Article  Google Scholar 

  27. Wang S, Yue Y. Protein subnuclear localization based on a new effective representation and intelligent kernel linear discriminant analysis by dichotomous greedy genetic algorithm. Plos ONE. 2018; 13(4):e0195636.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  28. Patel S, Tripathi R, Kumari V, et al.DeepInteract: Deep Neural Network based Protein-Protein Interaction prediction tool. Current Bioinformatics. 2017; 12(6). https://doi.org/10.2174/1574893611666160815150746.

  29. Veltri D, Kamath U, Shehu A. Deep Learning Improves Antimicrobial Peptide Recognition. Bioinformatics. 2018; 34(16):2740–7.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  30. Sun S, Wu Q, Peng Z, et al.Enhanced prediction of RNA solvent accessibility with long short-term memory neural networks and improved sequence profiles. Bioinformatics. 2019; 35(10):1686–91.

    Article  CAS  PubMed  Google Scholar 

  31. Zhang Z, Zhao Y, et al.Deep learning in omics: a survey and guideline. Brief Funct Genom. 2019; 8(1):41–57.

    Article  Google Scholar 

  32. Guohui C, Hanhui M, Jifang Y, et al.DeepCRISPR: optimized CRISPR guide RNA design by deep learning. Genome Biol. 2018; 19(1):88.

    Article  CAS  Google Scholar 

  33. Wei L, Ding Y, Su R, et al.Prediction of human protein subcellular localization using deep learning. J Parallel Distrib Comput. 2017; 117:212–7.

    Article  Google Scholar 

  34. Derevyanko G, Grudinin S, et al.Deep convolutional networks for quality assessment of protein folds. Bioinformatics. 2018; 34(23):4046–53.

    Article  CAS  PubMed  Google Scholar 

  35. Luo F, Wang M, et al.DeepPhos: prediction of protein phosphorylation sites with deep learning. Bioinformatics. 2019; 35(16):2766–73.

    Article  PubMed  PubMed Central  Google Scholar 

  36. Nie L, Deng L, Fan C, et al.Prediction of Protein S-Sulfenylation Sites Using a Deep Belief Network. Curr Bioinforma. 2018; 13(5):461–7.

    Article  CAS  Google Scholar 

  37. Shen Z, Bao W, et al.Recurrent neural network for predicting transcription factor binding sites. Sci Rep. 2018; 8:15270.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  38. Lecun Y, Bengio Y, Hinton G. Deep learning. Nature. 2015; 521(7553):436.

    Article  CAS  PubMed  Google Scholar 

  39. Graves A, Schmidhuber J. Framewise phoneme classification with bidirectional LSTM and other neural network architectures. Neural Netw. 2005; 18(5):602–10.

    Article  PubMed  Google Scholar 

  40. Sabour S, Frosst N, Hinton GE. Dynamic routing between capsules. Advances in Neural Information Processing Systems. 2017;:3859–69.

  41. Wan S, Mak MW, Kung SY. Benchmark data for identifying multi-functional types of membrane proteins. Data Brief. 2016; 8(C):105–7.

    Article  PubMed  PubMed Central  Google Scholar 

  42. Jones DT. Protein secondary structure prediction based on position-specific scoring matrices. J Mol Biol. 1999; 292(2):0–202.

    Article  CAS  Google Scholar 

  43. Bin L, Kai L, De-Shuang H, et al.iEnhancer-EL: Identifying enhancers and their strength with ensemble learning approach. Bioinformatics. 2018; 34(22):3835–42.

    Article  CAS  Google Scholar 

  44. Wang B, San Wong H, Huang DS. Inferring Protein-Protein Interacting Sites Using Residue Conservation and Evolutionary Information. Protein Pept Lett. 2006; 13(10):999–1005.

    Article  CAS  PubMed  Google Scholar 

  45. Xia J, Zhao X, et al.Predicting protein-protein interactions from protein sequences using meta predictor. Amino Acids. 2010; 39(5):1595–9.

    Article  CAS  PubMed  Google Scholar 

  46. Ting KM, Witten IH. Issues in stacked generalization. J Artif Intell Res. 2011; 10(1):271–89.

    Google Scholar 

  47. Powers DMW. Evaluation: from precision, recall and f-factor to roc, informedness, markedness & correlation. J Mach Learn Technol. 2011; 2(1):37–63.

    Google Scholar 

  48. Mikolov T, Chen K, Corrado G, Dean J. Efficient estimation of word representations in vector space. 2013. arXiv:1301.3781.

  49. Kawashima S, Ogata H, Kanehisa M. AAindex: Amino Acid Index Database. Nucleic Acids Res. 1999; 27(1):368–9.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  50. Sankari ES, Manimegalai D. Predicting membrane protein types using various decision tree classifiers based on various modes of general PseAAC for imbalanced datasets. J Theor Biol. 2017; 435:208–17.

    Article  CAS  PubMed  Google Scholar 

  51. Zhu L, Deng SP, You ZH, et al.Identifying Spurious Interactions in the Protein-Protein Interaction Networks Using Local Similarity Preserving Embedding. IEEE/ACM Trans Comput Biol Bioinforma. 2017; 14(2):345–52.

    Article  Google Scholar 

  52. Zou Q, Xing PW, Wei LY, et al.Gene2vec: gene subsequence embedding for prediction of mammalian N-6-methyladenosine sites from mRNA. RNA. 2019; 25(2):205–18.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  53. Bin L, Fan W, De-Shuang H, et al.iRO-3wPseKNC: identify DNA replication origins by three-window-based PseKNC. Bioinformatics. 2018; 34(18):3086–93.

    Article  CAS  Google Scholar 

  54. Liu B, Yang F, Huang DS, et al.iPromoter-2L: a two-layer predictor for identifying promoters and their types by multi-window-based PseKNC. Bioinformatics. 2018; 4(1):33–40.

    Article  CAS  Google Scholar 

  55. Guo WL, Huang DS. An efficient method to transcription factor binding sites imputation via simultaneous completion of multiple matrices with positional consistency. Mol BioSyst. 2017; 13(9):1827–37.

    Article  CAS  PubMed  Google Scholar 

  56. You ZH, Yin Z, Han K, et al.A semi-supervised learning approach to predict synthetic genetic interactions by combining functional and topological properties of functional gene network. BMC Bioinformatics. 2010; 11(1):343–0.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  57. Chen P, Han K, Li X, et al.Predicting Key Long-Range Interaction Sites by B-Factors. Protein Pept Lett. 2008; 15(5):478–83.

    Article  CAS  PubMed  Google Scholar 

  58. Huang DS, Zheng CH. Independent component analysis-based penalized discriminant method for tumor classification using gene expression data. Bioinformatics. 2006; 22(15):1855–62.

    Article  CAS  PubMed  Google Scholar 

  59. Huang DS, Zhang L, Han K, et al.Prediction of Protein-Protein Interactions Based on Protein-Protein Correlation Using Least Squares Regression. Curr Protein Pept Sci. 2014; 15(6):553–60.

    Article  CAS  PubMed  Google Scholar 

  60. You ZH, Lei YK, Gui J, et al.Using manifold embedding for assessing and predicting protein interactions from high-throughput experimental data. Bioinformatics. 2010; 26(21):2744–51.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  61. Shi MG, Xia JF, Li XL, et al.Predicting protein–protein interactions from sequence using correlation coefficient and high-quality interaction dataset. Amino Acids. 2010; 38(3):891–9.

    Article  CAS  PubMed  Google Scholar 

  62. Xia JF, Han K, Huang DS. Sequence-Based Prediction of Protein-Protein Interactions by Means of Rotation Forest and Autocorrelation Descriptor. 17. 2010; 1:137–45.

    Google Scholar 

  63. Cui G, Chen Y, Huang DS, et al.An Algorithm for Finding Functional Modules and Protein Complexes in Protein-Protein Interaction Networks. Journal of Biomedicine & Biotechnology. 2014; 2008(1110-7243):860270.

    Google Scholar 

Download references

Funding

This research is supported by grants from National Natural Science Foundation of China (11661081), Natural Science Foundation of Yunnan Province (2017FA032) and Training Plan for Young and Middle-aged Academic Leaders of Yunnan Province (2018HB031).

Author information

Authors and Affiliations

Authors

Contributions

LG and SW conceived the algorithm, prepared the datasets, carried out experiments. LG wrote the manuscript. SW, ML, ZC provided some suggestions for the manuscript writing and revised it critically for important intellectual content. All authors read and approved the final manuscript.

Corresponding author

Correspondence to Shunfang Wang.

Ethics declarations

Ethics approval and consent to participate

Not applicable.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver(http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Guo, L., Wang, S., Li, M. et al. Accurate classification of membrane protein types based on sequence and evolutionary information using deep learning. BMC Bioinformatics 20 (Suppl 25), 700 (2019). https://doi.org/10.1186/s12859-019-3275-6

Download citation

  • Published:

  • DOI: https://doi.org/10.1186/s12859-019-3275-6

Keywords