Using genome-wide protein sequence data to predict amino acid conservation.
For most proteins, multiple sequence alignments are a viable method to identify functionally and structurally important amino acids, but for most organisms, there is a subset of proteins that are unique or found in a few closely related organisms. For these proteins, it is not possible to produce sequence alignments that are useful in identifying functionally or structurally important amino acids. We have investigated the relationship between amino acid conservation and five factors (the amino acid's identity, N-terminal neighbor, C-terminal neighbor, the local hydropathy of surrounding amino acids, and the local expected net charge of the surrounding amino acids based on the primary sequence) in Escherichia coli proteins. For four of the factors examined (all but the amino acid's identity), there is a significant relationship with conservation for some of the standard 20 amino acids. Using the combination of all five factors, we show that it is possible to calculate a score based on the primary sequences of a subset of E. coli proteins that has statistically significant predictive value with respect to predicting conserved amino acids in other E. coli proteins and Saccharomyces cerevisiae proteins. As these five variables show significant relationships with conservation, we have termed them conservation factors.
Palenchar
Peter
P
Department of Chemistry, Rutgers University, 315 Penn St, Camden, NJ 08102-1411, USA. peterpal@crab.rutgers.edu
Mount
Mathew
M
Cusato
Douglas
D
Dougherty
Jeffery
J
D000595:Amino Acid Sequence; D000818:Animals; D019295:Computational Biology; D017124:Conserved Sequence; D030562:Databases, Protein; D004926:Escherichia coli / Q000737:chemistry*; D029968:Escherichia coli Proteins / Q000737:chemistry*; D008958:Models, Molecular; D011487:Protein Conformation; D012441:Saccharomyces cerevisiae / Q000737:chemistry*; D029701:Saccharomyces cerevisiae Proteins / Q000737:chemistry*; D016415:Sequence Alignment
D029968:Escherichia coli Proteins; D029701:Saccharomyces cerevisiae Proteins
Adv Enzymol Relat Areas Mol Biol. 1978;47:45-148
364941
BMC Bioinformatics. 2003 Sep 11;4:41
12969510
EXS. 2000;88:215-31
10803381
Protein J. 2008 Aug;27(5):283-91
18478318
Bioinformatics. 2001 May;17(5):455-60
11331240
Proc Natl Acad Sci U S A. 1989 Jan;86(2):699-703
2521391
Comput Biol Med. 2006 Oct;36(10):1143-54
16253222
J Mol Biol. 1982 May 5;157(1):105-32
7108955
Biochim Biophys Acta. 1976 Aug 9;439(2):261-73
952963
BMC Bioinformatics. 2007 Sep 21;8:356
17888146