Benchmarking inverse statistical approaches for protein structure and design with exactly solvable models

Abstract

Inverse statistical approaches to determine protein structure and function from Multiple Sequence Alignments (MSA) are emerging as powerful tools in computational biology. However the underlying assumptions of the relationship between the inferred effective Potts Hamiltonian and real protein structure and energetics remain untested so far. Here we use lattice protein model (LP) to benchmark those inverse statistical approaches. We build MSA of highly stable sequences in target LP structures, and infer the effective pairwise Potts Hamiltonians from those MSA. We find that inferred Potts Hamiltonians reproduce many important aspects of 'true' LP structures and energetics. Careful analysis reveals that effective pairwise couplings in inferred Potts Hamiltonians depend not only on the energetics of the native structure but also on competing folds; in particular, the coupling values reflect both positive design (stabilization of native conformation) and negative design (destabilization of competing folds). In addition to providing detailed structural information, the inferred Potts models used as protein Hamiltonian for design of new sequences are able to generate with high probability completely new sequences with the desired folds, which is not possible using independent-site models. Those are remarkable results as the effective LP Hamiltonians used to generate MSA are not simple pairwise models due to the competition between the folds. Our findings elucidate the reasons for the success of inverse approaches to the modelling of proteins from sequence data, and their limitations.

Supplementary Information available at http://journals.plos.org/ploscompbiol/article?id=10.1371%2Fjournal.pcbi.1004889

References in corpus (11)

• 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
• Why Do Proteins Look Like Proteins?
• Coevolutionary landscape inference and the context-dependence of mutations in beta-lactamase TEM-1
• Adaptive Cluster Expansion for Inferring Boltzmann Machines with Noisy Data
• Positive and negative design in stability and thermal adaptation of natural proteins
• Structural Determinant of Protein Designability
• Adaptive cluster expansion for the inverse Ising problem: convergence, algorithm and tests
• The network of stabilizing contacts in proteins studied by coevolutionary data
• On the entropy of protein families

Cited by in corpus (15)

• Inverse Statistical Physics of Protein Sequences: A Key Issues Review
• Inferring interaction partners from protein sequences
• Influence of Multiple Sequence Alignment Depth on Potts Statistical Models of Protein Covariation
• The role of evolutionary selection in the dynamics of protein structure evolution
• Parsimonious evolutionary scenario for the origin of allostery and coevolution patterns in proteins
• On the entropy of protein families
• Direct Coupling Analysis of Epistasis in Allosteric Materials
• Barriers and Dynamical Paths in Alternating Gibbs Sampling of Restricted Boltzmann Machines
• Selection originating from protein stability/foldability: Relationships between protein folding free energy, sequence ensemble, and fitness
• Resummed mean-field inference for strongly coupled data
• Evolutionary Dynamics of a Lattice Dimer: a Toy Model for Stability vs. Affinity Trade-offs in Proteins
• Interpretable Neural Networks based classifiers for categorical inputs
• Configurational entropy, transition rates, and optimal interactions for rapid folding in coarse-grained model proteins
• Prediction of Structures and Interactions from Genome Information
• Boltzmann machine learning and regularization methods for inferring evolutionary fields and couplings from a multiple sequence alignment