Predicting protein contact map using evolutionary and physical constraints by integer programming (extended version)

TL;DR

This paper introduces PhyCMAP, a novel contact map prediction method that combines evolutionary and physical constraints using machine learning and integer linear programming, significantly improving accuracy and speed over existing methods.

Contribution

The paper presents a new approach integrating evolutionary and physical restraints via ILP, reducing solution space and enhancing contact map prediction accuracy.

Findings

Outperforms popular existing methods in accuracy.

Predicts contacts within minutes after sequence homolog search.

Effective regardless of the number of available sequence homologs.

Abstract

Motivation. Protein contact map describes the pairwise spatial and functional relationship of residues in a protein and contains key information for protein 3D structure prediction. Although studied extensively, it remains very challenging to predict contact map using only sequence information. Most existing methods predict the contact map matrix element-by-element, ignoring correlation among contacts and physical feasibility of the whole contact map. A couple of recent methods predict contact map based upon residue co-evolution, taking into consideration contact correlation and enforcing a sparsity restraint, but these methods require a very large number of sequence homologs for the protein under consideration and the resultant contact map may be still physically unfavorable. Results. This paper presents a novel method PhyCMAP for contact map prediction, integrating both evolutionary and physical restraints by machine learning and integer linear programming (ILP). The evolutionary restraints include sequence profile, residue co-evolution and context-specific statistical potential. The physical restraints specify more concrete relationship among contacts than the sparsity restraint. As such, our method greatly reduces the solution space of the contact map matrix and thus, significantly improves prediction accuracy. Experimental results confirm that PhyCMAP outperforms currently popular methods no matter how many sequence homologs are available for the protein under consideration. PhyCMAP can predict contacts within minutes after PSIBLAST search for sequence homologs is done, much faster than the two recent methods PSICOV and EvFold. See http://raptorx.uchicago.edu for the web server.