Next Article in Journal
The Increase in the Drug Resistance of Acute Myeloid Leukemia THP-1 Cells in High-Density Cell Culture Is Associated with Inflammatory-like Activation and Anti-Apoptotic Bcl-2 Proteins
Next Article in Special Issue
Melatonin Inhibits hIAPP Oligomerization by Preventing β-Sheet and Hydrogen Bond Formation of the Amyloidogenic Region Revealed by Replica-Exchange Molecular Dynamics Simulation
Previous Article in Journal
Distinct Chemokine Receptor Expression Profiles in De Novo DLBCL, Transformed Follicular Lymphoma, Richter’s Trans-Formed DLBCL and Germinal Center B-Cells
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Identify Bitter Peptides by Using Deep Representation Learning Features

1
College of Biomedical Engineering, Sichuan University, Chengdu 610065, China
2
College of Software Engineering, Sichuan University, Chengdu 610065, China
3
West China School of Medicine, Sichuan University, Chengdu 610065, China
4
Key Laboratory of Coarse Cereal Processing, Ministry of Agriculture and Rural Affairs, College of Food and Biological Engineering, Chengdu University, Chengdu 610106, China
*
Author to whom correspondence should be addressed.
Int. J. Mol. Sci. 2022, 23(14), 7877; https://doi.org/10.3390/ijms23147877
Submission received: 31 May 2022 / Revised: 1 July 2022 / Accepted: 14 July 2022 / Published: 17 July 2022
(This article belongs to the Special Issue State-of-the-Art Molecular Biophysics in China)

Abstract

:
A bitter taste often identifies hazardous compounds and it is generally avoided by most animals and humans. Bitterness of hydrolyzed proteins is caused by the presence of bitter peptides. To improve palatability, bitter peptides need to be identified experimentally in a time-consuming and expensive process, before they can be removed or degraded. Here, we report the development of a machine learning prediction method, iBitter-DRLF, which is based on a deep learning pre-trained neural network feature extraction method. It uses three sequence embedding techniques, soft symmetric alignment (SSA), unified representation (UniRep), and bidirectional long short-term memory (BiLSTM). These were initially combined into various machine learning algorithms to build several models. After optimization, the combined features of UniRep and BiLSTM were finally selected, and the model was built in combination with a light gradient boosting machine (LGBM). The results showed that the use of deep representation learning greatly improves the ability of the model to identify bitter peptides, achieving accurate prediction based on peptide sequence data alone. By helping to identify bitter peptides, iBitter-DRLF can help research into improving the palatability of peptide therapeutics and dietary supplements in the future. A webserver is available, too.

1. Introduction

Humans and most animals instinctively dislike bitter substances, as the taste often identifies toxic compounds. However, some beneficial nutrients, such as soy products, endive, and other Asteraceae vegetables, as well as certain therapeutic peptides, are often bitter. Proteins can be enzymatically digested into shorter polypeptides that have certain beneficial biological activities. Studies have shown that hydrolyzed polypeptides have good nutritional properties and can be easily absorbed and utilized. However, hydrolysis often produces peptides with varying degrees of bitterness that can be detected even at very low concentrations [1]. The bitter taste of protein hydrolysates is caused by the presence of peptides containing hydrophobic amino acids. Most of these peptides are typically composed of no more than eight amino acids and few contain more than ten. However, bitter peptides containing up to 39 amino acids have been described [2]. The bitter taste of protein hydrolysates is the result of a variety of factors. Hydrophobic amino acids within the polypeptide tend to become exposed, stimulating the taste buds, causing the bitterness. Generally, the more hydrophobic amino acids are exposed, the stronger the bitter taste. In addition, the length, overall hydrophobicity, sequence, and amino composition of a polypeptide chain also have a significant impact on bitterness [3].
Identifying bitter peptides using conventional laboratory approaches is expensive and time-consuming. With the availability of a large number of peptide sequence databases in the post-genomic era [4], the development of automated computational models to identify novel bitter peptides has far-reaching practical implications. The success of machine learning-based peptide prediction methods led to an increased interest in bioinformatics [5,6,7,8]. A number of computational methods based on quantitative structure-activity relationship modeling (QSAR) and machine learning have been developed to predict the bitterness of polypeptides [9,10,11,12,13,14]. For example, BitterX [15], BitterPredict [16], iBitter-SCM [17], and iBitter-Fuse [18] using traditional sequence features were proposed to identify bitter peptides and showed increasing performance. Now the methods for identifying bitter peptides sequences are focused on feature engineering. The more meaningful the representation features for sequence are, the better the accuracy of the models will be. In 2021, BERT4Bitter [19] was proposed, which used natural language processing (NLP) heuristic signature coding methods to represent peptide sequences as feature descriptors, and it displayed better accuracy. While in 2022, He et al. [20] proposed mutual information-based meta learning (MIMML) to discover the best feature combination for bitter peptides and attained an independent accuracy of 93.8%. Although great progress has been made in this field, there is still much room for improvement in the performance of machine learning-based bitter peptide identification models using sequences only.
Deep learning refers to an algorithm in machine learning that does not require pre-processing or prior characterization of data, allowing the transformation of raw protein sequences into a form that machine learning can effectively utilize [21]. Inspired by deep learning application in the NLP area, sequence-based deep representation learning of proteins and peptides has emerged as an efficient and accurate method for predicting the characteristics of proteins [22,23,24], such as UniRep [25] and TAPE [26]. The models with deep representation learning features have been applied for identifying anticancer peptides [23], protein subcellular location [27], and post-translational modifications of proteins [28], showing significant improvement of models using traditional sequence features.
Here, we used two protein/peptide sequence-based deep representation learning feature to develop a machine learning based model named iBitter-DRLF for bitter peptides identification. It achieved impressive 10-fold cross-validation results, and independent test accuracy compared favorably with existing methods and traditional, non-deep representation learning methods. Furthermore, iBitter-DRLF could be used more generally and showed superior discrimination in detecting bitter peptides. In addition, based on the results of independent tests, iBitter-DRLF (ACC = 0.944, MCC = 0.889, Sp = 0.977, F1 = 0.952, auPRC = 0.984, auROC = 0.977) outperformed the most advanced predictors currently available. Apart from the feature analysis methods, unified manifold approximation and projection dimensionality reduction (UMAP) was also used to explore the effects of different feature analysis approaches and different depth characterization learning features on model performance. The user-friendly webserver is available at https://www.aibiochem.net/servers/iBitter-DRLF/ (accessed on 1 May 2022). A diagram describing the steps taken during the development of iBitter-DRLF is shown in Figure 1.

2. Result and Discussion

2.1. Results of Preliminary Optimization

