Published January 9, 2025 | Version v1
Journal article Open

Investigating Statistical Conditions of Coevolutionary Signals that Enable Algorithmic Predictions of Protein Partners

  • 1. Laboratório de Biologia Teórica e Computacional (LBTC), Universidade de Brasília

Description

The directory includes a comprehensive set of scripts for running Genetic Algorithm simulations and performing parameter calculations, including I', I*, I₀, σ₀², and TP rate. All analyses are demonstrated using both orthologous and paralogous protein families.

Notes

The work described here  was supported in part by the National Council of Technological and Scientific Development CNPq [grant no 302089/2019-5 (WT)] and Fundação de Apoio a Pesquisa do Distrito Federal FAPDF [grant no 00193-00001721/2024-66 (WT)]. WT thanks CAPES for doctoral fellowship to JA (grant no 88887.826533/2023-00).

Abstract

This study examines the statistical conditions of coevolutionary signals that allow algorithmic predictions of protein partners based on amino acid sequences, rather than 3D structures. It introduces a Markov stochastic model that predicts the number of correct protein partners based on coevolutionary information. The model defines state probabilities using a Poisson mixture of normal distributions, with key parameters including the total number of protein sequences M , the coevolutionary information gap α , and variance σ₀² . The model suggests that algorithmic approaches that maximize coevolutionary information cannot effectively resolve partners in protein families with a large number of sequences M ≥100 . The model shows that true positive (TP) rates can be enhanced by disregarding mismatches among similar sequences. This approach allows a distinction, in terms of { α , σ₀² } , between optimized solutions with trivial errors and other degenerate solutions. Our findings enable the a priori classification of protein families where partners can be reliably predicted by ignoring trivial errors between similar sequences, advancing the understanding of coevolutionary models for large protein datasets.

Files

GA_simulations_and_parameters.zip

Files (4.8 MB)

Name Size Download all
md5:e3c917ed934068e05e90b848bb61b0ef
4.8 MB Preview Download

Additional details

Software

Programming language
Python