Correlation-compressed direct-coupling analysis

Research output: Contribution to journalArticleScientificpeer-review

Details

Original languageEnglish
Article number032407
Pages (from-to)1-15
JournalPhysical Review E
Volume98
Issue number3
Publication statusPublished - 11 Sep 2018
MoE publication typeA1 Journal article-refereed

Researchers

Research units

  • CAS - Institute of Theoretical Physics
  • Chinese Academy of Sciences
  • Hunan Normal University
  • KTH Royal Institute of Technology

Abstract

Learning Ising or Potts models from data has become an important topic in statistical physics and computational biology, with applications to predictions of structural contacts in proteins and other areas of biological data analysis. The corresponding inference problems are challenging since the normalization constant (partition function) of the Ising or Potts distribution cannot be computed efficiently on large instances. Different ways to address this issue have resulted in a substantial amount of methodological literature. In this paper we investigate how these methods could be used on much larger data sets than studied previously. We focus on a central aspect, that in practice these inference problems are almost always severely undersampled, and the operational result is almost always a small set of leading predictions. We therefore explore an approach where the data are prefiltered based on empirical correlations, which can be computed directly even for very large problems. Inference is only used on the much smaller instance in a subsequent step of the analysis. We show that in several relevant model classes such a combined approach gives results of almost the same quality as inference on the whole data set. It can therefore provide a potentially very large computational speedup at the price of only marginal decrease in prediction quality. We also show that the results on whole-genome epistatic couplings that were obtained in a recent computation-intensive study can be retrieved by our approach. The method of this paper hence opens up the possibility to learn parameters describing pairwise dependences among whole genomes in a computationally feasible and expedient manner.

Download statistics

No data available

ID: 28341522