Fast and accurate multivariate Gaussian modeling of protein families: Predicting residue contacts and protein-interaction partners

Abstract

In the course of evolution, proteins show a remarkable conservation of their three-dimensional structure and their biological function, leading to strong evolutionary constraints on the sequence variability between homologous proteins. Our method aims at extracting such constraints from rapidly accumulating sequence data, and thereby at inferring protein structure and function from sequence information alone. Recently, global statistical inference methods (e.g. direct-coupling analysis, sparse inverse covariance estimation) have achieved a breakthrough towards this aim, and their predictions have been successfully implemented into tertiary and quaternary protein structure prediction methods. However, due to the discrete nature of the underlying variable (amino-acids), exact inference requires exponential time in the protein length, and efficient approximations are needed for practical applicability. Here we propose a very efficient multivariate Gaussian modeling approach as a variant of direct-coupling analysis: the discrete amino-acid variables are replaced by continuous Gaussian random variables. The resulting statistical inference problem is efficiently and exactly solvable. We show that the quality of inference is comparable or superior to the one achieved by mean-field approximations to inference with discrete variables, as done by direct-coupling analysis. This is true for (i) the prediction of residue-residue contacts in proteins, and (ii) the identification of protein-protein interaction partner in bacterial signal transduction. An implementation of our multivariate Gaussian approach is available at the website http://areeweb.polito.it/ricerca/cmp/code

arXiv comment24 pages, 7 pdf figures, 2 tables, plus supporting informations. Published on PLOS ONE

References in corpus (7)

• Direct-coupling analysis of residue co-evolution captures native contacts across many protein families
• Identification of direct residue contacts in protein-protein interaction by message passing
• Improved contact prediction in proteins: Using pseudolikelihoods to infer Potts models
• Maximum entropy models for antibody diversity
• From principal component to direct coupling analysis of coevolution in proteins: Low-eigenvalue modes are needed for structure prediction
• Dissecting the Specificity of Protein-Protein Interaction in Bacterial Two-Component Signaling: Orphans and Crosstalks
• Using Sequence Alignments to Predict Protein Structure and Stability With High Accuracy

Cited by in corpus (30)

• Distance-based Protein Folding Powered by Deep Learning
• Coevolutionary landscape inference and the context-dependence of mutations in beta-lactamase TEM-1
• Inverse Statistical Physics of Protein Sequences: A Key Issues Review
• Fast pseudolikelihood maximization for direct-coupling analysis of protein structure from many homologous amino-acid sequences
• Inferring interaction partners from protein sequences
• Simultaneous identification of specifically interacting paralogs and inter-protein contacts by Direct-Coupling Analysis
• Inter-protein sequence co-evolution predicts known physical interactions in bacterial ribosomes and the trp operon
• Inferring interaction partners from protein sequences using mutual information
• Selection of sequence motifs and generative Hopfield-Potts models for protein familiesilies
• Phylogenetic correlations can suffice to infer protein partners from sequences
• Revealing evolutionary constraints on proteins through sequence analysis
• Modeling sequence-space exploration and emergence of epistatic signals in protein evolution
• Inference of compressed Potts graphical models
• Identifying relevant positions in proteins by Critical Variable Selection
• A generative model for protein contact networks
• Correlations from structure and phylogeny combine constructively in the inference of protein partners from sequences
• Aligning biological sequences by exploiting residue conservation and coevolution
• Estimating probabilistic context-free grammars for proteins using contact map constraints
• Resummed mean-field inference for strongly coupled data
• Statistical Genetics in and out of Quasi-Linkage Equilibrium (Extended)
• Unified framework for modeling multivariate distributions in biological sequences
• Improved Pseudolikelihood Regularization and Decimation methods on Non-linearly Interacting Systems with Continuous Variables
• Inference for interacting linear waves in ordered and random media
• Graph Coding for Model Selection and Anomaly Detection in Gaussian Graphical Models
• Prediction of Structures and Interactions from Genome Information
• Ancestral Sequence Reconstruction for Co-evolutionary models
• Boltzmann machine learning and regularization methods for inferring evolutionary fields and couplings from a multiple sequence alignment
• Correlated evolution: models and methods
• Graph Compression with Application to Model Selection
• Small Coupling Expansion for Multiple Sequence Alignment