Next Article in Journal
Selective Membrane Sensor for Aluminum Determination in Food Products, Real Samples and Standard Alloys
Next Article in Special Issue
Rush Hour of LATs towards Their Transport Cycle
Previous Article in Journal
Apneic Tracheostomy in COVID-19 Patients on Veno-Venous Extracorporeal Membrane Oxygenation
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

COMTOP: Protein Residue–Residue Contact Prediction through Mixed Integer Linear Optimization

1
School of Computer Science and Technology, University of Chinese Academy of Sciences, Beijing 100049, China
2
Centre for High Performance Computing, Joint Engineering Research Center for Health Big Data Intelligent Analysis Technology, Shenzhen Institutes of Advanced Technology, Chinese Academy of Sciences, Shenzhen 518055, China
3
Department of Computer Science and Technology, School of Computer Science and Technology, Harbin University of Science and Technology, 52 Xuefu Road, Nangang District, Harbin 150080, China
*
Author to whom correspondence should be addressed.
These authors contribute equally to this work.
Membranes 2021, 11(7), 503; https://doi.org/10.3390/membranes11070503
Submission received: 25 April 2021 / Revised: 24 June 2021 / Accepted: 25 June 2021 / Published: 30 June 2021
(This article belongs to the Collection Membrane Protein Structure and Functions)

Abstract

:
Protein contact prediction helps reconstruct the tertiary structure that greatly determines a protein’s function; therefore, contact prediction from the sequence is an important problem. Recently there has been exciting progress on this problem, but many of the existing methods are still low quality of prediction accuracy. In this paper, we present a new mixed integer linear programming (MILP)-based consensus method: a Consensus scheme based On a Mixed integer linear opTimization method for prOtein contact Prediction (COMTOP). The MILP-based consensus method combines the strengths of seven selected protein contact prediction methods, including CCMpred, EVfold, DeepCov, NNcon, PconsC4, plmDCA, and PSICOV, by optimizing the number of correctly predicted contacts and achieving a better prediction accuracy. The proposed hybrid protein residue–residue contact prediction scheme was tested in four independent test sets. For 239 highly non-redundant proteins, the method showed a prediction accuracy of 59.68%, 70.79%, 78.86%, 89.04%, 94.51%, and 97.35% for top-5L, top-3L, top-2L, top-L, top-L/2, and top-L/5 contacts, respectively. When tested on the CASP13 and CASP14 test sets, the proposed method obtained accuracies of 75.91% and 77.49% for top-L/5 predictions, respectively. COMTOP was further tested on 57 non-redundant α-helical transmembrane proteins and achieved prediction accuracies of 64.34% and 73.91% for top-L/2 and top-L/5 predictions, respectively. For all test datasets, the improvement of COMTOP in accuracy over the seven individual methods increased with the increasing number of predicted contacts. For example, COMTOP performed much better for large number of contact predictions (such as top-5L and top-3L) than for small number of contact predictions such as top-L/2 and top-L/5. The results and analysis demonstrate that COMTOP can significantly improve the performance of the individual methods; therefore, COMTOP is more robust against different types of test sets. COMTOP also showed better/comparable predictions when compared with the state-of-the-art predictors.

1. Introduction

Protein contact prediction aims at predicting which residues of a protein are in contact. Two non-local residues are far away from each other in the protein primary structure, but they are close to each other in the 3D structure. Protein contact prediction is helpful in determining protein structure, model ranking, selection, and evaluation [1,2] and is also important for other fields in evolutionary biology and biotechnology, such as protein function prediction and drug design [3]. A protein contact map is a 2D representation of a protein’s 3D structure. Contact map information can be used as distance restraints to guide protein structure modeling [4,5,6,7,8,9,10]. This creates a new direction for solving the grand challenge of the de novo protein structure. The idea of residue–residue contact prediction and using it to predict 3D models was introduced around two decades ago [11,12]; the realization of that idea has only recently gained much attention by the community and has come into practice as many authors have shown how residue contacts can be predicted with reasonable accuracy [13,14,15,16,17,18,19,20].
The last decade has witnessed significant progress in the development of algorithms for protein contact map prediction. The existing residue–residue contact prediction methods can be broadly classified into three categories: (1) coevolution-based, (2) machine-learning-based, and (3) mathematical optimization technique-based. Coevolution information-based analysis predicts contacts by identifying co-evolved residues in a protein, such as plmDCA [13], COLORS [17], EVfold [4], PSICOV [21], CCMpred [15], GREMLIN [22], MetaPSICOV [16], FreeContact [23], and gDCA [24]. These coevolution-based methods use MSA (multiple sequence alignment) as input, which can be generated using methods like PSI-Blast, HHblits, or Jackhmmer. Protein contact map prediction is regarded as a pattern recognition problem in machine-learning-based approaches, and it is solved using machine learning algorithms. Hidden Markov model, support vector machines, and/or artificial neural networks were used in the early development of such methods; falling into this category of methods are FragHMMent [25], SVMcon [26], SVMSEQ [27], and NNcon [28]. In recent years, with the advancement of deep learning techniques, the precision of the predicted contact maps has increased significantly. Methods in this category include PconsC4 [20], DeepCov [19], RaptorX-contact [29], DNCON2 [30], SPOT-contact [31], DeepContact [32], DeepConPred [33], plmConv [34], DeepFold [35], and so on. In the third group of methods, such as recent ones include COMSAT [18], which is based on support vector machines (SVM), and mixed integer linear programming (MILP) for residue–residue contact prediction of TM proteins. Other optimization methods also solved the problem through satisfying categories of constraints in the MILP models to maximize the probability of the sum of residue contacts [36,37,38,39].
By combining multiple methods, prediction accuracy can be improved because non-systematic errors can be removed by correctly combining them. This idea has been broadly used in different areas of computational biology, such as protein contact prediction [16,18,40,41], secondary structure prediction [42,43,44,45], fold recognition and threading [46,47,48,49], protein disorder prediction [50,51,52], prediction of gene coding sequences [53], clustering and functional explanation of gene-expression data [54], computational studies of protein folding [55,56,57,58], etc. One example of consensus methods for protein contact prediction is COMSAT. It is based on support vector machines (SVM) and mixed integer linear programming (MILP) for contact prediction of TM proteins. It is claimed that this combined method overcomes the cross-validation prediction accuracy barrier of 64.5% and a specificity of 99.4% [18]. LOMETS is an online server for protein structure prediction using the consensus method [47]. CONCORD is a consensus method for secondary structure prediction and uses seven selected secondary structure prediction methods. It achieved an average prediction accuracy of about 83.04% with 6-fold cross-validation on a PDBselect25 protein set [44]. Another consensus method for secondary structure prediction is the consensus data mining (CDM) method. It combines two complementary methods having different strengths: GOR and fragment database mining (FDM) methods. It is claimed that this combined method overcomes the cross-validation prediction accuracy barrier of 80% [42]. Kumar and Carugo [50] developed a consensus method for predicting protein conformational disorder with an accuracy of 81.4% based on 12 methods, using least-squares optimization.
Recently exciting progress has been made on this problem, but contact prediction for proteins without many sequence homologs is still of low quality and insufficient for accurate contact-assisted protein folding [59,60]. Jumper [61] recently established AlphaFold2, a novel strategy that used a different deep learning technique than CASP13 AlphaFold to simulate protein 3D structures. However, it had some targets where the prediction was not very good. The protein targets set at the CASP14 conference do not fully represent all the proteins with many unique structural prediction issues. Thus, the algorithm could not be universally applicable to all proteins.
COMTOP is the first consensus method for protein contact prediction using MILP to maximize the probability of the sum of residue contacts, based on the previous work [44]. This method uses seven selected residue–residue contact prediction methods, including CCMpred [15], EVfold [4], DeepCov [19], NNcon [28], PconsC4 [20], plmDCA [13], and PSICOV [21]. COMTOP maximally combines the strengths of seven protein contact prediction methods by optimizing the number of correctly predicted pairs of residues in the training set. A consensus prediction score based on the confidence scores of the seven individual methods is initiated to assess the likelihood of a residue pair being at one of the protein contact states. Our method performed well compared with seven individual methods when tested by 239 proteins, and a prediction accuracy of about 89.04%, 94.51%, and 97.35% for top-L, top-L/2, and top-L/5 predicted contacts, respectively, was obtained. When tested on CASP13, CASP14, and 57 non-redundant TM proteins, the consensus method achieved accuracies of 75.91%, 77.49%, and 73.91% for top-L/5 predictions, which was better than the seven individual methods and could achieve state-of-the-art prediction performance.

2. Materials and Methods

2.1. Data Description

