<it>H2r</it>: Identification of evolutionary important residues by means of an entropy based analysis of multiple sequence alignments

• Main Authors: Zwick Matthias, Merkl Rainer
• Format: Article
• Language: English
• Published: BMC 2008-03-01
• Series: BMC Bioinformatics
• Online Access: http://www.biomedcentral.com/1471-2105/9/151

<p>Abstract</p> <p>Background</p> <p>A multiple sequence alignment (MSA) generated for a protein can be used to characterise residues by means of a statistical analysis of single columns. In addition to the examination of individual positions, the investigation of co-variation of amino acid frequencies offers insights into function and evolution of the protein and residues.</p> <p>Results</p> <p>We introduce <it>conn(k)</it>, a novel parameter for the characterisation of individual residues. For each residue <it>k</it>, <it>conn(k) </it>is the number of most extreme signals of co-evolution. These signals were deduced from a normalised mutual information (<it>MI</it>) value <it>U</it>(<it>k</it>, <it>l</it>) computed for all pairs of residues <it>k</it>, <it>l</it>. We demonstrate that <it>conn(k) </it>is a more robust indicator than an individual <it>MI</it>-value for the prediction of residues most plausibly important for the evolution of a protein. This proposition was inferred by means of statistical methods. It was further confirmed by the analysis of several proteins. A server, which computes <it>conn(k)</it>-values is available at <url>http://www-bioinf.uni-regensburg.de</url>.</p> <p>Conclusion</p> <p>The algorithms <it>H2r</it>, which analyses MSAs and computes <it>conn(k)</it>-values, characterises a specific class of residues. In contrast to strictly conserved ones, these residues possess some flexibility in the composition of side chains. However, their allocation is sensibly balanced with several other positions, as indicated by <it>conn(k)</it>.</p>