Abstract:Several methods are available for the detection of covarying positions from a multiple sequence alignment (MSA). If the MSA contains a large number of sequences, information about the proximities between residues derived from covariation maps can be sufficient to predict a protein fold. If the structure is already known, information on the covarying positions can be valuable to understand the protein mechanism. In this study we have sought to determine whether a multivariate extension of traditional mutual information (MI) can be an additional tool to study covariation. The performance of two multidimensional MI (mdMI) methods, designed to remove the effect of ternary/quaternary interdependencies, was tested with a set of 9 MSAs each containing <400 sequences, and was shown to be comparable to that of methods based on maximum entropy/pseudolikelyhood statistical models of protein sequences. However, while all the methods tested detected a similar number of covarying pairs among the residues separated by < 8 Å in the reference X-ray structures, there was on average less than 65% overlap between the top scoring pairs detected by methods that are based on different principles. We have also attempted to identify whether the difference in performance among methods is due to different efficiency in removing covariation originating from chains of structural contacts. We found that the reason why methods that derive partial correlation between the columns of a MSA provide a better recognition of close contacts is not because they remove chaining effects, but because they filter out the correlation between distant residues that originates from general fitness constraints. In contrast we found that true chaining effects are expression of real physical perturbations that propagate inside proteins, and therefore are not removed by the derivation of partial correlation between variables.

Analysis of the Co-evolving Residues in Proteins Using Mutual Information Theory

Gene Prediction by the Noise-Assisted MEMD and Wavelet Transform for Identifying the Protein Coding Regions

Underestimation-Assisted Global-Local Cooperative Differential Evolution and the Application to Protein Structure Prediction

Combining Evolutionary Information Extracted From Frequency Profiles With Sequence-Based Kernels For Protein Remote Homology Detection

Sequence alignment and mutual information

Multidimensional mutual information methods for the analysis of covariation in multiple sequence alignments

Algorithm‐based Coevolution Network Identification Reveals Key Functional Residues of the Α/β Hydrolase Subfamilies

Efficiently Predicting Mutational Effect on Homologous Proteins by Evolution Encoding

Functionally Important Residues from Graph Analysis of Coevolved Dynamic couplings

Predicting Pathology of Missense Mutations through Protein-Specific Evolutionary Pattern

Protein Remote Homology Detection by Combining Profile-based Protein Representation with Local Alignment Kernel

A New Method for Quantifying Residue Conservation and Its Applications to the Protein Folding Nucleus

Simultaneous identification of specifically interacting paralogs and inter-protein contacts by Direct-Coupling Analysis

Analysis and prediction of protein folding energy changes upon mutation by element specific persistent homology

Inter-protein sequence co-evolution predicts known physical interactions in bacterial ribosomes and the trp operon

Aligning biological sequences by exploiting residue conservation and coevolution

How pairwise coevolutionary models capture the collective residue variability in proteins

Identifying the Genetic Basis of Functional Protein Evolution Using Reconstructed Ancestors

Predicting Protein Interaction Sites from Residue Spatial Sequence Profile and Evolution Rate

Inferring interaction partners from protein sequences using mutual information

Combining protein sequence, structure, and dynamics: A novel approach for functional evolution analysis of PAS domain superfamily