High-quality training sets and test sets are crucial for the development and validation of prediction models. To train our proposed methods, a training set and a validation set were constructed as follows. We first downloaded a list of 3298 protein chains from PISCES website [62] with a maximum sequence identity of 20%, a maximum R-factor of 0.3, and resolutions better than 2.0 Å. Further, we removed the protein chains with less than 50 amino acids, and 3133 chains (Table S1) were left for our database construction. Then we generated the confidence score by the seven locally installed methods (CCMpred, EVfold, DeepCov, NNcon, PconsC4, PlmDCA, PSICOV) based on 3133 proteins, but some methods could not predict the confidence score for some protein chains. We deleted those proteins from our dataset because our method combines seven individual methods. Following this criterion, we had 1189 proteins in our dataset. The selected protein set was divided into two parts: a training set with 950 proteins (Table S2) and a testing set with 239 proteins (Table S3). Then we ranked the confidence score for the training set and test set, and prepared the dataset for MILP model by selecting the top-5L, top-3L, top-2L, top-L, top-L/2, and top-L/5 predictions where L was the length of a protein.
The second test set was based on the CASP13 protein domains (http://www.predictioncenter.org/download_area/CASP13/, accessed on 5 March 2021). There were 32 target domains, and certain methods were unable to predict the confidence score for two target domains, so they were excluded from the test set, leaving 30 target domains in the CASP13 dataset, which are listed in Table S4 of the Supplementary Materials.
The third test dataset was based on the CASP14 target domains (https://www.predictioncenter.org/casp14/index.cgi, accessed on 5 March 2021). There were 45 target domains, and 28 of them had PDB IDs. After excluding the targets that were unable to obtain confidence scores from all the seven methods, the CASP14 dataset contained 25 target domains, as shown in Table S5.
Finally, we assessed our method on a non-redundant α-helical TM protein test set consisting of 57 α-helical TM proteins through culling the α-helical TM proteins in the PDBTM database against the training test sets of COMTOP and the individual methods NNcon, DeepCov, and PconsC4 (EVfold, plmDCA, PSICOV, and CCMpred had no training sets) with a maximum sequence identity of 25%, a maximum R-factor of 0.3, and resolutions better than 2.0 Å. The 57 α-helical TM proteins are listed in Table S6 of the Supplementary Materials.

2.2. Method Description

2.2.1. List of the Sets, Parameters and Variables

This section lists the sets, parameters, and variables used in this method.
(I)
Indices and sets
(I,J): set of pairs of amino acid positions of a protein, iI, jJ;
P: set of training proteins, pP;
M: set of the methods used, [1,2,3,4,5,6,7], m∈M. Seven methods were used in this consensus method: m = 1 indicates the CCMpred method; m = 2 indicates the DeepCov method; m = 3 indicates the EVfold method; m = 4 indicates the NNcon method; m = 5 indicates the PconsC4 method; m = 6 indicates the plmDCA method; m = 7 indicates the PSICOV method; and subsetP(I,J)(p,(I,J)): subset indicates the number of pairs of residues for each protein p.
(II)
Parameters
confS(p,(i,j),m): the confidence score predicted by method m for pairs of residues (i,j) of a protein p, pP, (i,j)∈subsetP(I,J)(p,(I,J)), m∈M. Firstly the score was normalized for each protein; we used min-max normalization,
z i j = x i j min ( x i j max x i j   min ( x i j
where x i j is the ith and jth score for a residue pair (i,j) of a protein. Then we selected the top 4.5 L prediction score from the top 5 L prediction score for each protein and created a score matrix by using top 4.5 L scores of method M; however, the score matrix generated nan-value because all methods were unable to predict the same pair of residues of a protein. Thus, the predictions scoring matrix dealt with missing values by the following strategy;
Missing   value   = x p x p y p     1000 , when   the   pairs   is   indeed   contact z p     , when   the   pairs   is   not   a   contact  
where x p , y p , and z p are the minimum value, maximum value, and average value of the p-th protein. P: set of training proteins, [1…950], pP; predSS(p,(i,j),m): the training model label matrix depends on the real prediction, 1 for a true prediction and 0 for a false prediction, pP, (i,j)∈subset P(I,J)(p,(I,J)), m∈M. So, we handled the missing nan value by this strategy for label matrix;
missing   value   = 1 0 when   the   pairs   are   not   a   contact     when   the   pairs   are   indeed   contact  
(III)
Binary variables
y(p, (i,j)): equals to 1 if the sum of the scores of the correct contact predictions is higher than the sum of the incorrect ones for residue pair (i,j) of a protein p by at least ϵ p , i , j , pP, i,j∈subset(I,J) (P, (I,J));
y2(p): equals to 1 if the sum of the scores of the correct contact predictions of all pairs of amino acid (i,j) of a protein p is higher than the sum of the score of the incorrect prediction of the same protein p by at least ϵ 2 p , pP, i,j∈subset(I,J) (P, (I,J));
(IV)
Positive variables
λ m : the weight variables for different methods, 0 λ m 1 , m∈M;
ϵ p , i , j : a soft margin variable for the binary variable y((p, (i,j)), pP, i,j∈subset(I,J) (P, (I,J)) (see Section 2.2.1 (III)); and
ϵ 2 p : a soft margin variable for the binary variable y2(p), pP (see Section 2.2.1 (III))

2.2.2. The Training Objective Function

For protein contact prediction, the training objective function of the MILP model takes the following format:
M a x p , i , j y ( p , i , j p , i , j ϵ p , i , j ,   p , i , j s u b s e t P I , J P ,   I , J
where y(p,(i,j)) is a set of binary variables, and it equals to 1 if the sum of the scores of the correct contact predictions is higher than the sum of the incorrect ones for pairs of amino acids (i,j) of a protein p by at least ϵ(p,(i,j)), and this objective function is to maximize the total number of pair of residues. ϵ p , i , j is included here to minimize the sum of soft margins.
The training objective function was conducted on the individual contact/pairs of residues of a protein. The principle of this function is that some protein contact prediction approaches have better prediction performance in some native regions of a protein than other approaches. The consensus approach aims to identify the correct contact prediction for proteins from various approaches by relying on confidence scores for each contact of residues in a protein.

2.2.3. The Model Constraints

For the protein contact prediction, there are two basic constraints in the consensus scheme. The first constraint makes sure that the binary variable y(p,(i,j)) is equal to zero for each contact of residues of a protein if the difference between the sum of the scores of correct contact predictions and the sum of the scores of incorrect predictions from different methods is lower than e(p,(i,j)); this constraint is expressed as:
m λ m     c o n f S p , i , j , m     1 p r e d S S p ,   i , j , m m λ m     c o n f S p , i , j , m     p r e d S S p , i , j , m + ϵ p , i , j < 1 y p ,   i , j ,   p ,   i , j s u b s e t P I , J P , I , J ,   m M
The second type of constraint used in the model normalizes the weights terms λ m of the seven methods.
m λ m = 1 , λ m 0 ,   m M

2.2.4. The Prediction Score and Prediction Label

By carefully selecting the weight variables λ m of each individual method through the MILP optimization-based approach, the consensus method was developed to score higher for the correct contact predictions than for the incorrect contact predictions from different methods. The consensus method ensures that pairs of residues have contact if the sum of scores of correct predictions is higher than the sum of scores of incorrect predictions of each pairs of residues for different methods. It is expressed as
S p , i , j ,   m   =   m λ m     c o n f S p , i , j , m   ,             p ,   i , j s u b s e t P I , J P , I , J ,   m M
where S(p,(i,j),m) is the consensus contact confidence score for (i,j) residues of pth protein and mth method, conf(p,(i,j), m) is the confidence score for the mth method, λ(m) is the weighting factor for the mth method.
The consensus label matrix is the prediction result of method m for the (i,j) residues of a protein p (a value of 1 corresponds to a true prediction, a value of 0 corresponds to a false prediction). It is expressed by
Label   matrix = 1 ,   returned   prediction   by   methods   0 ,   not   returned   prediction   by   methods

2.2.5. The Training and Prediction Procedure

The MILP based training system uses CPLEX (ILOG CPLEX 8.0 reference manual) to optimize the training objective function, from which the weight parameters λ m are attained. The training system for each fold takes around two/four weeks. For the prediction procedure, the seven individual programs run in parallel rather than serially. Once we get the results from the individual methods, the prediction from the MILP model of COMTOP can finish in one second. The running time taken by COMTOP depends on the slowest time of the seven methods, and each protein contact prediction takes around 5–25 min, depending on the size of the database used for sequence-profile/MSA analysis.

2.2.6. Evaluation Measures for Prediction Performance

The effectiveness of our proposed method was evaluated by five widely used metrics: the prediction accuracy, coverage, specificity, negative predictive value (NPV), and Matthews’s correlation coefficient (MCC).
The accuracy is defined as the ratio of correct predictions to total predictions. Accuracy can also be written in terms of true positives (TP) and false positives (FP), as shown in Equation (9). A higher value of accuracy means a better contact prediction model.
Accuracy = N corr N pred = TP TP + FP
where N corr is the number of correctly predicted protein contacts, N pred is the number of total predicted contacts, TP is the number of true positive contacts, and FP is the number of false positive contacts.
Coverage, also called the true positive rate or referred to as “sensitivity” is defined as the ratio of correct predictions to the number of protein contacts in the native structure, as shown in Equation (10). A higher value of coverage means a better contact prediction model.
Coverage = N corr N native = TP TP + FN
where, N corr is the number of correctly predicted protein contacts, N native   is the number of protein contacts in the native structure, and FN is the number of false negative contacts.
Specificity, also called the true negative rate, is the percentage of predicted contacts that are present in the native structure, as shown in Equation (11). It denotes how good the test is at identifying negative conditions.
Specificity = TN TN + FP
where TN is the number of true negative contacts.
When a screening test returns a negative result, the negative predictive value (NPV) is the probability of properly detecting all pairs of residues that do not have contact from among pairs of residues that might or might not have contact and can be calculated using the following formula:
NPV = TN TN + FN
The last metric used to measure the performance of the contact prediction method was Matthew’s correlation coefficient (MCC), a measure of the quality of two-class classifications, which can be calculated using the following formula:
MCC = TP × TN FP × FN TP + FP     TP + FN     TN + FP     TN + FN
The proposed method was developed with the purpose of producing higher accuracy, so we consider that accuracy should have a higher weight than the other metrics.

3. Results

3.1. Performance Evaluation Based on the Training Set

The overall workflow of protein contact prediction is illustrated in Figure 1. The measure of COMTOP performance depends on the weight values given by seven individual methods, which balances the accuracy, coverage, specificity, and MCC. For the training process, a set of optimal parameters from the MILP model were generated (listed in Supplementary Table S7); these parameters were the weights for the seven individual methods that should be used in the consensus prediction model.
Figure 2 shows the weight values and overall performance with different numbers of sample sizes for the training dataset. Figure 2A shows that PconsC4 generated the highest weight value and DeepCov generated the second highest weight value, while the PSICOV and NNcon generated the lowest weight values at different sample sizes. The fact that different systems reveal different weights is owing to the different prediction accuracies of each method. Similarity between seven different approaches plays an important role in determining the weight values (Table 1). Jaccard’s similarity coefficient value was 0.585 between PconsC4 and DeepCov, so PconsC4 and DeepCov generated the highest weight values. CCMpred, PSICOV, and EVfold are similar type methods (see Table 1) and generated lowest weight values. This explains why CCMpred, plmDCA, NNcon, and PISICOV methods had the smallest weights. On the other hand, note that even though the CCMpred method has the third highest accuracy after PconsC4 and DeepCov for contact prediction, its weight for two cases was very small. This is because many of these seven contact prediction methods use PSI-Blast, HHblits, or Jackhmmer (e.g., CCMpred, EVfold, plmDCA, and PISICOV), artificial neural networks use NNcon and deep learning approaches (e.g., DeepCov, PconsC4), and the prediction results of the different methods may correlate with each other in some fashion.
Figure 2B shows overall performance based on the training dataset, the average accuracies, coverages, specificities, and MCCs for COMTOP are plotted against different sample sizes. The prediction accuracy, coverages, and specificities were highest when the sample sizes were small, such as top-L/5 and top-L/2 predictions, respectively, while the accuracy and coverage decreased monotonically with the increasing sample size. This represents the classic trade-off phenomenon common to many prediction problems. Although MCC is also an important estimator in protein contact prediction evaluation, this value increased with increasing sample size. The prediction accuracy of the different training models based on the different sample size was 98.99%, 96.88%, 92.18%, 89.91%, 85.77%, and 79.09%, respectively.

3.2. Performance Evaluation Based on the Independent Set

The measure of COMTOP’s performance for residue contact prediction depends on the best range that balances the accuracy, coverage, specificity, NPV, and MCC. We used four datasets to evaluate the performance of our method (see Section 2.1). Figure 3 and Figure 4 and Table 2 summarize the performance of COMTOP in terms of accuracy (alternatively known as positive predictive value) compared with the seven individual methods when applied to the 239 test proteins.
We have evaluated the accuracy of the top L/k (k = 5, 2, 1) and top K*L (K = 5, 3, 2, 1) predicted contacts where L is the length of a protein. The prediction accuracy is defined as the percentage of native contacts among the top L/k and KL predicted contacts. On 239 proteins, we evaluated our method using the suggested weight values. The prediction results of our method, together with the other seven individual methods, are shown in Table 2 and Figure 3. We observe that our method generated the highest prediction scores compared with the seven individual methods.
The COMTOP model has six sub-models (see in Table S7), and the accuracy of the COMTOP model was 59.68%, 70.79%, 78.86%, 89.04%, 94.51%, and 97.35% for top-5L, top-3L, top-2L, top-L, top-L/2, and top-L/5 predictions, respectively. Among the seven methods, the PconsC4 method had the highest score of accuracy with 34.79%, 47.35%, 59.85%, 78.29%, 89.41%, and 95.76% for top-5L, top-3L, top-2L, top-L, top-L/2, and top-L/5 predictions, respectively. The DeepCov method had the second highest score of accuracy with 29.89%, 41.27%, 51.92%, 70.21%, 83.28%, and 91.61% for top-5L, top-3L, top-2L, top-L, top-L/2, and top-L/5 predictions, respectively. On the other hand, the NNcon method had the lowest score of accuracy with 15.66%, 18.27%, 21.24%, 27.44%, 34.34%, and 43.18% for top-5L, top-3L, top-2L, top-L, top-L/2, and top-L/5 predictions, respectively. The prediction accuracies of all the methods were ranked as follows: PconsC4 > DeepCov > CCMpred > EVfold > plmDCA > PSICOV > NNcon, whereas the ranking of all methods by weight values was as follows: PconsC4 > DeepCov > CCMpred > plmDCA > EVfold > NNcon > PSICOV. The weight values are highly dependent on an individual method’s performance. The improvement of COMTOP in accuracy over the seven individual methods increased with the increasing number of predicted contacts. For example, the improvements of COMTOP in accuracy over the best individual method PconsC4 were 24.89%, 23.34%, 19.01%, 10.75%, 5.1%, and 1.59% for top-5L, top-3L, top-2L, top-L, top-L/2, and top-L/5 predictions, respectively.
Figure 3A shows the overall performance of the COMTOP model in terms of its average accuracy. For comparison, results for CCMpred, EVfold, DeepCov, NNcon, PconsC4, PlmDCA, and PSICOV are also shown. Clearly, the COMTOP model significantly outperforms the other individual methods. To see this in more detail, Figure 3B–E shows the Bland–Altman plot indicating the relationship of COMTOP against PconsC4 and DeepCov. For Figure 3B–E, the majority of the points are above the zero line, and about 95%, 95.8%, 95.4%, and 97.5% fall within the confidence limit, respectively. The mean/bias values of the differences are also all positive, indicating that COMTOP outperformed DeepCov and PconsC4 for test datasets [63].
In Figure 4, the average accuracies, coverages, specificities, NPVs, and MCCs for COMTOP are plotted against the different sample sizes on the 239 non-redundant proteins. The prediction accuracies and specificities are highest when the sample sizes are small, such as top-L/5 predictions, while the accuracy and specificities decrease monotonically with the increasing sample size. The accuracies of the COMTOP model for top-5L, top-3L, top-2L, top-L, top-L/2, and top-L/5 predictions were 59.68%, 70.79%, 78.86%, 89.04%, 94.51%, and 97.35%, respectively, while the corresponding NPVs were 99.22%, 95.41%, 91.66%, 82.34%, 70.10%, and 59.85%, respectively [64]. The coverage improves monotonically with the increase of the sample size. This represents the classic trade-off phenomenon common to many prediction problems. Although MCC is also an important estimator in protein contact prediction evaluation, this value is highest when the sample size is top-L and top-2L. More importantly, top-L is the best range for protein contact prediction. Concerning the number of contacts required for accurate folding, the top-L contacts have been shown to produce good results [65,66]; nevertheless, the researchers have recommended that the number of contacts required be specific to the prediction methods.

3.3. Testing on CASP13 Targets

The critical assessment of protein structure prediction (CASP) is a biennial worldwide competition for protein structure prediction, identifying what progress has been made and highlighting where future effort may be most productively focused. The competition unfolds in a double-blind fashion: The structures of the target domains are unknown to the predictors and the organizers (http://predictioncenter.org/download_area/CASP13/, accessed on 25 April 2021).

3.3.1. Comparison of COMTOP’s Performance with the Seven Individual Methods

In this section, we tested COMTOP on 30 domains in CASP13 using the proposed parameter sets listed in Table S7. The prediction results of COMTOP, together with the other seven individual methods, are shown in Table 3 and Figure 5A, and the COMTOP model significantly beat the other individual techniques. As shown in Table 3, the improvement of COMTOP in accuracy over the seven individual methods increases with the increasing number of predicted contacts. For example, the improvement of COMTOP in accuracies over the top individual method DeepCov were 27.1%, 23.40%, 22.0%, 14.1%, 9.5%, and 0.81% for top-5L, top-3L, top-2L, top-L, top-L/2, and top-L/5 predictions, respectively.
Figure 5B–E as Bland–Altman plots shows the relationship of COMTOP against PconsC4 and DeepCov. For Figure 5B–E, the majority of the points are above the zero line, and about 93.4%, 93.4%, 93.4%, and 96.7% fall within the confidence limit, respectively. The mean/bias of the differences are also shown to be all positive, indicating that COMTOP outperformed DeepCov and PconsC4 for CASP13 datasets. In Figure 6, the average accuracies, coverages, specificities, NPVs, and MCCs for COMTOP are plotted against the different number of contact predictions for 30 domains on CASP13 datasets. For top-L/5 predictions, the accuracy and specificity of the COMTOP model were 75.91% and 91.24%, respectively, while NPVs and coverage were 64.66% and 61.45%. Although MCC is also an important estimator in protein contact prediction evaluation, this value was highest when the sample size was top-5L and top-2L, respectively.
Table 4 shows the accuracies achieved by COMTOP on the domains classified as FM/TBM-easy/TBM-hard/FM/TBM based on CASP13 datasets. Over these target domains, COMTOP achieved average accuracies of 75.91% and 73.90% when considering the top-L/5 and top-L/2 predictions. For top-L/5 predictions, COMTOP showed prediction accuracies larger than 90% for 18 domains and accuracies of about 100% for 16 of these domains. Remarkably, our system obtained high accuracies for TBM-easy, TBM-hard, and FM/TBM classification domain. Among all the individual methods, DeepCov performed best with an accuracy value of 21.40%, 29.80%, 38.10%, 52.60%, 64.40%, and 75.10% for top-5L, top-3L, top-2L, top-L, top-L/2, and top-L/5 predictions, respectively. CCMpred, EVfold, and plmDCA methods had the same ranking as the test set, with 239 proteins, but DeepCov, PconsC4, PSICOV, and NNcon showed different rankings. The prediction accuracies of all the methods were ranked as follows: DeepCov > PconsC4 > CCMpred > EVfold > plmDCA > NNcon > PSICOV, whereas the ranks of all methods based on the test set with 239 proteins were as follows: PconsC4 > DeepCov > CCMpred > EVfold > plmDCA > PSICOV > NNcon.
The prediction accuracy of COMTOP was reduced from an accuracy of 97.35% on the 239 proteins to 75.91% on the CASP13 dataset for the top-L/5 predictions. Moreover, the prediction accuracy of COMTOP was reduced about 11.18%, 17.59%, 18.76%, 22.34%, 20.61%, and 21.45% in CASP13 for top-5L, top-3L, top-2L, top-L, top-L/2, and top-L/5 predictions, respectively. Among the seven individual methods, PconsC4 method had the largest decrease in accuracy, about 13.29%, 17.25%, 29.75%, 28.59%, 30.51%, and 31.86% in CASP13 for top-5L, top-3L, top-2L, top-L, top-L/2, and top-L/5 predictions, respectively. The CCMpred method had the second largest decrease in accuracy, about 8.91%, 12.14%, 19.56%, 22.13%, 27.88%, and 31.70% in CASP13 for top-5L, top-3L, top-2L, top-L, top-L/2, and top-L/5 predictions, respectively. On the other hand, PSICOV performed better for CASP13 targets than the test dataset only for top-L/5 predicted contacts; the prediction accuracy improved from 40.45% to 42.06%.

3.3.2. Comparison of COMTOP’s Performance with a Few State-of-the-Art Schemes

Finally, we compared our COMTOP system with a few state-of-the-art systems that were not used for developing COMTOP; another set of five contact prediction systems was chosen: RapterX [29], Yang_Server [67], TripletRes [68], ResTriplet [68], and DNCON3 [69]. For this assessment, 20 domains were selected from the CASP13 dataset for which the native structures were publicly available, and all 7 methods generated the results. These 20 domains included both FM (free-modeling) and TBM (template-based modeling) domains. To get the predicted contacts for other methods, we evaluated the contact predictions of the top-L/2 and top-L/5 groups in CASP13 over these 20 targets from the webserver at https://predictioncenter.org/casp13/rrc_results.cgi, accessed on 25 April 2021.
The performance of the COMTOP scheme, together with five state-of-the-art systems in terms of accuracy of top-L/2 and top-L/5 predictions is shown in Table 5. From the table, we can see that COMTOP performed better than Yang_server, TripletRes, ResTriplet, and DNCON3 schemes for the top-L/2 contacts and better than TripletRes, ResTriplet, and DNCON3 schemes for the top-L/5 predictions. On the other hand, our model had accuracies of 84.02% and 88.87% for top-L/2 and top-L/5 predictions, respectively, after RapterX. In addition, Yang_server has a slightly better performance than our method for top-L/5 predictions.

3.4. Testing on CASP14 Targets

In this section, we tested COMTOP on 25 CASP14 target domains using the proposed parameter sets. These 25 domains in the CASP14 dataset included both FM and TBM domains, for which the native structures were publicly available.

3.4.1. Performance Comparison of COMTOP against the Seven Individual Methods

Table 6 shows the COMTOP prediction results alongside the seven individual methods, and the COMTOP model significantly outperformed the other methods.
For top-L/2 and top-L/5 predictions, COMTOP achieved average accuracies of 68.33% and 77.49%, respectively. Our prediction accuracies were more than 90% for 13 domains, with 11 of these domains showing 100% accuracies for top-L/5 predictions. Among all the individual methods, DeepCov performed best with accuracies of 67.65% and 71.86% for top_L/2 and top_L/5 contacts, respectively. The overall ranking differed from the CASP13 dataset. DeepCov, PconsC4, PSICOV, and plmDCA methods had the same ranking as the CASP14 dataset, but the ranking of EVfold, NNcon, and CCMpred changed. The prediction accuracies of all the methods were ranked as follows: DeepCov > PconsC4 > NNcon > CCMpred > plmDCA > EVfold > PSICOV. As shown in Table 6, the improvement of COMTOP in accuracy over the seven individual methods increased with the increasing number of predicted contacts.
For the CASP14 dataset, COMTOP’s prediction accuracy for top-L/5 contacts increased by 1.58% compared with the CASP13 dataset; overall it reduced by 7.23%, 6.28%, 7%, 5.52%, and 5.57% for the top-5L, top-3L, top-2L, top-L, and top-L/2 predictions. In addition, COMTOP’s prediction accuracy for the top-L/2, and top-L/5 contacts decreased by 26.18% and 19.86% in the CASP14 dataset from the 239 proteins test dataset. Other methods also showed a decreasing trend in accuracy for CASP14 than for CASP13 and the 239 proteins test dataset, indicating that CASP14 was more difficult than others. Among the seven individual methods, the best performing method, DeepCov, decreased in accuracy for top-L/5 contacts, approximately by 3.24% and 19.75% in the CASP14 from the CASP13 and the test set with 239 proteins, respectively. The second-best method, PconsC4, showed a 5.56% improvement in accuracy for top-L/5 contacts in the CASP14 from CASP13 dataset but a 26.3% decrease from the 239 proteins test dataset.

3.4.2. Performance Comparison of COMTOP against State-of-the-Art Methods

We compared our COMTOP model against a group of seven state-of-the-art contact prediction systems that were recently developed: MULTICOM-AI [70], Kiharalab_Contact [71], tFold, DeepPotential, trfold, RapterX [29], and TripletRes [68]. To get the prediction accuracy for state-of-the-art systems, we took the contact predictions of the top-L/2 and top-L/5 predictions of these methods from the webserver at https://predictioncenter.org/casp14/rrc_results.cgi, accessed on 25 April 2021. Table 7 compares the prediction accuracy of the COMTOP method with state-of-the-art systems on the CASP14 dataset. From Table 7, COMTOP performed better than tFold, MULTICOM-AI, RaptorX, trfold, and Kiharalab_Contact schemes for the top-L/2 and top-L/5 predictions, respectively. On the other hand, TripletRes and DeepPotential had better performance for top-L/2 and top-L/5 predictions.

3.5. Testing on the Independent TM Test Set

In this section, we tested COMTOP on 57 α-helical TM proteins using the proposed parameter sets. Table 8 compares the COMTOP prediction results to the seven individual approaches, revealing that the COMTOP model beat the others by a substantial margin. COMTOP obtained an average accuracy of 41.56%, 53.38%, 64.34%, and 73.91% for the top-2L, top-L, top-L/2, and top-L/5 predictions, respectively.
Among the seven methods, PconsC4 had the highest accuracy with 37.32%, 45.44%, 55.16%, and 65.26% for top-2L, top-L, top-L/2, and top-L/5 predictions, respectively. The DeepCov method had second-highest score of accuracy with 35.45%, 42.34%, 53.24%, and 63.03% for top-2L, top-L, top-L/2, and top-L/5 predictions, respectively. On the other hand, the PSICOV method had the lowest accuracy with 24.75%, 25.21%, 26.83%, and 29.17% for top-2L, top-L, top-L/2, and top-L/5 predictions, respectively. The prediction accuracies of all the methods were ranked as follows: PconsC4 > DeepCov > CCMpred > plmDCA > EVfold > NNcon > PSICOV, whereas the ranking of all methods based on test dataset was as follows: PconsC4 > DeepCov > CCMpred > EVfold > plmDCA > PSICOV > NNcon. COMTOP’s prediction accuracy for top-L/5 contacts decreased by 23.44%, 2%, and 3.58% in the TM protein dataset from the 239 non-redundant proteins, CAS13, and CASP14 datasets, respectively.

4. Discussion

In this paper, the proposed hybrid framework COMTOP model used information from the seven individual methods that were different from each other in terms of both methodology and input features. The seven methods can be roughly classified into three different categories: traditional machine learning, evolutionary coupling analysis, and deep learning. These methods also rely on different input data types. As shown in our previous work [72], the prediction results of these methods show certain degrees of similarity and difference, and the differences of prediction results from methods in the different categories are larger than that in the same category. COMTOP selects the weight variables for each individual method through the MILP optimization-based approach. Our consensus method scores higher for the correct contact predictions than for the incorrect contact predictions from the different methods; thus, individual methods with higher accuracy usually obtain higher weights. The seven methods can complement each other in prediction performance, so our method produces higher accuracy compared with seven individual methods and shows better or close prediction performance when compared with other state-of-the-art methods.
These seven methods were selected among the available methods based on two criteria: methods that (1) have comparatively better prediction accuracy and (2) belong to different method categories and complement each other. These seven contact prediction methods are classified into two categories: (1) coevolution-derived information-based, and (2) machine-learning-based. Most of these coevolution-derived approaches have been used MSA as input, which can be generated by approaches such as PSI-Blast, HHblits, or Jackhmmer. Most of these machine-learning-based methods have accepted a wide range of features as input, including features involved with the local window of the amino acids, amino acid type information (polarity and acidic properties), and the protein itself. This includes features such as mutual information of sequence profiles, information about the amino acid type (polarity and acidic properties), sequence profiles, sequence separation length between the amino acids under consideration, secondary structure, solvent accessibility, and pairwise information between all the amino acids involved [73]. The commonly used machine learning techniques for contact prediction are hidden Markov model, support vector machines, shallow neural networks and deep learning techniques.
COMTOP uses the confidence scores of the predicted contact and the weights of the seven individual methods to determine the protein residue–residue contact. The contact prediction of COMTOP is based on the sum of the products between the confidence score and the weight term over all methods; thus, the consensus method can also provide the confidence scores of the prediction for each pair of residues. Although the results for COMTOP use seven individual methods, COMTOP consistently shows better performance than the seven individual methods for top-5L, top-3L, top-2L, top-L, top-L/2, and top-L/5 contacts. It can be seen in Table 2 that the prediction accuracy of our methods is about 24.89%, 23.44%, 19.01%, 10.75%, 5.10%, and 1.59% better than the best of the individual methods, PconsC4, for top-5L, top-3L, top-2L, top-L, top-L/2, and top-L/5 contacts, respectively.
We also tested the performance of our method on CASP13, CASP14, and a non-redundant TM protein test set. The prediction accuracies were 75.91%, 77.49%, and 73.91%, respectively, for top-L/5 contacts. The prediction accuracies of COMTOP for all test sets significantly outperformed those of the seven individual methods. Furthermore, when we compared our method with a few state-of-the-art methods, it can be observed that the RapterX method showed better performance than our model for the CASP13 dataset. The COMTOP model outperformed Yang_server, TripletRes, ResTriplet, and DNCON3 schemes for the top-L/2 and top-L/5 contact predictions for CASP13 targets. In the CASP14 dataset, COMTOP performs better than tFold, MULTICOM-AI, RaptorX, trfold, and Kiharalab_Contact schemes for the top-L/2 and top-L/5 contacts, respectively. On the other hand, TripletRes and DeepPotential have better performance for top-L/2 and top-L/5 contacts. The results and analysis demonstrate that COMTOP can significantly improve the performance of the individual methods; therefore, COMTOP is more robust against different types of test sets. COMTOP also showed better/comparable predictions when compared with the state-of-the-art predictors.

5. Conclusions

In this paper, we presented a novel hybrid consensus method named as COMTOP and based on seven methods, aiming to predict high-quality protein contacts that can be used for 3D structure prediction. This consensus contact prediction method is based on a MILP model that produces the parameters for protein residue–residue contact prediction. The test on the 239 targets showed that COMTOP performed well compared with seven individual methods according to the prediction accuracy. COMTOP achieved a prediction accuracy of 75.91%, 77.49%, and 73.91% for top-L/5 contacts test on the CASP13, CASP14, and 50 TM target proteins, respectively, and showed satisfactory results compared with the state-of-the-art predictors. For all test datasets, the improvement of COMTOP in accuracy over the seven individual methods increases with the increasing number of predicted contacts. For example, COMTOP performs much better with a large number of contact predictions (such as top-5L and top-3L) than for a small number of contact predictions, such as top-L/2 and top-L/5.

Supplementary Materials

The following are available online at https://www.mdpi.com/article/10.3390/membranes11070503/s1, Table S1: The non-redundant list with 3133 proteins, Table S2: The training set of COMTOP, Table S3: 239 non-redundant proteins used for testing, Table S4: 30 CASP13 target domains used for testing, Table S5: 25 CASP14 target domains list used for testing, Table S6: 57 non-redundant TM proteins used for testing, Table S7: Parameters λ m 1 i (for the seven method COMTOP) obtained from the overall training process are shown in columns 1 to 6, respectively.

Author Contributions

M.S.R., H.Z. and Y.W. designed the study. M.S.R., S.F. and Hossain, M.T.H. developed the method, collected data, performed computations, and tested the method. H.Z. and L.J. collected data, calculated the scores for seven methods, and analyzed the data. Y.W. supervised the work. All authors have read and agreed to the published version of the manuscript.

Funding

This work was partly supported by the National Key Research and Development Program of China under Grant No. 2018YFB0204403; Strategic Priority CAS Project XDB38000000; National Science Foundation of China under grant no. U1813203; the Shenzhen Basic Research Fund under grant no. JCYJ20200109114818703, RCYX2020071411473419, and JSGG20201102163800001; CAS Key Lab under grant no. 2011DP173015 (YW). We would also like to thank the funding support from the Youth Innovation Promotion Association, CAS to Yanjie Wei, and Outstanding Youth Innovation Fund (CAS-SIAT to Huiling Zhang).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

Abbreviations

MILPMixed integer linear programming
CASPCritical assessment of protein structure prediction
PPVPositive predictive value
LLength of the protein sequence
PDBProtein Data Bank
TPTrue positive
FPFalse positive
MSAMultiple sequence alignment
TMTransmembrane protein

References

  1. Miller, C.S.; Eisenberg, D. Using inferred residue contacts to distinguish between correct and incorrect protein models. Bioinformatics 2008, 24, 1575–1582. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  2. Wang, Z.; Eickholt, J.; Cheng, J. APOLLO: A quality assessment service for single and multiple protein models. Bioinformatics 2011, 27, 1715–1716. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  3. Kliger, Y.; Levy, O.; Oren, A.; Ashkenazy, H.; Tiran, Z.; Novik, A.; Rosenberg, A.; Amir, A.; Wool, A.; Toporik, A.; et al. Peptides modulating conformational changes in secreted chaperones: From in silico design to preclinical proof of concept. Proc. Natl. Acad. Sci. USA 2009, 106, 13797–13801. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  4. Marks, D.S.; Colwell, L.J.; Sheridan, R.; Hopf, T.A.; Pagnani, A.; Zecchina, R.; Sander, C. Protein 3D Structure Computed from Evolutionary Sequence Variation. PLoS ONE 2011, 6, e28766. [Google Scholar] [CrossRef] [Green Version]
  5. Nugent, T.; Jones, D.T. Accurate de novo structure prediction of large transmembrane protein domains using fragment-assembly and correlated mutation analysis. Proc. Natl. Acad. Sci. USA 2012, 109, E1540–E1547. [Google Scholar] [CrossRef] [Green Version]
  6. Sułkowska, J.I.; Morcos, F.; Weigt, M.; Hwa, T.; Onuchic, J.N. Genomics-aided structure prediction. Proc. Natl. Acad. Sci. USA 2012, 109, 10340–10345. [Google Scholar] [CrossRef] [Green Version]
  7. Kim, D.E.; Di Maio, F.; Wang, R.Y.-R.; Song, Y.; Baker, D. One contact for every twelve residues allows robust and accurate topology-level protein structure modeling. Prot. Struct. Funct. Bioinform. 2014, 82, 208–218. [Google Scholar] [CrossRef] [Green Version]
  8. Kosciolek, T.; Jones, D.T. De Novo Structure Prediction of Globular Proteins Aided by Sequence Variation-Derived Contacts. PLoS ONE 2014, 9, e92197. [Google Scholar] [CrossRef] [Green Version]
  9. Yang, J.; Yan, R.; Roy, A.; Xu, D.; Poisson, J.; Zhang, Y. The I-TASSER Suite: Protein structure and function prediction. Nat. Methods 2015, 12, 7–8. [Google Scholar] [CrossRef] [Green Version]
  10. Ovchinnikov, S.; Kim, D.E.; Wang, R.Y.-R.; Liu, Y.; Di Maio, F.; Baker, D. Improved de novo structure prediction in CASP11 by incorporating coevolution information into Rosetta. Proteins Struct. Funct. Bioinform. 2016, 84, 67–75. [Google Scholar] [CrossRef] [Green Version]
  11. Vendruscolo, M.; Domany, E. Protein folding using contact maps. Vitam. Hormon. 2000, 58, 171–212. [Google Scholar] [CrossRef]
  12. Mirny, L.; Domany, E. Protein fold recognition and dynamics in the space of contact maps. Prot. Struct. Funct. Genet. 1996, 26, 391–410. [Google Scholar] [CrossRef]
  13. Ekeberg, M.; Lövkvist, C.; Lan, Y.; Weigt, M.; Aurell, E. Improved contact prediction in proteins: Using pseudolikelihoods to infer Potts models. Phys. Rev. E 2013, 87, 012707. [Google Scholar] [CrossRef] [Green Version]
  14. Eickholt, J.; Cheng, J. A study and benchmark of DNcon: A method for protein residue-residue contact prediction using deep networks. BMC Bioinform. 2013, 14, S12. [Google Scholar] [CrossRef] [Green Version]
  15. Seemayer, S.; Gruber, M.; Söding, J. CCMpred—Fast and precise prediction of protein residue–residue contacts from correlated mutations. Bioinformatics 2014, 30, 3128–3130. [Google Scholar] [CrossRef] [Green Version]
  16. Jones, D.T.; Singh, T.; Kosciolek, T.; Tetchner, S. MetaPSICOV: Combining coevolution methods for accurate prediction of contacts and long range hydrogen bonding in proteins. Bioinformatics 2015, 31, 999–1006. [Google Scholar] [CrossRef]
  17. Zhang, H.; Gao, Y.; Deng, M.; Wang, C.; Zhu, J.; Li, S.C.; Zheng, W.-M.; Bu, D. Improving residue–residue contact prediction via low-rank and sparse decomposition of residue correlation matrix. Biochem. Biophys. Res. Commun. 2016, 472, 217–222. [Google Scholar] [CrossRef]
  18. Zhang, H.; Huang, Q.; Bei, Z.; Wei, Y.; Floudas, C.A. COMSAT: Residue contact prediction of transmembrane proteins based on support vector machines and mixed integer linear programming. Proteins Struct. Funct. Bioinform. 2016, 84, 332–348. [Google Scholar] [CrossRef]
  19. Jones, D.T.; Kandathil, S.M. High precision in protein contact prediction using fully convolutional neural networks and minimal sequence features. Bioinformation 2018, 34, 3308–3315. [Google Scholar] [CrossRef]
  20. Michel, M.; Hurtado, D.M.; Elofsson, A. PconsC4: Fast, accurate and hassle-free contact predictions. Bioinformation 2019, 35, 2677–2679. [Google Scholar] [CrossRef]
  21. Jones, D.T.; Buchan, D.W.A.; Cozzetto, D.; Pontil, M. PSICOV: Precise structural contact prediction using sparse inverse covariance estimation on large multiple sequence alignments. Bioinformation 2011, 28, 184–190. [Google Scholar] [CrossRef] [Green Version]
  22. Kamisetty, H.; Ovchinnikov, S.; Baker, D. Assessing the utility of coevolution-based residue-residue contact predictions in a sequence- and structure-rich era. Proc. Natl. Acad. Sci. USA 2013, 110, 15674–15679. [Google Scholar] [CrossRef] [Green Version]
  23. Kaján, L.; Hopf, T.A.; Kalaš, M.; Marks, D.S.; Rost, B. FreeContact: Fast and free software for protein contact prediction from residue co-evolution. BMC Bioinform. 2014, 15, 85. [Google Scholar] [CrossRef] [Green Version]
  24. Baldassi, C.; Zamparo, M.; Feinauer, C.; Procaccini, A.; Zecchina, R.; Weigt, M.; Pagnani, A. Fast and Accurate Multivariate Gaussian Modeling of Protein Families: Predicting Residue Contacts and Protein-Interaction Partners. PLoS ONE 2014, 9, e92721. [Google Scholar] [CrossRef]
  25. Björkholm, P.; Daniluk, P.; Kryshtafovych, A.; Fidelis, K.; Andersson, R.; Hvidsten, T.R. Using multi-data hidden Markov models trained on local neighborhoods of protein structure to predict residue–residue contacts. Bioinformation 2009, 25, 1264–1270. [Google Scholar] [CrossRef] [Green Version]
  26. Cheng, J.; Baldi, P. Improved residue contact prediction using support vector machines and a large feature set. BMC Bioinform. 2007, 8, 113. [Google Scholar] [CrossRef] [Green Version]
  27. Wu, S.; Zhang, Y. A comprehensive assessment of sequence-based and template-based methods for protein contact prediction. Bioinformation 2008, 24, 924–931. [Google Scholar] [CrossRef] [Green Version]
  28. Tegge, A.N.; Wang, Z.; Eickholt, J.; Cheng, J. NNcon: Improved protein contact map prediction using 2D-recursive neural networks. Nucleic Acids Res. 2009, 37, W515–W518. [Google Scholar] [CrossRef] [PubMed]
  29. Wang, S.; Sun, S.; Li, Z.; Zhang, R.; Xu, J. Accurate De Novo Prediction of Protein Contact Map by Ultra-Deep Learning Model. PLoS Comput. Biol. 2017, 13, e1005324. [Google Scholar] [CrossRef] [Green Version]
  30. Adhikari, B.; Hou, J.; Cheng, J. DNCON2: Improved protein contact prediction using two-level deep convolutional neural networks. Bioinformation 2018, 34, 1466–1472. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  31. Hanson, J.; Paliwal, K.; Litfin, T.; Yang, Y.; Zhou, Y. Accurate prediction of protein contact maps by coupling residual two-dimensional bidirectional long short-term memory with convolutional neural networks. Bioinformation 2018, 34, 4039–4045. [Google Scholar] [CrossRef]
  32. Liu, Y.; Palmedo, P.; Ye, Q.; Berger, B.; Peng, J. Enhancing Evolutionary Couplings with Deep Convolutional Neural Networks. Cell Syst. 2018, 6, 65–74.e3. [Google Scholar] [CrossRef] [Green Version]
  33. Xiong, D.; Zeng, J.; Gong, H. A deep learning framework for improving long-range residue–residue contact prediction using a hierarchical strategy. Bioinformation 2017, 33, 2675–2683. [Google Scholar] [CrossRef] [Green Version]
  34. Golkov, V.; Skwark, M.J.; Golkov, A.; Dosovitskiy, A.; Brox, T.; Meiler, J.; Cremers, D. Protein contact prediction from amino acid co-evolution using convolutional networks for graph-valued images. In Proceedings of the 30th International Conference on Neural Information Processing Systems (NIPS’16); Curran Associates Inc.: Red Hook, NY, USA, 2016; pp. 4222–4230. [Google Scholar]
  35. Liu, Y.; Ye, Q.; Wang, L.; Peng, J. Learning structural motif representations for efficient protein structure search. Bioinformation 2018, 34, i773–i780. [Google Scholar] [CrossRef] [Green Version]
  36. McAllister, S.R.; Floudas, C.A. An improved hybrid global optimization method for protein tertiary structure prediction. Comput. Optim. Appl. 2010, 45, 377–413. [Google Scholar] [CrossRef] [Green Version]
  37. Rajgaria, R.; Wei, Y.; Floudas, C.A. Contact prediction for beta and alpha-beta proteins using integer linear optimization and its impact on the first principles 3D structure prediction method ASTRO-FOLD. Prot. Struct. Funct. Bioinform. 2010, 78, 1825–1846. [Google Scholar] [CrossRef] [Green Version]
  38. Wei, Y.; Floudas, C.A. Enhanced inter-helical residue contact prediction in transmembrane proteins. Chem. Eng. Sci. 2011, 66, 4356–4369. [Google Scholar] [CrossRef] [Green Version]
  39. Subramani, A.; Wei, Y.; Floudas, C.A. ASTRO-FOLD 2.0: An enhanced framework for protein structure prediction. AIChE J. 2011, 58, 1619–1637. [Google Scholar] [CrossRef] [Green Version]
  40. Zheng, W.; Li, Y.; Zhang, C.; Pearce, R.; Mortuza, S.M.; Zhang, Y. Deep-learning contact-map guided protein structure prediction in CASP. Proteins Struct. Funct. Bioinform. 2019, 87, 1149–1164. [Google Scholar] [CrossRef] [Green Version]
  41. Kandathil, S.M.; Greener, J.G.; Jones, D.T. Prediction of inter-residue contacts with DeepMetaPSICOV in CASP13. Proteins Struct. Funct. Bioinform. 2019, 87, 1092–1099. [Google Scholar] [CrossRef] [Green Version]
  42. Cheng, H.; Sen, T.Z.; Jernigan, R.L.; Kloczkowski, A. Consensus Data Mining (CDM) Protein Secondary Structure Prediction Server: Combining GOR V and Fragment Database Mining (FDM). Bioinformation 2007, 23, 2628–2630. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  43. Gupta, A.; Deshpande, A.; Sabarinathan, R.; Senthilkumar, R.; Sekar, K.; Amburi, J. CSSP(Consensus Secondary Structure Prediction): A web-based server for structural biologists. J. Appl. Crystallogr. 2009, 42, 336–338. [Google Scholar] [CrossRef]
  44. Wei, Y.; Thompson, J.; Floudas, C.A. CONCORD: A consensus method for protein secondary structure prediction via mixed integer linear optimization. Proc. R. Soc. A 2011, 468, 831–850. [Google Scholar] [CrossRef] [Green Version]
  45. Toussi, C.A.; Haddadnia, J. Improving protein secondary structure prediction: The evolutionary optimized classification algorithms. Struct. Chem. 2019, 30, 1257–1266. [Google Scholar] [CrossRef]
  46. Xu, J.; Yu, L.; Li, M.; Chen, Y.-P.P.; Wong, L. Consensus fold recognition by predicted model quality. Series on Advances in Bioinformatics and Computational Biology; Imperial College Press: London, UK, 2005; Volume 1, pp. 73–83. [Google Scholar]
  47. Wu, S.; Zhang, Y. LOMETS: A local meta-threading-server for protein structure prediction. Nucleic Acids Res. 2007, 35, 3375–3382. [Google Scholar] [CrossRef] [Green Version]
  48. McGuffin, L.J.; Adiyaman, R.; Maghrabi, A.H.A.; Shuid, A.N.; Brackenridge, D.A.; Nealon, J.O.; Philomina, L.S. IntFOLD: An integrated web resource for high performance protein structure and function prediction. Nucleic Acids Res. 2019, 47, W408–W413. [Google Scholar] [CrossRef]
  49. Liu, B.; Li, C.-C.; Yan, K. DeepSVM-fold: Protein fold recognition by combining support vector machines and pairwise sequence similarity scores generated by deep learning networks. Briefings Bioinform. 2019, 21, 1733–1741. [Google Scholar] [CrossRef]
  50. Kumar, S.; Carugo, O. Consensus Prediction of Protein Conformational Disorder from Amino Acidic Sequence. Open Biochem. J. 2008, 2, 1–5. [Google Scholar] [CrossRef]
  51. Xue, B.; Dunbrack, R.L.; Williams, R.W.; Dunker, A.K.; Uversky, V.N. PONDR-FIT: A meta-predictor of intrinsically disordered amino acids. Biochim. Biophys. Acta Prot. Proteom. 2010, 1804, 996–1010. [Google Scholar] [CrossRef] [Green Version]
  52. Liu, Y.; Chen, S.; Wang, X.; Liu, B. Identification of Intrinsically Disordered Proteins and Regions by Length-Dependent Predictors Based on Conditional Random Fields. Mol. Ther. Nucleic Acids 2019, 17, 396–404. [Google Scholar] [CrossRef] [Green Version]
  53. Kang, S.; Yang, S.-J.; Kim, S.; Bhak, J. CONSORF: A consensus prediction system for prokaryotic coding sequences. Bioinformatics 2007, 23, 3088–3090. [Google Scholar] [CrossRef]
  54. Swift, S.; Tucker, A.; Vinciotti, V.; Martin, N.; Orengo, C.; Liu, X.; Kellam, P. Consensus clustering and functional interpretation of gene-expression data. Genome Biol. 2004, 5, R94. [Google Scholar] [CrossRef] [Green Version]
  55. Wei, Y.; Nadler, W.; Hansmann, U.H.E. Side-chain and backbone ordering in a polypeptide. J. Chem. Phys. 2006, 125, 164902. [Google Scholar] [CrossRef] [Green Version]
  56. Kar, P.; Wei, Y.; Hansmann, U.H.E.; Höfinger, S. Systematic study of the boundary composition in Poisson Boltzmann calculations. J. Comput. Chem. 2007, 28, 2538–2544. [Google Scholar] [CrossRef] [Green Version]
  57. Wei, Y.; Nadler, W.; Hansmann, U.H.E. Backbone and side-chain ordering in a small protein. J. Chem. Phys. 2008, 128, 025105. [Google Scholar] [CrossRef] [Green Version]
  58. Zhang, H.; Xi, W.; Hansmann, U.H.E.; Wei, Y. Fibril–Barrel Transitions in Cylindrin Amyloids. J. Chem. Theory Comput. 2017, 13, 3936–3944. [Google Scholar] [CrossRef]
  59. Ma, J.; Wang, S.; Wang, Z.; Xu, J. Protein contact prediction by integrating joint evolutionary coupling analysis and supervised learning. Bioinformation 2015, 31, 3506–3513. [Google Scholar] [CrossRef] [Green Version]
  60. Shiwang, L.; Li, W.; Zhang, R.; Liu, S.; Xu, J. CoinFold: A web server for protein contact prediction and contact-assisted protein folding. Nucleic Acids Res. 2016, 44, W361–W366. [Google Scholar] [CrossRef] [Green Version]
  61. Jumper, J.; Evans, R.; Pritzel, A.; Green, T.; Figurnov, M.; Tunyasuvunakool, K.; Ronneberger, O.; Bates, R.; Žídek, A.; Bridgland, A.; et al. High Accuracy Protein Structure Prediction Using Deep Learning. In Proceedings of the Fourteenth Critical Assessment of Techniques for Protein Structure Prediction (Abstract Book), Online, 30 November–4 December 2020. [Google Scholar]
  62. Wang, G.; Dunbrack, R.L. PISCES: A protein sequence culling server. Bioinformatics 2003, 19, 1589–1591. [Google Scholar] [CrossRef] [Green Version]
  63. Rowe, P. Essential Statistics for the Pharmaceutical Sciences. Essent. Stat. Pharmaceut. Sci. 2015, 12, 89–91. [Google Scholar]
  64. Jessop, A. Let the Evidence Speak: Using Bayesian Thinking in Law, Medicine, Ecology and Other Areas; Springer: New York, NY, USA, 2018; ISBN 9783319713922. [Google Scholar]
  65. Michel, M.; Hayat, S.; Skwark, M.J.; Sander, C.; Marks, D.S.; Elofsson, A. PconsFold: Improved contact predictions improve protein models. Bioinformatics 2014, 30, i482–i488. [Google Scholar] [CrossRef]
  66. Adhikari, B.; Bhattacharya, D.; Cao, R.; Cheng, J. CONFOLD: Residue-residue contact-guidedab initioprotein folding. Proteins Struct. Funct. Bioinform. 2015, 83, 1436–1449. [Google Scholar] [CrossRef] [Green Version]
  67. Wu, Q.; Peng, Z.; Anishchenko, I.; Cong, Q.; Baker, D.; Yang, J. Protein contact prediction using metagenome sequence data and residual neural networks. Bioinformation 2019, 36, 41–48. [Google Scholar] [CrossRef]
  68. Li, Y.; Zhang, C.; Bell, E.W.; Yu, D.; Zhang, Y. Ensembling multiple raw coevolutionary features with deep residual neural networks for contact-map prediction in CASP13. Proteins: Struct. Funct. Bioinform. 2019, 87, 1082–1091. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  69. Wu, T.; Hou, J.; Adhikari, B.; Cheng, J. Analysis of several key factors influencing deep learning-based inter-residue contact prediction. Bioinformation 2020, 36, 1091–1098. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  70. Chen, X.; Liu, J.; Guo, Z.; Wu, T.; Hou, J.; Cheng, J. Protein model accuracy estimation empowered by deep learning and inter-residue distance prediction in CASP14. Sci. Rep. 2021, 11, 1–12. [Google Scholar] [CrossRef]
  71. Venkata Subramaniya, S.R.M.; Terashi, G.; Jain, A.; Kagaya, Y.; Kihara, D. Protein Contact Map Denoising Using Generative Adversarial Networks. bioRxiv 2020, 1–26. [Google Scholar] [CrossRef]
  72. Zhang, H.; Bei, Z.; Xi, W.; Hao, M.; Ju, Z.; Saravanan, K.M.; Zhang, H.; Guo, N.; Wei, Y. Evaluation of residue-residue contact prediction methods: From retrospective to prospective. PLoS Comput. Biol. 2021, 17, e1009027. [Google Scholar] [CrossRef]
  73. Adhikari, B.; Cheng, J. Protein Residue Contacts and Prediction Methods. In Methods in Molecular Biology; Springer Science and Business Media LLC: Berlin/Heidelberg, Germany, 2016; Volume 1415, pp. 463–476. [Google Scholar]
Figure 1. The workflow of COMTOP model for protein contact prediction.
Figure 1. The workflow of COMTOP model for protein contact prediction.
Membranes 11 00503 g001
Figure 2. The weight values and the overall performance for the training set. (A) The weight values of PconsC4, DeepCov, CCMpred, EVfold, plmDCA, PSICOV, and NNcon with the variation of the sample size. (B) The average accuracy, coverage, specificity, and MCC of COMTOP at different sample sizes.
Figure 2. The weight values and the overall performance for the training set. (A) The weight values of PconsC4, DeepCov, CCMpred, EVfold, plmDCA, PSICOV, and NNcon with the variation of the sample size. (B) The average accuracy, coverage, specificity, and MCC of COMTOP at different sample sizes.
Membranes 11 00503 g002
Figure 3. The overall performance of the COMTOP model. (A): the prediction accuracy of COMTOP with seven individual methods against predicted number of protein contacts for 239 proteins test dataset. (B,D): Bland–Altman plots between COMTOP and DeepCov, with the representation of the limits of agreement (dotted line), from +2SD to −2SD. (C,E): Bland–Altman plots between COMTOP and PconsC4, with the representation of the limits of agreement (dotted line), from +2SD to −2SD.
Figure 3. The overall performance of the COMTOP model. (A): the prediction accuracy of COMTOP with seven individual methods against predicted number of protein contacts for 239 proteins test dataset. (B,D): Bland–Altman plots between COMTOP and DeepCov, with the representation of the limits of agreement (dotted line), from +2SD to −2SD. (C,E): Bland–Altman plots between COMTOP and PconsC4, with the representation of the limits of agreement (dotted line), from +2SD to −2SD.
Membranes 11 00503 g003
Figure 4. The average accuracy, coverage, specificity, NPV, and MCC of COMTOP in terms of the predicted contact size on the 239 non-redundant proteins.
Figure 4. The average accuracy, coverage, specificity, NPV, and MCC of COMTOP in terms of the predicted contact size on the 239 non-redundant proteins.
Membranes 11 00503 g004
Figure 5. The overall performance of the COMTOP model. (A): the prediction accuracy of COMTOP with seven individual methods against predicted number of protein contacts size for CASP13 dataset. (B,D): Bland–Altman plots between COMTOP and DeepCov, with the representation of the limits of agreement (dotted line), from +2SD to −2SD. (C,E): Bland–Altman plots between COMTOP and PconsC4, with the representation of the limits of agreement (dotted line), from +2SD to −2SD.
Figure 5. The overall performance of the COMTOP model. (A): the prediction accuracy of COMTOP with seven individual methods against predicted number of protein contacts size for CASP13 dataset. (B,D): Bland–Altman plots between COMTOP and DeepCov, with the representation of the limits of agreement (dotted line), from +2SD to −2SD. (C,E): Bland–Altman plots between COMTOP and PconsC4, with the representation of the limits of agreement (dotted line), from +2SD to −2SD.
Membranes 11 00503 g005
Figure 6. The average accuracy, coverage, specificity, NPV, and MCC of COMTOP at predicted contact sizes on CASP13 dataset.
Figure 6. The average accuracy, coverage, specificity, NPV, and MCC of COMTOP at predicted contact sizes on CASP13 dataset.
Membranes 11 00503 g006
Table 1. Jaccard’s similarity coefficient based on training dataset of seven methods.
Table 1. Jaccard’s similarity coefficient based on training dataset of seven methods.
NNconEVfoldplmDCAPSICOVCCMpredDeepCovPconsC4
Nncon1
EVfold0.0761
plmDCA0.0640.3531
PSICOV0.0680.4880.3511
CCMpred0.0790.5170.3060.5401
DeepCov0.2100.2150.1760.2430.2851
PconsC40.1940.2910.2070.3210.3750.5851
Table 2. The prediction accuracies for COMTOP and the seven individual methods (PconsC4, DeepCov, CCMpred, EVfold, plmDCA, PSICOV, and NNcon) on the 239 non-redundant proteins.
Table 2. The prediction accuracies for COMTOP and the seven individual methods (PconsC4, DeepCov, CCMpred, EVfold, plmDCA, PSICOV, and NNcon) on the 239 non-redundant proteins.
MethodsTop-5LTop-3LTop-2LTop-LTop-L/2Top-L/5
NNcon15.6618.2721.2427.4434.3443.18
PSICOV16.7520.2623.2828.6034.2642.71
plmDCA17.1124.3930.1640.9850.7259.26
EVfold17.5824.0429.8941.8054.6267.57
CCMpred22.2129.3436.7651.4363.9874.00
DeepCov29.8941.2751.9270.2183.2891.61
PconsC434.7947.3559.8578.2989.4195.76
COMTOP59.6870.7978.8689.0494.5197.35
Table 3. The prediction accuracy on CASP13 target domains for COMTOP with seven individual methods (PconsC4, DeepCov, CCMpred, EVfold, plmDCA, PSICOV, and NNcon).
Table 3. The prediction accuracy on CASP13 target domains for COMTOP with seven individual methods (PconsC4, DeepCov, CCMpred, EVfold, plmDCA, PSICOV, and NNcon).
MethodsTop-5LTop-3LTop-2LTop-LTop-L/2Top-L/5
PSICOV14.5014.9017.9021.8026.5032.90
NNcon10.5013.9017.4023.5029.1036.70
plmDCA12.1015.4018.7025.0032.0040.10
EVfold12.2015.8019.3025.6032.1039.00
CCMpred13.3017.2021.7029.3036.1042.30
PconsC421.5030.1037.8049.7058.9063.90
DeepCov21.4029.8038.1052.6064.4075.10
COMTOP48.5053.2060.1066.7073.9075.91
Table 4. Performance of COMTOP scheme in CASP13. Top-L/2 and top-L/5 overall accuracy is shown for 30 FM/TBM and FM domains. Targets are ordered by domain classification.
Table 4. Performance of COMTOP scheme in CASP13. Top-L/2 and top-L/5 overall accuracy is shown for 30 FM/TBM and FM domains. Targets are ordered by domain classification.
DomainLength of DomainClassification of DomainAccuracy of Top-L/2 (%)Accuracy of Top-L/5 (%)
T0950-D1342FM63.275
T0953s1-D167FM92.9100
T0953s2-D2127FM68.833.3
T0953s2-D377FM28.60
T0957s1-D1108FM61.150
T0957s2-D1155FM7070
T0960-D284FM63.760
T0963-D282FM80100
T0968s1-D1118FM96.8100
T0968s2-D1115FM100100
T0951-D1266TBM-easy97.8100
T0960-D5105TBM-easy100100
T0963-D594TBM-easy100100
T1003-D1434TBM-easy8889.1
T1016-D1202TBM-easy97100
T0954-D1336TBM-hard100100
T0957s1-D254TBM-hard33.350
T0960-D389TBM-hard100100
T0963-D393TBM-hard100100
T0966-D1492TBM-hard84.691.9
T1009-D1718TBM-hard92.7100
T1011-D1280TBM-hard95.8100
T0953s2-D144FM/TBM75100
T0958-D177FM/TBM84.6100
T1005-D1326FM/TBM89.791.3
T0960-D132not evaluated00
T0960-D464not evaluated38.533.3
T0963-D131not evaluated00
T0963-D464not evaluated3033.3
T1011-D2160not evaluated84.4100
Table 5. Comparison of COMTOP model performance with a few state-of-the-art systems for 20 domains.
Table 5. Comparison of COMTOP model performance with a few state-of-the-art systems for 20 domains.
MethodsTop-L/2Top-L/5
RapterX85.9293.37
COMTOP84.0288.87
Yang_server77.1692.24
TripletRes76.9888.61
ResTriplet75.4787.63
DNCON352.2464.87
Table 6. The prediction accuracy of COMTOP with seven individual methods on CASP14 target domains.
Table 6. The prediction accuracy of COMTOP with seven individual methods on CASP14 target domains.
MethodsTop-5LTop-3LTop-2LTop-LTop-L/2Top-L/5
PSICOV08.5809.2309.6111.7612.8814.46
plmDCA08.0608.3509.0011.7814.3819.95
EVfold09.9010.9011.4810.8613.4917.36
CCMpred09.0210.2911.8914.9718.9224.34
NNcon17.6118.8819.7327.5434.4743.99
PconsC422.3224.7630.8741.5065.6269.46
DeepCov22.0032.1938.2448.6467.6571.86
COMTOP41.2746.9253.1061.1868.3377.49
Table 7. Comparison of COMTOP’s performance with a few state-of-the-art methods on the CASP14 dataset.
Table 7. Comparison of COMTOP’s performance with a few state-of-the-art methods on the CASP14 dataset.
MethodsTop-L/2Top-L/5
TripletRes76.0783.45
DeepPotential73.3181.23
COMTOP68.3377.49
tFold67.5677.11
MULTICOM-AI66.6774.93
RaptorX62.3470.12
trfold56.7466.04
Kiharalab_Contact53.1159.87
Table 8. The prediction accuracy of COMTOP with seven individual methods for 57 α-helical TM proteins dataset.
Table 8. The prediction accuracy of COMTOP with seven individual methods for 57 α-helical TM proteins dataset.
Methodstop-2Ltop-Ltop-L/2top-L/5
PSICOV24.7525.2126.8329.17
NNcon23.6924.6126.4129.22
EVfold27.8229.8932.9437.62
PlmDCA27.9830.9535.5942.83
CCMpred30.6534.8340.8847.74
DeepCov35.4542.3453.2463.03
PconsC437.3245.4455.1665.26
COMTOP41.5653.3864.3473.91
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Reza, M.S.; Zhang, H.; Hossain, M.T.; Jin, L.; Feng, S.; Wei, Y. COMTOP: Protein Residue–Residue Contact Prediction through Mixed Integer Linear Optimization. Membranes 2021, 11, 503. https://doi.org/10.3390/membranes11070503

AMA Style

Reza MS, Zhang H, Hossain MT, Jin L, Feng S, Wei Y. COMTOP: Protein Residue–Residue Contact Prediction through Mixed Integer Linear Optimization. Membranes. 2021; 11(7):503. https://doi.org/10.3390/membranes11070503

Chicago/Turabian Style

Reza, Md. Selim, Huiling Zhang, Md. Tofazzal Hossain, Langxi Jin, Shengzhong Feng, and Yanjie Wei. 2021. "COMTOP: Protein Residue–Residue Contact Prediction through Mixed Integer Linear Optimization" Membranes 11, no. 7: 503. https://doi.org/10.3390/membranes11070503

APA Style

Reza, M. S., Zhang, H., Hossain, M. T., Jin, L., Feng, S., & Wei, Y. (2021). COMTOP: Protein Residue–Residue Contact Prediction through Mixed Integer Linear Optimization. Membranes, 11(7), 503. https://doi.org/10.3390/membranes11070503

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