To explore embedded features that are useful in identifying bitter peptides, we first used three deep representation learning feature extraction methods, soft symmetric alignment (SSA), unified representation (UniRep), and bidirectional long short-term memory (BiLSTM). To generate each of these, we used three distinct machine learning methods, SVM, LGBM, and RF, to develop models and carry out their initial optimization. Table 1 shows the results of 10-fold cross-validation and independent tests for the three models developed based on the above assumptions. The values in the table represent model performance measures after the optimization of model parameters. The best values achieved for individual features are shown in bold and are underlined.
As shown in Table 1, the 10-fold cross-validation results of the UniRep feature vector, developed using SVM, performed the best of all tested feature/model combinations (accuracy (ACC) = 0.865, Matthews correlation coefficient (MCC) = 0.730, sensitivity (Sn) = 0.867, specificity (Sp) = 0.863, F1 = 0.865, area under the PRC curve (auPRC) = 0.937, and area under the ROC curve (auROC) = 0.931). The ACC of this feature vector exceeded other options by 1.17–9.91%. MCC was better by 2.67–26.69%, Sn by 0.46–6.25%, Sp by 1.77–14.46%, F1 by 0.98–9.08%, and auROC by 0.54–7.63%. Regarding its performance in independent tests (ACC = 0.867, MCC = 0.735, Sn = 0.844, Sp = 0.891, F1 = 0.864, auPRC = 0.952, auROC = 0.948), ACC was 1.85% lower, MCC was reduced by 4.22%, Sn was reduced by 7.35%, F1 was reduced by 2.49%, and auPRC by 0.47% compared to the BiLSTM feature vector developed based on SVM. It can be concluded that, for the identification of bitter peptides, UniRep features were superior to BiLSTM features.

2.2. The Effects of Feature Fusion on the Automatic Identification of Bitter Peptides

At this stage of the development work, we evaluated the use of pairwise combinations of features to generate fusion features. Features were combined in all possible pairings, namely, SSA + UniRep, SSA + BiLSTM, and UniRep + BiLSTM. In addition, SSA, UniRep, and BiLSTM were combined to obtain a triple fusion feature, SSA + UniRep + BiLSTM. These fusion feature combinations were used as input into SVM, LGBM, and RF algorithms to train predictive models and to optimize model performance (Figure 2). Table 2 shows the 10-fold cross-validation and independent test results for all the developed models. The values in this table represent the results after optimizing model parameters. Again, the best performance metrics values are shown in bold and are underlined.
Comparing Table 1 (see Section 3.1) and Table 2, it is immediately apparent that the optimal performance values of the models using fusion features are better than the best values obtained with non-combined features. As an example, the performance of the combination of the 121-dimensional SSA feature and the 5505-dimensional BiLSTM, SSA + BiLSTM, developed using RF, showed an ACC value of 0.898, while the ACC of the SSA feature alone was 0.820, representing a 9.51% better performance of the fusion feature.

2.3. The Effects of Feature Selection on the Automatic Identification of Bitter Peptides

As described in previous section, fused feature encoding was clearly superior to non-fused feature encoding. The sequence vector used in the training set had 512 dimensions, while the feature vectors based on the combined fusion feature scheme had 2021, 3726, 5505, and 5626 dimensions. This high number of dimensions increases the risk of redundancy and the overfitting of feature information. To resolve this problem, we used the LGBM algorithm for feature selection, while using an incremental feature strategy and a hyperparametric mesh search method. For the latter, we selected the scikit-learn GridSearchCV module to perform the hyperparameter search for each model. The performance metrics of each individual, double, and triple fused feature developed using all three machine learning models (SVM, LGBM, RF) are summarized in Table 3, while a visual representation of the outcomes is shown in Figure 3.
The outcome measures shown in Figure 3 and Table 3 clearly indicate that the selected fusion feature sets performed significantly better than unselected fusion features. It is apparent from these results that the overall performance of the 106D UniRep + BiLSTM feature vector was better than outcomes achieved with any other feature vector. The results of 10-fold cross-validation (ACC = 0.899, MCC = 0.777, Sn = 0.891, Sp = 0.887, F1 = 0.889, auPRC = 0.947, auROC = 0.952) outperformed the ACC of other feature sets by 0.91–5.83%. MCC was better by 1.97–14.26%, Sn by 0.45–5.57%, Sp by 1.37–6.61%, and F1 by 0.79–5.70%. The independent test results for the 106D UniRep + BiLSTM fusion feature (ACC = 0.944, MCC = 0.889, Sn = 0.922, Sp = 0.977, F1 = 0.952, auPRC = 0.984, auROC = 0.977) were also better by 0.64–9.90% (ACC), 1.60–23.64% (MCC), 1.77–7.33% (Sp), 2.52–15.76% (Sn), 1.65–10.78% (F1), 1.23–3.42% (auPRC), and 0.62–3.17% (auROC).
These results clearly show that selecting feature descriptors is an effective way of resolving problems with information redundancy and was beneficial in optimizing the prediction performance of the bitter peptide prediction model.

2.4. The Effect of Machine Learning Model Parameter Optimization on the Automated Identification of Bitter Peptides

It is apparent from Table 3 that the overall performance of the SSA + UniRep_106 feature set was superior to all other options. There were only two isolated exceptions to this statement. The ten-fold cross-validation result developed based on LGBM with an auROC = 0.957 and the Sn of SSA + UniRep + BiLSTM_336 feature developed based on SVM (Sn = 0.953) were marginally better. Although these measures are 0.52% and 3.36% higher than those achieved using the UniRep + BiLSTM_106 feature, we believe that the UniRep + BiLSTM_106 feature developed using LGBM provided the best overall performance. Therefore, the UniRep + BiLSTM feature set was selected for final development, using three different machine learning methods to build models.
We utilized the scikit-learn GridSearchCV module to perform a hyperparameter search on each model, recording the corresponding optimal hyperparameters for each, and comparing them with default parameters. The observed values are shown in Figure 4 and in Supplemental Table S1.
As shown in Figure 4 and Supplemental Figure S1, when the UnRep + BiLSTM feature prediction was run using different algorithms and hyperparameters, the best performance was seen with the RF (Nleaf = 2, n_estimators = 300) model and the LGBM (depth = 3, n_estimators = 75) model. Although in the independent tests the Sn = 0.938 of the RF model was marginally better (Sn was 1.73% higher), in every other respect, including ACC, MCC, Sp, F1, auPRC, and auROC, the LGBM-based model showed clearly superior performance in both independent testing and 10-fold cross-validation.
Based on the analysis above, we selected the first 106D features of UniRep + BiLSTM to build the iBitter-DRLF predictor based on the LGBM model and selected the parameter depth = 3 and n_estimator = 75 values for further use. Although in 10-fold cross-validation the auPRC = 0.955 of LGBM model was marginally better (auPRC was 0.84% higher), in every other respect, the LGBM (depth = 3, n_estimators = 75) model showed clearly superior performance in both independent testing and 10-fold cross-validation. The 10-fold cross-validation results of this final prediction model were ACC = 0.889, MCC = 0.777, Sn = 0.891, Sp = 0.887, F1 = 0.889, auPRC = 0.947, and auROC = 0.952. The corresponding independent test results were ACC = 0.944, MCC = 0.889, Sn = 0.922, Sp = 0.977, F1 = 0.952, auPRC = 0.984, and auROC = 0.977.

