Identification of thermophilic species by the amino acid compositions deduced from their genomes

The global amino acid compositions as deduced from the complete genomic sequences of six thermophilic archaea, two thermophilic bacteria, 17 mesophilic bacteria and two eukaryotic species were analysed by hierarchical clustering and principal components analysis. Both methods showed an influence of several factors on amino acid composition. Although GC content has a dominant effect, thermophilic species can be identified by their global amino acid compositions alone. This study presents a careful statistical analysis of factors that affect amino acid composition and also yielded specific features of the average amino acid composition of thermophilic species. Moreover, we introduce the first example of a 'compositional tree' of species that takes into account not only homologous proteins, but also proteins unique to particular species. We expect this simple yet novel approach to be a useful additional tool for the study of phylogeny at the genome level.     

Molecular cloning and nucleotide sequences of the complementary DNAs to chicken skeletal muscle myosin two alkali light chain mRNAs.

We report here the molecular cloning and sequence analysis of DNAs complementary to mRNAs for myosin alkali light chain of chicken embryo and adult leg skeletal muscle. pSMA2-1 contained an 818 base-pair insert that includes the entire coding region and 5' and 3' untranslated regions of A2 mRNA. pSMA1-1 contained a 848 base-pair insert that included the 3' untranslated region and almost all of the coding region except for the N-terminal 13 amino acid residues of the A1 light chain. The 741 nucleotide sequences of A1 and A2 mRNAs corresponding to C-terminal 141 amino acid residues and 3' untranslated regions were identical. The 5' terminal nucleotide sequences corresponding to N-terminal 35 amino acid residues of A1 chain were quite different from the sequences corresponding to N-terminal 8 amino acid residues and of the 5' untranslated region of A2 mRNA. These findings are discussed in relation to the structures of the genes for A1 and A2 mRNA.     

Does spiralin has a cleavable N-terminal signal sequence?

The amino acid sequence of spiralin deduced from the nucleotide sequence of its gene was fictitiously shortened by 1 to 50 residues from each terminus and the compositions of both series of theoretical polypeptides were calculated. The two series of compositions thus obtained were compared to that of the purified protein, with the use of the Marchalonis and Weltman index (S delta Q). The results of this analysis, which permits the difficulty resulting from the blocking of the N-terminal amino acid to be overcome, show that spiralin is probably synthesized as a 241-residue precursor containing an N-terminal signal sequence cleaved close to cysteine-24. Since spiralin is acylated and since the sequence Val-Val-Ala-Cys24 shares some similarity with the consensus sequence of bacterial lipoprotein modification/processing site, the hypothesis of a cleavage just before cysteine-24 seems plausible.     

Comparative analysis of the base biases at the gene terminal portions in seven eukaryote genomes

Adenine nucleotides have been found to appear preferentially in the regions after the initiation codons or before the termination codons of bacterial genes. Our previous experiments showed that AAA and AAT, the two most frequent second codons in Escherichia coli, significantly enhance translation efficiency. To determine whether such a characteristic feature of base frequencies exists in eukaryote genes, we performed a comparative analysis of the base biases at the gene terminal portions using the proteomes of seven eukaryotes. Here we show that the base appearance at the codon third positions of gene terminal regions is highly biased in eukaryote genomes, although the codon third positions are almost free from amino acid preference. The bias changes depending on its position in a gene, and is characteristic of each species. We also found that bias is most outstanding at the second codon, the codon after the initiation codon. NCN is preferred in every genome; in particular, GCG is strongly favored in human and plant genes. The presence of the bias implies that the base sequences at the second codon affect translation efficiency in eukaryotes as well as bacteria.     

Terminal disorder: a common structural feature of the axial proteins of bacterial flagellum?

We report, based on proteolytic experiments and high resolution 1H nuclear magnetic resonance studies that the terminal regions of the monomeric hook protein are highly mobile and exposed to the solvent. The disordered parts of the hook protein span approximately the first 70 and the last 30 amino acid residues. Although the amino acid sequences of flagellin and hook protein do not resemble each other at all, both proteins have now been shown to contain large disordered terminal regions. Sequential similarities of flagellin and hook protein, especially near the NH2 and COOH termini, to other axial components of bacterial flagellum suggest that terminal disorder may be a common structural feature of the axial proteins of the bacterial flagellum.     

