Rate4Site: an algorithmic tool for the identification of functional regions in proteins by surface mapping of evolutionary determinants within their homologues

Abstract

Motivation: A number of proteins of known three-dimensional (3D) structure exist, with yet unknown function. In light of the recent progress in structure determination methodology, this number is likely to increase rapidly. A novel method is presented here: ‘Rate4Site’, which maps the rate of evolution among homologous proteins onto the molecular surface of one of the homologues whose 3D-structure is known. Functionally important regions often correspond to surface patches of slowly evolving residues. Results: Rate4Site estimates the rate of evolution of amino acid sites using the maximum likelihood (ML) principle. The ML estimate of the rates considers the topology and branch lengths of the phylogenetic tree, as well as the underlying stochastic process. To demonstrate its potency, we study the Src SH2 domain. Like previously established methods, Rate4Site detected the SH2 peptide-binding groove. Interestingly, it also detected inter-domain interactions between the SH2 domain and the rest of the Src protein that other methods failed to detect. Availability: Rate4Site can be downloaded at: http://ashtoret.tau.ac.il/ It is implemented as a web server at: bioinfo.tau.ac.il/ConSurf Contact: tal@ism.ac.jprebell@ashtoret.tau.ac.ilfabian@ashtoret.tau.ac.ilbental@ashtoret.tau.ac.il Supplementary Information: Multiple sequence alignment of homologous SH2 domains, the corresponding phylogenetic tree and additional examples are available at http://ashtoret.tau.ac.il/~rebell Keywords: rate variation among sites; evolutionary conservation; protein evolution; maximum likelihood; SH2 domains. *To whom correspondence should be addressed. †These authors contributed equally.