Identification and analysis of novel tandem repeats in the cell surface proteins of archaeal and bacterial genomes using computational tools.

We have identified four novel repeats and two domains in cell surface proteins encoded by the Methanosarcina acetivorans genome and in some archaeal and bacterial genomes. The repeats correspond to a certain number of amino acid residues present in tandem in a protein sequence and each repeat is characterized by conserved sequence motifs. These correspond to: (a) a 42 amino acid (aa) residue RIVW repeat; (b) a 45 aa residue LGxL repeat; (c) a 42 aa residue LVIVD repeat; and (d) a 54 aa residue LGFP repeat. The domains correspond to a certain number of aa residues in a protein sequence that do not comprise internal repeats. These correspond to: (a) a 200 aa residue DNRLRE domain; and (b) a 70 aa residue PEGA domain. We discuss the occurrence of these repeats and domains in the different proteins and genomes analysed in this work.

Combinations of turns in proteins.

We observed that beta- and gamma-turns in protein structure may be associated as peptides representing combinations of turns that span between nine and 26 amino acid residues along the polypeptide backbone chain and often correspond to loops in the protein structure. Around 475 peptides resulted from the analysis of a non-redundant data set corresponding to 248 protein crystal structures selected from the Protein Data Bank. Nearly 40% protein chains are associated with two or more peptides and the peptides with nine and 10 amino acid residues are more frequent. A maximum of four distinct peptides varying in number of amino acid residues were observed in at least 10 proteins along the same protein chain. Nearly 80% peptides comprise type IV beta-turns that are associated with irregular dihedral angle values suggesting this may be important for the conformational diversity associated with the loops in proteins. In general, predominant interactions that possibly stabilize these peptides involve main-chain and side-chain interactions with solvent, in addition to hydrogen bond, salt-bridge and non-bonded interactions. Majority of the peptides were observed in hydrolase, oxidoreductase, transferase, serine proteinase/inhibitor complex, electron transport/electron transfer and lyase proteins.

Prediction of gamma-turns from amino acid sequences.

We predicted gamma-turns from amino acid sequences using the first-order Markov chain theory and enlarged representative data sets corresponding to protein chains selected from the Protein Data Bank (PDB). The following data sets were used for training and deriving the probability values: (1) an initial data set containing 315 protein chains comprising 904 gamma-turns and (2) a later data set in order to include new entries in the PDB, containing 434 protein chains and comprising 1053 gamma-turns. By excluding 93 protein chains that were common to these two training data sets, we generated two mutually exclusive data sets containing 222 and 341 protein chains for testing our predictions. Applying amino acid probability values derived from training data sets on to testing data sets yielded overall prediction accuracies in the range 54-57%. We recommend the use of probability values derived from the data set comprising 315 protein chains that represents more gamma-turns and also provides better predictions.