Comparative sequence analysis of acid sensitive/resistance proteins in Escherichia coli and Shigella flexneri

The molecular basis for the survival of bacteria under extreme conditions in which growth is inhibited is a question of great current interest. A preliminary study was carried out to determine residue pattern conservation among the antiporters of enteric bacteria, responsible for extreme acid sensitivity especially in Escherichia coli and Shigella flexneri. Here we found the molecular evidence that proved the relationship between E. coli and S. flexneri. Multiple sequence alignment of the gadC coded acid sensitive antiporter showed many conserved residue patterns at regular intervals at the N-terminal region. It was observed that as the alignment approaches towards the C-terminal, the number of conserved residues decreases, indicating that the N-terminal region of this protein has much active role when compared to the carboxyl terminal. The motif, FHLVFFLLLGG, is well conserved within the entire gadC coded protein at the amino terminal. The motif is also partially conserved among other antiporters (which are not coded by gadC) but involved in acid sensitive/resistance mechanism. Phylogenetic cluster analysis proves the relationship of Escherichia coli and Shigella flexneri. The gadC coded proteins are converged as a clade and diverged from other antiporters belongs to the amino acid-polyamine-organocation (APC) superfamily.     

Mapping of the DNA linking tyrosine residue of the PRD1 terminal protein.

DNA replication of PRD1, a lipid-containing phage, is initiated by a protein-priming mechanism. The terminal protein encoded by gene 8 acts as a protein primer in DNA synthesis by forming an initiation complex with the 5'-terminal nucleotide, dGMP. The linkage between the terminal protein and the 5' terminal nucleotide is a tyrosylphosphodiester bond. The PRD1 terminal protein contains 13 tyrosine residues in a total of 259 amino acids. By site-directed mutagenesis of cloned PRD1 gene 8, we replaced 12 of the 13 tyrosine residues in the terminal protein with phenylalanine and the other tyrosine residue with asparagine. Functional analysis of these mutant terminal proteins suggested that tyrosine-190 is the linking amino acid that forms a covalent bond with dGMP. Cyanogen bromide cleavage studies also implicated tyrosine-190 as the DNA-linking amino acid residue of the PRD1 terminal protein. Our results further show that tyrosine residues at both the amino-terminal and the carboxyl-terminal regions are important for the initiation complex forming activity. Predicted secondary structures for the regions around the DNA linking amino acid residues were compared in three terminal proteins (phi 29, adenovirus-2, and PRD1). While the linking amino acids serine-232 (phi 29) and serine-577 (adenovirus-2) are found in beta-turns in hydrophilic regions, the linking tyrosine-190 of the PRD1 terminal protein is found in a beta-sheet in a hydrophobic region.     

Complete Amino Acid Sequence of Globin Chains and Biological Activity of Fragmented Crocodile Hemoglobin (Crocodylus siamensis)

Hemoglobin, ??-chain, ??-chain and fragmented hemoglobin of Crocodylus siamensis demonstrated both antibacterial and antioxidant activities. Antibacterial and antioxidant properties of the hemoglobin did not depend on the heme structure but could result from the compositions of amino acid residues and structures present in their primary structure. Furthermore, thirteen purified active peptides were obtained by RP-HPLC analyses, corresponding to fragments in the ??-globin chain and the ??-globin chain which are mostly located at the N-terminal and C-terminal parts. These active peptides operate on the bacterial cell membrane. The globin chains of Crocodylus siamensis showed similar amino acids to the sequences of Crocodylus niloticus. The novel amino acid substitutions of ??-chain and ??-chain are not associated with the heme binding site or the bicarbonate ion binding site, but could be important through their interactions with membranes of bacteria.     

The interacting domains of three MutL heterodimers in man: hMLH1 interacts with 36 homologous amino acid residues within hMLH3, hPMS1 and hPMS2