2.5. Comparison with Existing Methods

We compared the predictive performance of iBitter-DRLF with existing methods, including iBitter-Fuse [18], MIMML [20], iBitter-SCM [17], and BERT4Bitter [19] to assess the effectiveness and utility of our method against its competitors. Independent test results for iBitter-DRLF and the existing methods are compared in Table 4. These results clearly demonstrate that iBitter-DRLF has significantly better ACC, MCC, Sp, and auROC, than existing methods. ACC outperformed other methods by 0.64–11.85%. MCC was 1.60–29.22% better, Sp was 4.16–15.76% higher, while auROC values were up by 1.35–8.08%. These comparisons show that iBitter-DRLF is more reliable and stable than existing algorithms in predicting the bitterness of peptides.

2.6. Feature Visualization of the Picric Peptide Automatic Recognition Effect

Feature visualizations can communicate key data and features through graphics and colors to enable better insight into complex datasets. UMAP is a consistent manifold approximation and projection for the reduction of dimensionality. This algorithm is also suitable for the visual analysis of peptide characteristics. Feature visualization analysis of the automatic recognition of bitter peptides carried out by the UMAP algorithm preserved the characteristics of the original data well while greatly reducing the dimensions of characteristics. Through UMAP feature visualization, differences in feature representation can be clearly shown. Furthermore, the reasons for performance improvements of the model after feature optimization can be explained. The visualization of dimension-reduced features achieved using UMAP is shown in Figure 5. Apparently as compared to Figure 5 A–C, the first 106 features of the UniRep and BiLSTM fusion function, shown in Figure 5D, can better discriminate bitter peptides from non-bitter ones.

2.7. iBitter-DRLF Webserver

To facilitate the widespread use of our algorithm, we developed an iBitter-DRLF webserver that is freely available online at https://www.aibiochem.net/servers/iBitter-DRLF/ (accessed on 1 May 2022) to other investigators for the prediction of bitter peptides. The webserver is easy to use. Just pasting the peptide sequences into the text box, clicking the run button and waiting for a few minutes, the results will be displayed in the web pages. Please see the webserver interface at the website or in Supplementary Figures S2–S4.

3. Materials and Methods

3.1. Benchmark Dataset

The updated benchmark dataset from iBitter-SCM [17] is utilized here for modeling and to make future comparisons easier. Both peptides constructed as non-bitter using the BIOPEP database [29] and those previously experimentally confirmed to be bitter were included in the datasets used in this study. There are 320 bitter peptides and 320 non-bitter peptides in the BTP640 benchmark dataset. The dataset was randomly split into a training subset known as BTP-CV and an independent subset of test peptides known as BTP-TS in order to prevent overfitting the prediction model. The BTV-CV and BTS-TS groups had peptide ratios of 4:1. As a result, while the BTP-TS dataset contains 64 peptides in each category, the BTP-CV dataset contains 256 bitter peptides and 256 non-bitter peptides. Users can obtain both datasets from https://www.aibiochem.net/servers/iBitter-DRLF/ (accessed on 1 May 2022).

3.2. Feature Extraction

To explore the effects of different features on bitter peptide recognition, we used three deep representation learning feature extraction methods, namely, SSA [30], UniRep [25], and BiLSTM [31]. Models were trained on an alternate dataset for the identification of bitter peptides. Different feature encoding schemes were compared to build more comprehensive predictive models.

3.2.1. Pre-Trained SSA Embedding Model

SSA defines a novel measure of similarity between sequences of arbitrary lengths embedded in vectors. First, a peptide sequence was used as input to a pre-trained model and encoded through a three-tier stacked BiLSTM encoder output. The final embedding matrix of each peptide sequence was obtained through a linear layer, R L × 121 , where L is the length of the peptide. Such a model trained and optimized using SSA is referred to as an SSA-embedded model.
Suppose there are two embedded matrix of R L × 121 named F 1 and F 2 for two different peptide sequences of differing lengths, called L 1 and L 2 .
F 1 = [ x 1 , x 2 , , x L 1 ] ,
where x i is a vector of 121D.
F 2 = [ y 1 , y 2 , , y L 2 ]
where y i is also a vector of 121D.
To calculate the similarity between two amino acid sequences represented as F 1 and F 2 separately, a soft symmetry alignment mechanism was developed, in which the similarity between the two sequences was calculated based on their embedded vector as follows:
s ^ = 1 A i = 1 L 1 j = 1 L 2 a ij x i y j 1
a ij is determined by the following Formulas (4)–(7).
φ ij = exp ( x k y j 1 ) k = 1 L 1 exp ( x k y j 1 )
ω ij = exp ( x i y k 1 ) k = 1 L 2 exp ( x i y k 1 )
a ij = ω ij + φ ij ω ij φ ij
A = i = 1 L 1 j = 1 L 2 a ij
These parameters are backfitted with the parameters of the sequence encoder by a fully differentiated SSA. The trained model converts the peptide sequence into an embedding matrix, R L × 121 , and a 121D SSA feature vector is generated by averaging pooling operations.

3.2.2. Pre-Trained UniRep Embedding Model

The UniRep model was trained on 24 million UniRef50 primary amino acid sequences. The model performs next amino acid prediction by minimizing cross-entropy losses, thus learning how to represent proteins internally in the process. Using the trained model, a single fixed-length vector representation of the input sequence was generated by mLSTM (hidden state). The output vector representation was then trained into the best machine learning model. This characterizes the input sequence and enables supervised learning during different bioinformatics tasks.
First, the sequence with L amino acid residues was embedded into a matrix using a single thermal code, R L × 10 . The matrix was then fed into the mLSTM encoder to obtain a hidden state output of R 1900 × L as an embedding matrix. Finally, by an averaging pooling operation, the UniRep feature vector of 1900D was derived.
The calculation of the mLSTM encoder involves the following Equations (8)–(14).
m t = ( X t W xm ) ( h t 1 W hm )
h ^ t = tan h ( X t W xh + m t W mh )
f t = σ ( X t W xf + m t W mf )
i t = σ ( X t W xi + m t W mi )
o t = σ ( X t W xo + m t W mo )
C t = f t C t 1 + i t h ^ t
h t = o t tan h ( C t )
In these equations, represents element-by-element multiplication, h t 1 represents the previous hidden state, X t is the current input, and m t is the current intermediate multiplication state. h ^ t represents the input before the hidden state, f t is the forgotten gate, i t is the input gate, and o t is the output gate. C t 1 is the previous unit state, C t is the current unit state, and h t is the hide state for output. σ is a sigmoid function, while tan h it is a tangent function.

