2014/10/31 by Benedikt Obermayer, Erel Levine
Biochemistry, Genetics and Molecular Biology · Physics and Astronomy · #A priori and a posteriori #Inference #Inverse #Ising model #Phylogenetic tree #Protein Structure and Dynamics #Sequence (biology) #Simple (philosophy) #Statistical Mechanics and Entropy #Statistical inference #Theoretical and Computational Physics #cond-mat.dis-nn #cond-mat.stat-mech #q-bio.PE
paper · pdf · doi:10.1088/1367-2630/16/12/123017
published as New J. Phys. 16:123017 (2014) · 18 pages, 6 figures; accepted at New J Phys
arxiv created 2014/10/31 · openalex publication_date 2014/12/08 · arxiv updated 2014/12/10 · openalex created_date 2016/06/24 · openalex updated_date 2026/08/06
Correlations between two variables of a high-dimensional system can be indicative of an underlying interaction, but can also result from indirect effects. Inverse Ising inference is a method to distinguish one from the other. Essentially, the parameters of the least constrained statistical model are learned from the observed correlations such that direct interactions can be separated from indirect correlations. Among many other applications, this approach has been helpful for protein structure prediction, because residues which interact in the 3D structure often show correlated substitutions in a multiple sequence alignment. In this context, samples used for inference are not independent but share an evolutionary history on a phylogenetic tree. Here, we discuss the effects of correlations between samples on global inference. Such correlations could arise due to phylogeny but also via other slow dynamical processes. We present a simple analytical model to address the resulting inference biases, and develop an exact method accounting for background correlations in alignment data by combining phylogenetic modeling with an adaptive cluster expansion algorithm. We find that popular reweighting schemes are only marginally effective at removing phylogenetic bias, suggest a rescaling strategy that yields better results, and provide evidence that our conclusions carry over to the frequently used mean-field approach to the inverse Ising problem.