In human cells, hMLH1, hMLH3, hPMS1 and hPMS2 are four recognised and distinctive homologues of MutL, an essential component of the bacterial DNA mismatch repair (MMR) system. The hMLH1 protein forms three different heterodimers with one of the other MutL homologues. As a first step towards functional analysis of these molecules, we determined the interacting domains of each heterodimer and tried to understand their common features. Using a yeast two-hybrid assay, we show that these MutL homologues can form heterodimers by interacting with the same amino acid residues of hMLH1, residues 492–742. In contrast, three hMLH1 partners, hMLH3, hPMS1 and hPMS2 contain the 36 homologous amino acid residues that interact strongly with hMLH1. Contrary to the previous studies, these homologous residues reside at the N-terminal regions of three subdomains conserved in MutL homologues in many species. Interestingly, these residues in hPMS2 and hMLH3 may form coiled-coil structures as predicted by the MULTICOIL program. Furthermore, we show that there is competition for the interacting domain in hMLH1 among the three other MutL homologues. Therefore, the quantitative balance of these three MutL heterodimers may be important in their functions.     

Structural proteins of Western equine encephalitis virus: amino acid compositions and N-terminal sequences

The structural proteins of Western equine encephalitis virus, a member of the alphavirus group, have been characterized by the determination of their amino acid compositions and by N-terminal sequence analysis. More than 60 residues of the N-terminal sequences of each of the envelope glycoproteins have been determined. A comparison of these sequences with the previously determined sequences of two related alphaviruses. Sindbis virus and Semliki Forest virus, strongly supports the view that all three viruses have evolved from a common ancestor and provides information on the pattern of this evolution. The analysis of the capsid proteins of Western equine encephalitis virus shows that the nucleocapsid of this virus can accommodate a considerable degree of variability in its protein component and that at least some regions of alphavirus capsid proteins show more extensive differences between different viruses than do the envelope glycoproteins.     

A large family of eukaryotic-like protein Ser/Thr kinases of Myxococcus xanthus, a developmental bacterium

Myxococcus xanthus is a gram-negative bacterium that forms multicellular fruiting bodies upon starvation. Here, we demonstrate that it contains at least 13 eukaryotic-like protein Ser/Thr kinases (Pkn1 to Pkn13) individually having unique features. All contain the kinase domain of approximately 280 residues near the N-terminal end, which share highly conserved features in eukaryotic Ser/Thr kinases. The kinase domain is followed by a putative regulatory domain consisting of 185 to 692 residues. These regulatory domains share no significant sequence similarities. The C-terminal regions of 11 kinases contain at least 1 transmembrane domain, suggesting that they function as transmembrane sensor kinases. From the recent genomic analysis, protein Ser/Thr kinases were found in various pathogenic bacteria and coexist with protein His kinases. Phylogenetic analysis of these Ser/Thr kinases reveals that all bacterial Ser/Thr kinases were evolved from a common ancestral kinase together with eukaryotic Tyr and Ser/Thr kinases. Coexistence of both Ser/Thr and His kinases in some organisms may be significant in terms of functional differences between the two kinases. We argue that both kinases are essential for some bacteria to adapt optimally to severe environmental changes.     

Characterization of Tannase Protein Sequences of Bacteria and Fungi: An In Silico Study

The tannase protein sequences of 149 bacteria and 36 fungi were retrieved from NCBI database. Among them only 77 bacterial and 31 fungal tannase sequences were taken which have different amino acid compositions. These sequences were analysed for different physical and chemical properties, superfamily search, multiple sequence alignment, phylogenetic tree construction and motif finding to find out the functional motif and the evolutionary relationship among them. The superfamily search for these tannase exposed the occurrence of proline iminopeptidase-like, biotin biosynthesis protein BioH, O-acetyltransferase, carboxylesterase/thioesterase 1, carbon–carbon bond hydrolase, haloperoxidase, prolyl oligopeptidase, C-terminal domain and mycobacterial antigens families and alpha/beta hydrolase superfamily. Some bacterial and fungal sequence showed similarity with different families individually. The multiple sequence alignment of these tannase protein sequences showed conserved regions at different stretches with maximum homology from amino acid residues 389–469 and 482–523 which could be used for designing degenerate primers or probes specific for tannase producing bacterial and fungal species. Phylogenetic tree showed two different clusters; one has only bacteria and another have both fungi and bacteria showing some relationship between these different genera. Although in second cluster near about all fungal species were found together in a corner which indicates the sequence level similarity among fungal genera. The distributions of fourteen motifs analysis revealed Motif 1 with a signature amino acid sequence of 29 amino acids, i.e. GCSTGGREALKQAQRWPHDYDGIIANNPA, was uniformly observed in 83.3 % of studied tannase sequences representing its participation with the structure and enzymatic function.     

Purification and properties of a double-stranded ribonuclease from the yeast Saccharomyces cerevisiae