3.2.3. Pre-Trained BiLSTM Embedding Model

BiLSTM is a combination of a forward LSTM and a backward LSTM that captures bidirectional sequence features better than either LSTM model. LSTM obtains the ability to calculate by forgetting and memorizing information. This propagates information that is useful for subsequent computation moments to pass through while discarding useless information and outputting hidden states at each time point. The forgetting memory and output are controlled by the forgetting gate, memory gate, and output gate. These gates are calculated from the hidden state of the previous moment and the current input.
The calculation of BiLSTM involves the following Formulas (15)–(20).
f t = σ ( W f [ h t 1 , X t ] + b f )
i t = σ ( W i [ h t 1 , X t ] + b i )
C ˜ t = tan h ( W C [ h t 1 , X t ] + b C )
o t = σ ( W o [ h t 1 , X t ] + b o )
C t = f t C t 1 + i t C ˜ t
h t = o t tan h ( C t )
Here X t is the current input, h t 1 represents the previous hidden state, C ˜ t is the current cell state, f t is the forgotten gate, i t is the input gate, o t is the output gate, C t 1 is the previous cell state, and h t is the output hidden state. Again, σ is a sigmoid function, while tan h is a tangent function.

3.2.4. Feature Fusion

To establish the best feature combination, first we combined the SSA eigenvector of 121D with the UniRep eigenvector of 1900D, obtaining the SSA + UniRep fusion feature vector, 2021D. Second, the SSA eigenvector of 121D was combined with the BiLSTM eigenvector of 3605D to obtain the 3726D SSA + BiLSTM fusion eigenvector. Third, the 1900D UnIrep eigenvector was combined with the 3605D BiLSTM eigenvector giving the 5505D UniRep + BiLSTM fusion feature vector. Finally, the 121D SSA eigenvector, the 1900D UniRep eigenvector, and the 3605D BiLSTM eigenvector were combined to obtain the 5626D SSA + UniRep + BiLSTM fusion eigenvector.

3.3. Feature Selection Method

LGBM is a gradient boosting framework that abandons the level-wise decision tree growth strategy used by most gradient boosting tools in favor of a leaf-wise algorithm with depth restrictions. In this project, LGBM was utilized to identify the optimal feature space and sort features based on their importance values. Data and data labels were entered into the LGBM model to fit the model before using the built-in functions of LGBM to obtain the importance value for each feature. Features were ranked from ‘largest’ to ‘smallest’ based on feature importance values and those with an importance value greater than the critical value (the average feature importance values) were selected.

3.4. Machine Learning Methods

This work utilized three widely used high performance machine learning models, SVM [32,33], RF [34,35], and LGBM [23,36].
SVM is a typical machine learning (ML) algorithm for dealing with binary classification problems in bioinformatics. We chose gamma and C in the range: a row vector of 30 elements logarithmically divided from 10−4 to 104. ‘rbf’ is the default kernel.
RF is a bagging-based algorithm that not only randomly selects samples, but also randomly selects features during the node splitting process. We selected the range of n estimators as (25, 550) and the range of Nleaf as (2, 12).
LGBM is a gradient boosting framework that uses tree-based learning algorithms. We selected the range of n estimators as (25, 750) and the range of max_depth as (1, 12).

3.5. Evaluation Metrics and Methods

We utilized the following five widely used measures [37,38] to evaluate the performance of specific models, calculated as follows: (21)–(25).
ACC = TP + TN ( TP + TN + FP + FN )
MCC = TP × TN FP × FN ( TP + FP ) ( TP + FN ) ( TN + FP ) ( TN + FN )
Sn = TP ( TP + FN )
Sp = TN ( TN + FP )
F 1 = 2 × TP ( 2 × TP + FN + FP )
In these equations, TP represents the number of bitter peptides correctly predicted to be bitter, TN is the number of non-bitter peptides correctly predicted as non-bitter. FP represents the number of non-bitter peptides erroneously predicted as bitter, while FN is the number of bitter peptides falsely predicted as non-bitter. Using the auROC, the proposed models could be compared with each other and with previously described models. The precision-recall curve is the line connecting the points of precision and recall. Using the auPRC to represent the area enclosed by the precision-recall curve and the x-axis. The area under the ROC curve was also used to evaluate predicted performance, where AUC values ranging from 0.5 to 1 represent stochastic and perfect models, respectively.
K-fold cross-validation and independent testing are widely used to evaluate machine learning models. K-fold cross-validation divides the raw data into K-folds. One is used for the validation of each data subset, while the remaining K − 1 subset is used as the training set. In the validation set, K models are evaluated separately, and the final values of the model measures are averaged to obtain cross-validated values. During the work presented here, we used the 10-fold (K = 10) cross-validation method. For stand-alone testing, a completely different dataset from the training set was used, where all the samples were new to the trained model.

4. Conclusions

Here we describe the development of a new computational model called iBitter-DRLF that can accurately identify bitter peptides based on sequence data alone. It uses a deep representation learning feature embedding method to predict potential bitter peptides. As a result of extensive testing and optimization of multiple feature extraction approaches, using three distinct algorithms, we identified UniRep + BiLSTM_106 as the best fusion feature set. Additional feature selection, using LGBM classifier input, allowed us to develop a robust model. Results of the 10-fold cross-validation and the analysis of the results obtained through independent testing showed that iBitter-DRLF can effectively predict bitter peptides in protein hydrolysates or amongst artificially synthesized peptide therapeutics. Based on independent test results, iBitter-DRLF significantly outperformed existing predictors. Finally, to facilitate the use of the algorithm by other scientists, we built an iBitter-DRLF webserver. We hope that the use of iBitter-DRLF prediction of bitter peptides can improve adherence with taking nutritional supplements and peptide therapeutics in the future and advance drug development and nutrition research.
This work used deep representation learning [39,40] features to improve the predictive performance of the model. Although the exact physicochemical relevance of these features is unclear, this does not prevent the successful use of this method for computational predictions in peptide and protein sequence analysis.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/ijms23147877/s1.

Author Contributions

Conceptualization and methodology, J.J. and Z.L.; Z.L. conceived the ideas; writing—original draft preparation, validation and visualization, J.J. and Z.L.; writing—review and editing, J.J., X.L., Y.J., L.J. and Z.L. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the National Natural Science Foundation of China (No. 62001090) and Fundamental Research Funds for the Central Universities of Sichuan University (No. YJ2021104).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The dataset for model training and testing is available at http://public.aibiochem.net/peptides/BitterP/ or https://github.com/Shoombuatong2527/Benchmark-datasets (accessed on 30 May 2022).