The amino acid sequence of a single polypeptide chain, B-4, from fowl feather barbs has been determined. The B-4 chain was found to consist of 96 amino acid residues and to have a molecular weight of 10206 in the S-carboxymethylated form. The N terminus of this protein was an N-acetylserine residue. The B-4 protein contained seven S-carboxymethylcysteine residues, six of which are located in the N-terminal region (residues 1-26), and other one in C terminus. The central region of the peptide chain was rich in hydrophobic residues. There were homologous amino acids at 66 positions in the sequences of the feather keratins of fowl, emu and silver gull. The variation (substitution, deletion and insertion) in sequence was found to be localized in both terminal sections of the polypeptide chain. The B-4 protein structure was predicted to contain beta-sheet (about 30%), turn and random-coil-like structure, and no alpha-helix. beta-Sheet structure is mostly located in the central region (residues 22-70). On the other hand, both terminal regions are almost devoid of secondary structure.     

Multiple functional domains of Tat, the trans-activator of HIV-1, defined by mutational analysis.

The tat gene of HIV-1 is a potent trans-activator of gene expression from the HIV long terminal repeat (LTR). To define the functionally important regions of the product of the tat gene (Tat) of HIV-1, deletion, linker insertion and single amino acid substitution mutants within the Tat coding region of strain SF2 were constructed. The effect of these mutations on trans-activation was assessed by measuring the expression of the bacterial chloramphenicol acetyltransferase (CAT) reporter gene linked to the HIV-LTR. These studies have revealed that four different domains of the protein that map within the N-terminal 56 amino acid region are essential for Tat function. In addition to the essential domains, an auxiliary domain that enhances the activity of the essential region has also been mapped between amino acid residues 58 and 66. One of the essential domains maps in the N-terminal 20 amino acid region. The other three essential domains are highly conserved among the various strains of HIV-1 and HIV-2 as well as simian immunodeficiency virus (SIV). Of the conserved domains, one contains seven Cys residues and single amino acid substitutions for several Cys residues indicate that they are essential for Tat function. The second conserved domain contains a Lys X Leu Gly Ile X Tyr motif in which the Lys residue is essential for trans-activation and the other residues are partially essential. The third conserved domain is strongly basic and appears to play a dual role. Mutants lacking this domain are deficient in trans-activation and in efficient targeting of Tat to the nucleus and nucleolus. The combination of the four essential domains and the auxiliary domain contribute to the near full activity observed with the 101 amino acid Tat protein.     

Structural gene and complete amino acid sequence of Pseudomonas aeruginosa IFO 3455 elastase.

The DNA encoding the elastase of Pseudomonas aeruginosa IFO 3455 was cloned, and its complete nucleotide sequence was determined. When the cloned gene was ligated to pUC18, the Escherichia coli expression vector, bacteria carrying the gene exhibited high levels of both elastase activity and elastase antigens. The amino acid sequence, deduced from the nucleotide sequence, revealed that the mature elastase consisted of 301 amino acids with a relative molecular mass of 32,926 daltons. The amino acid composition predicted from the DNA sequence was quite similar to the chemically determined composition of purified elastase reported previously. We also observed nucleotide sequence encoding a signal peptide and "pro" sequence consisting of 197 amino acids upstream from the mature elastase protein gene. The amino acid sequence analysis revealed that both the N-terminal sequence of the purified elastase and the N-terminal side sequences of the C-terminal tryptic peptide as well as the internal lysyl peptide fragment were completely identical to the deduced amino acid sequences. The pattern of identity of amino acid sequences was quite evident in the regions that include structurally and functionally important residues of Bacillus subtilis thermolysin.     

Sequence complexity of histone H1 subtypes

H1 subtypes are involved in chromatin higher-order structure and gene regulation. H1 has a characteristic three-domain structure. We studied the length variation of the available H1 subtypes and showed that the length of the N-terminal and C-terminal domains was more variable than that of the central domain. The N-terminal and C-terminal domains were of low sequence complexity both at the nucleotide and at the amino acid level, whereas the globular domain was of high complexity. In most subtypes, low complexity was due only to cryptic simplicity, which reflects the clustering of a number of short and often imperfect sequence motifs. However, a subset of subtypes from eubacteria, plants, and invertebrates contained tandem repeats of short amino acid motifs (four to 12 residues), which could amount to a large proportion of the terminal domains. In addition, some other subtypes, such as those of Drosophila and mammalian H1t, were only marginally simple. The coexistence of these three kinds of subtypes suggests that the terminal domains could have originated in the amplification of short sequence motifs, which would then have evolved by point mutation and further slippage.     

Cluster analysis of amino acid indices for prediction of protein structure and function

The relationship among 222 published indices representing various physicochemical and biochemical properties of amino acid residues has been investigated by hierarchical cluster analysis. The clustering result is illustrated by the minimum spanning tree, which is conveniently divided into four regions: alpha and turn propensities, beta propensity, hydrophobicity and other physicochemical properties including, among others, bulkiness of amino acid residues. In addition, several subclasses of hydrophobicity scales have been identified: preference of inside and outside, accessible surface area, surrounding hydrophobicity and other mostly experimental scales including transfer free energy, partition coefficients, HPLC parameters and polarity. Representative amino acid indices are identified in each of these groups. The collection of amino acid indices is a useful resource for empirical analyses correlating sequence information with structural and functional properties of proteins. As an example, the indices that best reproduce the amino acid mutation data matrix are searched against this collection.     

Comparative In-Vitro Functional Analysis of Synthetic Defensins and Their Corresponding Peptide Variants Against HIV-1NL4.3, E. coli, S. aureus and P. aeruginosa

Defensins found in mammals belong to mainly two subfamilies α- and β-defensins. Mammalian defensins are small molecules (18–45 residues) that are cysteine, arginine rich compounds. Antimicrobial activities of these peptides were shown against a wide variety of microbes including bacteria, fungi, viruses and protozoan parasites. To investigate the structure and activity relationship, amino acid substitutions that alter charge were introduced into synthetic defensin peptides by adding 2–2 Arg (RR) and Asp (DD) at both the terminal and tested their effects on HIV-1, E. coli, S. aureus, and P. aeruginosa. In the present study, we have chemically synthesized native defensin peptides and their variants with Arg (RR) and Asp (DD) amino acid residues at N- and C-termini. Later, we assayed their anti-HIV, anti-microbial activities, stability, cytotoxicity and hemolytic properties. We reported that anti-HIV and antimicrobial activities of native defensins is increased significantly by adding Arg (RR) residues at both the termini while the substitution of Arg (RR) with Asp (DD), eliminate anti-HIV and antimicrobial activity against all bacterial species tested. While other physical features i.e. stability, cell toxicity and hemolytic property were not affected by any of the changes in the sequence. The results suggest that the terminal residues in defensins are crucial functional elements that determine their microbicidal potency. The enhanced microbicidal activity observed for defensin peptides with Arg (RR) residues could be due to optimization of amphiphilicity of the structure, which could facilitate specific interactions with the microbial membranes.     

Species-specific protein sequence and fold optimizations

Background

An organism's ability to adapt to its particular environmental niche is of fundamental importance to its survival and proliferation. In the largest study of its kind, we sought to identify and exploit the amino-acid signatures that make species-specific protein adaptation possible across 100 complete genomes.

Results

Environmental niche was determined to be a significant factor in variability from correspondence analysis using the amino acid composition of over 360,000 predicted open reading frames (ORFs) from 17 archae, 76 bacteria and 7 eukaryote complete genomes. Additionally, we found clusters of phylogenetically unrelated archae and bacteria that share similar environments by amino acid composition clustering. Composition analyses of conservative, domain-based homology modeling suggested an enrichment of small hydrophobic residues Ala, Gly, Val and charged residues Asp, Glu, His and Arg across all genomes. However, larger aromatic residues Phe, Trp and Tyr are reduced in folds, and these results were not affected by low complexity biases. We derived two simple log-odds scoring functions from ORFs (C_G) and folds (C_F) for each of the complete genomes. C_F achieved an average cross-validation success rate of 85 ± 8% whereas the C_G detected 73 ± 9% species-specific sequences when competing against all other non-redundant C_G. Continuously updated results are available at http://genome.mshri.on.ca.

Conclusion

Our analysis of amino acid compositions from the complete genomes provides stronger evidence for species-specific and environmental residue preferences in genomic sequences as well as in folds. Scoring functions derived from this work will be useful in future protein engineering experiments and possibly in identifying horizontal transfer events.