Conflicts of Interest

The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript; or in the decision to publish the results.

References

  1. Maehashi, K.; Huang, L. Bitter peptides and bitter taste receptors. Cell. Mol. Life Sci. 2009, 66, 1661–1671. [Google Scholar] [CrossRef] [PubMed]
  2. Puying, Z.; Xiaoying, Z.; Rui, Q.; Jing, W.; Xionghai, R. Research Progress on Flavour Protease for Debittering. Food Nutr. China 2021, 27, 29–34. [Google Scholar]
  3. Maehashi, K.; Matsuzaki, M.; Yamamoto, Y.; Udaka, S. Isolation of Peptides from an Enzymatic Hydrolysate of Food Proteins and Characterization of Their Taste Properties. Biosci. Biotechnol. Biochem. 1999, 63, 555–559. [Google Scholar] [CrossRef] [PubMed]
  4. Ayana, W.; Marina, S.; Anat, L.; Niv, M.Y. BitterDB: A database of bitter compounds. Nucleic Acids Res. 2012, 40, D413–D419. [Google Scholar]
  5. Cao, C.; Wang, J.; Kwok, D.; Cui, F.; Zhang, Z.; Zhao, D.; Li, M.J.; Zou, Q. webTWAS: A resource for disease candidate susceptibility genes identified by transcriptome-wide association study. Nucleic Acids Res. 2021, 50, D1123–D1130. [Google Scholar] [CrossRef]
  6. Le, N.Q.K.; Hung, T.N.K.; Do, D.T.; Lam, L.H.T.; Dang, L.H.; Huynh, T.-T. Radiomics-based machine learning model for efficiently classifying transcriptome subtypes in glioblastoma patients from MRI. Comput. Biol. Med. 2021, 132, 104320. [Google Scholar] [CrossRef]
  7. Ramzan, Z.; Hassan, M.A.; Asif, H.M.S.; Farooq, A. A Machine Learning-based Self-risk Assessment Technique for Cervical Cancer. Curr. Bioinform. 2021, 16, 315–332. [Google Scholar] [CrossRef]
  8. Su, Q.; Chen, Q. Application of Machine Learning in Animal Disease Analysis and Prediction. Curr. Bioinform. 2021, 16, 972–982. [Google Scholar] [CrossRef]
  9. Hyun-Ock, K.; Li-Chan, E.C.Y. Quantitative Structure—Activity Relationship Study of Bitter Peptides. J. Agric. Food Chem. 2006, 54, 10102–10111. [Google Scholar]
  10. Wu, J.; Aluko, R.E. Quantitative structure-activity relationship study of bitter di- and tri-peptides including relationship with angiotensin I-converting enzyme inhibitory activity. J. Pept. Sci. 2006, 13, 63–69. [Google Scholar] [CrossRef]
  11. Yin, J.; Diao, Y.; Wen, Z.; Wang, Z.; Li, M. Studying Peptides Biological Activities Based on Multidimensional Descriptors (E) Using Support Vector Regression. Int. J. Pept. Res. Ther. 2010, 16, 111–121. [Google Scholar] [CrossRef]
  12. Tong, J.; Liu, S.; Zhou, P.; Wu, B.; Li, Z. A novel descriptor of amino acids and its application in peptide QSAR. J. Theor. Biol. 2008, 253, 90–97. [Google Scholar] [CrossRef] [PubMed]
  13. Lin, Z.-H.; Long, H.-X.; Bo, Z.; Wang, Y.-Q.; Wu, Y.-Z. New descriptors of amino acids and their application to peptide QSAR study. Peptides 2008, 29, 1798–1805. [Google Scholar] [CrossRef] [PubMed]
  14. Liang, G.; Yang, L.; Kang, L.; Mei, H.; Li, Z. Using multidimensional patterns of amino acid attributes for QSAR analysis of peptides. Amino Acids 2008, 37, 583–591. [Google Scholar] [CrossRef]
  15. Huang, W.; Shen, Q.; Su, X.; Ji, M.; Liu, X.; Chen, Y.; Lu, S.; Zhuang, H.; Zhang, J. BitterX: A tool for understanding bitter taste in humans. Sci. Rep. 2016, 6, 23450. [Google Scholar] [CrossRef] [Green Version]
  16. Dagan-Wiener, A.; Nissim, I.; Ben Abu, N.; Borgonovo, G.; Bassoli, A.; Niv, M.Y. Bitter or not? BitterPredict, a tool for predicting taste from chemical structure. Sci. Rep. 2017, 7, 12074. [Google Scholar] [CrossRef] [Green Version]
  17. Charoenkwan, P.; Yana, J.; Schaduangrat, N.; Nantasenamat, C.; Hasan, M.M.; Shoombuatong, W. iBitter-SCM: Identification and characterization of bitter peptides using a scoring card method with propensity scores of dipeptides. Genomics 2020, 112, 2813–2822. [Google Scholar] [CrossRef]
  18. Charoenkwan, P.; Nantasenamat, C.; Hasan, M.M.; Moni, M.A.; Lio, P.; Shoombuatong, W. iBitter-Fuse: A Novel Sequence-Based Bitter Peptide Predictor by Fusing Multi-View Features. Int. J. Mol. Sci. 2021, 22, 8958. [Google Scholar] [CrossRef]
  19. Charoenkwan, P.; Nantasenamat, C.; Hasan, M.; Manavalan, B.; Shoombuatong, W. BERT4Bitter: A bidirectional encoder representations from transformers (BERT)-based model for improving the prediction of bitter peptides. Bioinformatics 2021, 37, 2556–2562. [Google Scholar] [CrossRef]
  20. He, W.; Jiang, Y.; Jin, J.; Li, Z.; Zhao, J.; Manavalan, B.; Su, R.; Gao, X.; Wei, L. Accelerating bioactive peptide discovery via mutual information-based meta-learning. Brief. Bioinform. 2022, 23, bbab499. [Google Scholar] [CrossRef]
  21. Chen, X.; Li, C.; Bernards, M.T.; Shi, Y.; Shao, Q.; He, Y. Sequence-based peptide identification, generation, and property prediction with deep learning: A review. Mol. Syst. Des. Eng. 2021, 6, 406–428. [Google Scholar] [CrossRef]
  22. Arif, M.; Kabir, M.; Ahmad, S.; Khan, A.; Ge, F.; Khelifi, A.; Yu, D.-J. DeepCPPred: A deep learning framework for the discrimination of cell-penetrating peptides and their uptake efficiencies. IEEE/ACM Trans. Comput. Biol. Bioinform. 2021, 1. [Google Scholar] [CrossRef] [PubMed]
  23. Lv, Z.; Cui, F.; Zou, Q.; Zhang, L.; Xu, L. Anticancer peptides prediction with deep representation learning features. Brief. Bioinform. 2021, 22, bbab008. [Google Scholar] [CrossRef] [PubMed]
  24. Sun, Y.-Y.; Lin, T.-T.; Cheng, W.-C.; Lu, I.-H.; Lin, C.-Y.; Chen, S.-H. Peptide-Based Drug Predictions for Cancer Therapy Using Deep Learning. Pharmaceuticals 2022, 15, 422. [Google Scholar] [CrossRef]
  25. Alley, E.C.; Khimulya, G.; Biswas, S.; AlQuraishi, M.; Church, G.M. Unified rational protein engineering with sequence-based deep representation learning. Nat. Methods 2019, 16, 1315–1322. [Google Scholar] [CrossRef]
  26. Rao, R.; Bhattacharya, N.; Thomas, N.; Duan, Y.; Chen, X.; Canny, J.; Abbeel, P.; Song, Y.S. Evaluating Protein Transfer Learning with TAPE. Adv. Neural Inf. Process. Syst. 2019, 32, 9689–9701. [Google Scholar]
  27. Lv, Z.; Wang, P.; Zou, Q.; Jiang, Q. Identification of sub-Golgi protein localization by use of deep representation learning features. Bioinformatics 2020, 36, 5600–5609. [Google Scholar] [CrossRef]
  28. Zhao, Q.; Ma, J.; Wang, Y.; Xie, F.; Lv, Z.; Xu, Y.; Shi, H.; Han, K. Mul-SNO: A Novel Prediction Tool for S-Nitrosylation Sites Based on Deep Learning Methods. IEEE J. Biomed. Health Inform. 2021, 26, 2379–2387. [Google Scholar] [CrossRef]
  29. Minkiewicz, P.; Dziuba, J.; Iwaniak, A.; Dziuba, M.; Darewicz, M. BIOPEP Database and Other Programs for Processing Bioactive Peptide Sequences. J. AOAC Int. 2019, 91, 965–980. [Google Scholar] [CrossRef] [Green Version]
  30. Bepler, T.; Berger, B. Learning protein sequence embeddings using information from structure. arXiv 2019, arXiv:1902.08661v2. [Google Scholar]
  31. Graves, A.; Schmidhuber, J. Framewise phoneme classification with bidirectional LSTM and other neural network architectures. Neural Netw. 2005, 18, 602–610. [Google Scholar] [CrossRef] [PubMed]
  32. Lv, Z.; Wang, D.; Ding, H.; Zhong, B.; Xu, L. Escherichia coli DNA N-4-Methycytosine Site Prediction Accuracy Improved by Light Gradient Boosting Machine Feature Selection Technology. IEEE Access 2020, 8, 14851–14859. [Google Scholar] [CrossRef]
  33. Lv, Z.; Ding, H.; Wang, L.; Zou, Q. A Convolutional Neural Network Using Dinucleotide One-hot Encoder for identifying DNA N6-Methyladenine Sites in the Rice Genome. Neurocomputing 2021, 422, 214–221. [Google Scholar] [CrossRef]
  34. Lv, Z.; Jin, S.; Ding, H.; Zou, Q. A Random Forest Sub-Golgi Protein Classifier Optimized via Dipeptide and Amino Acid Composition Features. Front. Bioeng. Biotechnol. 2019, 7, 215. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  35. Malebary, S.; Rahman, S.; Barukab, O.; Ash’ari, R.; Khan, S.A. iAcety-SmRF: Identification of Acetylation Protein by Using Statistical Moments and Random Forest. Membranes 2022, 12, 265. [Google Scholar] [CrossRef]
  36. Bao, W.; Cui, Q.; Chen, B.; Yang, B. Phage_UniR_LGBM: Phage Virion Proteins Classification with UniRep Features and LightGBM Model. Comput. Math. Methods Med. 2022, 2022, 9470683. [Google Scholar] [CrossRef]
  37. Jiao, Y.; Du, P. Performance measures in evaluating machine learning based bioinformatics predictors for classifications. Quant. Biol. 2016, 4, 320–330. [Google Scholar] [CrossRef] [Green Version]
  38. Wei, L.; Xing, P.; Zeng, J.; Chen, J.; Su, R.; Guo, F. Improved prediction of protein-protein interactions using novel negative samples, features, and an ensemble classifier. Artif. Intell. Med. 2017, 83, 67–74. [Google Scholar] [CrossRef]
  39. Cui, F.; Zhang, Z.; Zou, Q. Sequence representation approaches for sequence-based protein prediction tasks that use deep learning. Brief. Funct. Genom. 2021, 20, 61–73. [Google Scholar] [CrossRef]
  40. Naseer, S.; Hussain, W.; Khan, Y.D.; Rasool, N. NPalmitoylDeep-pseaac: A predictor of N-Palmitoylation Sites in Proteins Using Deep Representations of Proteins and PseAAC via Modified 5-Steps Rule. Curr. Bioinform. 2021, 16, 294–305. [Google Scholar] [CrossRef]
Figure 1. Overview of model development. The pre-trained SSA sequence embedding model, UniRep sequence embedding model, and BiLSTM sequence embedding model were used to embed peptide sequences into eigenvectors. Peptide sequences were converted into 121-dimensional (D) SSA eigenvectors, 1900-dimensional UniRep eigenvectors, and 3605-dimensional (D) BiLSTM eigenvectors. Features were combined and fused to derive the following fusion features: SSA-UniRep (2021D), SSA-BiLSTM (3726D), UniRep-BiLSTM (5505D), and SSA-UniRep-BiLSTM (5626D). These fusion features were used as inputs to the SVM, RF, and LGBM predictor algorithms. The model was optimized by feature selection using the LGBM method. Selected feature sets were subjected to another round of analysis using three algorithms and various hyperparameters. Through 10-fold cross-validation and comparison of independent tests results, the optimized final model was developed. Here, the example like SSA-BiLSTM means two kind of features are combined.
Figure 1. Overview of model development. The pre-trained SSA sequence embedding model, UniRep sequence embedding model, and BiLSTM sequence embedding model were used to embed peptide sequences into eigenvectors. Peptide sequences were converted into 121-dimensional (D) SSA eigenvectors, 1900-dimensional UniRep eigenvectors, and 3605-dimensional (D) BiLSTM eigenvectors. Features were combined and fused to derive the following fusion features: SSA-UniRep (2021D), SSA-BiLSTM (3726D), UniRep-BiLSTM (5505D), and SSA-UniRep-BiLSTM (5626D). These fusion features were used as inputs to the SVM, RF, and LGBM predictor algorithms. The model was optimized by feature selection using the LGBM method. Selected feature sets were subjected to another round of analysis using three algorithms and various hyperparameters. Through 10-fold cross-validation and comparison of independent tests results, the optimized final model was developed. Here, the example like SSA-BiLSTM means two kind of features are combined.
Ijms 23 07877 g001
Figure 2. The accuracy of independent tests for predicting individual and fused features of bitter peptides. Three individual features, SSA, UniRep, and BiLSTM, three double fusion features, SSA + UniRep, SSA + BiLSTM, and UniRep + BiLSTM, and a triple fused feature, SSA + UniRep + BiLSTM, were tested with three distinct machine learning algorithms. The same color is used to signify the same features. The accuracy of the individual or fused feature/machine learning algorithm combinations is ranked from highest to lowest. The accuracy of one of the triple fusion features and four fused double features outperformed the best performing individual feature. In contrast, the four least accurate predictors were individual feature models, indicating the superiority of fused features. The accuracy of the SSA + UniRep + BiLSTM-LGBM) combination was 0.898. The best performing SSA-SVM and BiLSTM-SVM individual features had an accuracy of 0.883. In summary, combining different feature information sets helped to improve the predictive performance of the model. Please note SSA-SVM means that the SVM model with SSA feature vectors as input while and SSA + UniRep + BiLSTM-LGBM means that the LGMB model with the combination of SSA, UniRep and BiLSTM features as input. Other similar labels have similar meanings.
Figure 2. The accuracy of independent tests for predicting individual and fused features of bitter peptides. Three individual features, SSA, UniRep, and BiLSTM, three double fusion features, SSA + UniRep, SSA + BiLSTM, and UniRep + BiLSTM, and a triple fused feature, SSA + UniRep + BiLSTM, were tested with three distinct machine learning algorithms. The same color is used to signify the same features. The accuracy of the individual or fused feature/machine learning algorithm combinations is ranked from highest to lowest. The accuracy of one of the triple fusion features and four fused double features outperformed the best performing individual feature. In contrast, the four least accurate predictors were individual feature models, indicating the superiority of fused features. The accuracy of the SSA + UniRep + BiLSTM-LGBM) combination was 0.898. The best performing SSA-SVM and BiLSTM-SVM individual features had an accuracy of 0.883. In summary, combining different feature information sets helped to improve the predictive performance of the model. Please note SSA-SVM means that the SVM model with SSA feature vectors as input while and SSA + UniRep + BiLSTM-LGBM means that the LGMB model with the combination of SSA, UniRep and BiLSTM features as input. Other similar labels have similar meanings.
Ijms 23 07877 g002
Figure 3. The performance metrics of fusion features using a range of selected features and different algorithms. Panels (A,C,E) show 10-fold cross-validation results, and panels (B,D,F) are independent test results. Please note that the same colors represent different selected features in the different panels. Refer to the codes directly underneath a panel for the appropriately selected feature/model combination codes.
Figure 3. The performance metrics of fusion features using a range of selected features and different algorithms. Panels (A,C,E) show 10-fold cross-validation results, and panels (B,D,F) are independent test results. Please note that the same colors represent different selected features in the different panels. Refer to the codes directly underneath a panel for the appropriately selected feature/model combination codes.
Ijms 23 07877 g003
Figure 4. Performance metrics of the UniRep + BiLSTM features analyzed by different models using default parameters (light bars) or hyperparameters (dark bars). Results using selected hyperparameters invariably matched or outperformed default ones.
Figure 4. Performance metrics of the UniRep + BiLSTM features analyzed by different models using default parameters (light bars) or hyperparameters (dark bars). Results using selected hyperparameters invariably matched or outperformed default ones.
Ijms 23 07877 g004
Figure 5. UMAP was used to visualize the dimension-reduced features of fused features: (A) is the UniRep feature, (B) is the BiLSTM feature, (C) is the BiLSTM + UniRep fusion, and (D) represents data using the first 106 features selected from the BiLSTM_UniRep fusion feature set.
Figure 5. UMAP was used to visualize the dimension-reduced features of fused features: (A) is the UniRep feature, (B) is the BiLSTM feature, (C) is the BiLSTM + UniRep fusion, and (D) represents data using the first 106 features selected from the BiLSTM_UniRep fusion feature set.
Ijms 23 07877 g005
Table 1. Performance metrics of three different deep representation learning features using three machine learning models.
Table 1. Performance metrics of three different deep representation learning features using three machine learning models.
FeatureModelDim10-Fold Cross-ValidationIndependent Test
ACCMCCSnSpF1auPRCauROCACCMCCSnSpF1auPRCauROC
SSA bSVM c1210.8260.6520.8360.8160.8280.890.8980.883 a0.7660.8910.8750.8840.9510.944
LGBM c0.7870.5750.8160.7580.7930.8740.8860.8590.7220.9060.8120.8660.9490.941
RF c0.7910.5840.8280.7540.7980.8480.8650.820.6440.8750.7660.830.9340.922
UniRep bSVM c19000.8650.730.8670.8630.8650.9370.9310.8670.7350.8440.8910.8640.9520.948
LGBM c0.840.680.8280.8520.8380.9390.930.8670.7350.8440.8910.8640.9530.952
RF c0.8420.6840.8360.8480.8410.9270.920.8440.6880.8280.8590.8410.9460.943
BiLSTM bSVM c36050.8180.6370.820.8160.8190.910.9120.8830.7660.9060.8590.8850.9560.951
LGBM c0.8550.7110.8630.8480.8570.9240.9260.8360.6730.8120.8590.8320.950.95
RF c0.8180.6370.8280.8090.820.90.9080.8440.6880.8440.8440.8440.9540.949
a Best performance values are in bold and are underlined. b SSA: soft symmetric alignment; UniRep: unified representation; BiLSTM: bidirectional long short-term memory. c SVM: support vector machine; LGBM: light gradient boosting machine; RF: random forest.
Table 2. Performance metrics of fusion features developed using three machine learning models in 10-fold cross-validation and in independent tests.
Table 2. Performance metrics of fusion features developed using three machine learning models in 10-fold cross-validation and in independent tests.
FeatureModelDim10-Fold Cross-ValidationIndependent Test
ACCMCCSnSpF1auPRCauROCACCMCCSnSpF1auPRCauROC
SSA
+ UniRep b
SVM c20210.8610.7230.875 a0.8480.8630.9290.9270.8670.7340.8590.8750.8660.9540.952
LGBM c0.8400.6800.8480.8320.8410.9330.9240.8590.7190.8590.8590.8590.9600.958
RF c0.8380.6760.8400.8360.8380.9230.9170.8670.7350.8440.8910.8640.9550.954
SSA
+ BiLSTM b
SVM c37260.8360.6720.8480.8240.8380.9150.9170.8830.7660.8590.9060.8800.9430.947
LGBM c0.8480.6960.8590.8360.8490.9270.9270.8750.7510.9060.8440.8790.9610.957
RF c0.8240.6490.8320.8160.8260.9060.9110.8980.7970.8910.9060.8980.9590.951
UniRep
+ BiLSTM b
SVM c55050.8440.6880.8590.8280.8460.9210.9260.8910.7830.9220.8590.8940.9660.962
LGBM c0.8630.7270.8710.8550.8640.9320.9350.8700.7370.8590.8860.8870.9720.958
RF c0.8320.6640.8440.8200.8340.9320.9300.8750.7500.8590.8910.8730.9630.960
SSA
+ UniRep
+ BiLSTM b
SVM c56260.8710.7420.8630.8790.8700.9430.9410.8910.7830.9220.8590.8940.9400.943
LGBM c0.8550.7110.8440.8670.8540.9450.9420.8980.7970.8910.9060.8980.9710.971
RF c0.8400.6800.8480.8320.8410.9260.9250.8980.7990.8590.9370.8940.9630.957
a Values representing the best performance values are in bold and are underlined. b SSA + UniRep: SSA features combined with UniRep features; SSA + BiLSTM: SSA features combined with BiLSTM features; UniRep + BiLSTM: UniRep features combined with BiLSTM features; SSA + UniRep + BiLSTM: all the above features are combined. c SVM: support vector machine; LGBM: light gradient boosting machine; RF: random forest.
Table 3. Performance metrics of individual features and fused features, according to the machine learning methods used to derive them.
Table 3. Performance metrics of individual features and fused features, according to the machine learning methods used to derive them.
FeatureModelDim10-Fold Cross-ValidationIndependent Test
ACCMCCSnSpF1auPRCauROCACCMCCSnSpF1auPRCauROC
SSA bSVM c530.8200.6410.8400.8010.8240.9100.9090.9140.8290.9370.8910.9160.9480.941
LGBM c770.8160.6340.8480.7850.8220.8770.8920.8830.7680.9220.8440.8870.9470.940
RF c160.8050.6100.8200.7890.8080.8600.8810.8670.7340.8750.8590.8680.8880.894
UniRep bSVM c650.8750.7500.8750.8750.8750.9460.9430.9060.8130.8910.9220.9050.9520.952
LGBM c3130.8540.7070.8550.8520.8540.9460.9380.9140.8290.8910.9370.9120.9540.948
RF c3290.8360.6720.8240.8480.8340.9180.9080.8910.7850.8440.9370.8850.9580.957
BiLSTM bSVM c3440.8200.6410.8240.8160.8210.9130.9150.9220.8440.9370.9060.9230.9550.956
LGBM c3390.8710.7420.8830.8590.8730.9250.9290.9060.8130.9060.9060.9060.9690.966
RF c4340.8300.6600.8360.8240.8310.9060.9140.8980.7970.9060.8910.8990.9570.950
SSA
+ UniRep b
SVM c620.8650.7300.8630.8670.8650.9440.9420.9140.8280.9060.9220.9130.9580.957
LGBM c1060.8810.7620.8870.8750.8820.9610.9570.8910.7830.8590.9220.8870.9520.947
RF c470.8380.6760.8590.8160.8410.9370.9310.9060.8160.8590.9530.9020.9560.947
SSA
+ BiLSTM b
SVM c2670.8360.6720.8360.8360.8360.9100.9110.9140.8280.9060.9220.9130.9560.952
LGBM c3170.8610.7230.8750.8480.8630.9240.9290.9060.8130.9060.9060.9060.9620.958
RF c1760.8320.6640.8480.8160.8350.9220.9250.9060.8130.9060.9060.9060.9590.952
UniRep
+ BiLSTM b
SVM c1860.8730.7460.8870.8590.8750.9320.9340.9140.8290.9370.8910.9160.9610.965
LGBM c1060.889a0.7770.8910.8870.8890.9470.9520.9440.8890.9220.9770.9520.9840.977
RF c450.8710.7420.8710.8710.8710.9370.9410.9380.8750.9380.9380.9380.9760.971
SSA
+ UniRep
+ BiLSTM b
SVM c3360.8810.7620.8830.8790.8810.9400.9420.9220.8450.9530.8910.9240.9420.946
LGBM c2850.8810.7620.8910.8710.8820.9510.9470.9380.8750.9220.9530.9370.9690.969
RF c1920.8630.7270.8590.8670.8630.9320.9320.9220.8440.9060.9370.9210.9700.967
a Best performance values are in bold and are underlined. b SSA: soft symmetric alignment; UniRep: unified representation; BiLSTM: bidirectional long short-term memory. SSA + UniRep: SSA features combined with UniRep features; SSA + BiLSTM: SSA features combined with BiLSTM features; UniRep + BiLSTM: UniRep features combined with BiLSTM features; SSA + UniRep + BiLSTM: all the above features are combined. c SVM: support vector machine; LGBM: light gradient boosting machine; RF: random forest.
Table 4. Comparison of performance metrics of bitterness prediction using iBitter-DRLF and alternative state-of-the-art methods while testing independent samples.
Table 4. Comparison of performance metrics of bitterness prediction using iBitter-DRLF and alternative state-of-the-art methods while testing independent samples.
ClassifierACCMCCSnSpauROC
iBitter-DRLF0.944 a0.8890.922 0.9770.977
iBitter-Fuse0.9300.8590.9380.9220.933
BERT4Bitter0.9220.8440.9380.9060.964
iBitter-SCM0.8440.6880.8440.8440.904
MIMML0.9380.8750.9380.9380.955
a Best performance metrics are shown in bold and are underlined.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Jiang, J.; Lin, X.; Jiang, Y.; Jiang, L.; Lv, Z. Identify Bitter Peptides by Using Deep Representation Learning Features. Int. J. Mol. Sci. 2022, 23, 7877. https://doi.org/10.3390/ijms23147877

AMA Style

Jiang J, Lin X, Jiang Y, Jiang L, Lv Z. Identify Bitter Peptides by Using Deep Representation Learning Features. International Journal of Molecular Sciences. 2022; 23(14):7877. https://doi.org/10.3390/ijms23147877

Chicago/Turabian Style

Jiang, Jici, Xinxu Lin, Yueqi Jiang, Liangzhen Jiang, and Zhibin Lv. 2022. "Identify Bitter Peptides by Using Deep Representation Learning Features" International Journal of Molecular Sciences 23, no. 14: 7877. https://doi.org/10.3390/ijms23147877

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop