Biased Amino Acid Composition Research Articles

Evolution of protein lipograms: A bioinformatics problem

A protein lacking one of the 20 common amino acids is a protein lipogram. This open-ended problem-based learning assignment deals with the evolution of proteins with biased amino acid composition. It has students query protein and metabolic databases to test the hypothesis that natural selection has reduced the frequency of each amino acid specifically in the enzymes required for its biosynthesis. Student groups work in parallel on different amino acids and share strategies. Aside from content objectives that integrate knowledge of protein structure, function, synthesis, and evolution, this problem incorporates oral and written presentations, statistical analysis, and substantial decision making. The point of the problem described here is that a deficiency or absence of a particular amino acid in a protein may be more than a chance occurrence and may be driven by natural selection. The challenge is to demonstrate the difference.

The Population Biology of Mitochondrial DNA and Its Phylogenetic Implications

▪ Abstract The reconstruction of evolutionary trees from mitochondrial DNA (mtDNA) data is a common tool with which to infer the relationships of living organisms. The wide use of mtDNA stems from the ease of getting new sequence data for a set of orthologus genes and from the availability of many existing mtDNA sequences for a wide array of species. In this review we argue that developing a fuller understanding of the biology of mitochondria is essential for the rigorous application of mtDNA to inferences about the evolutionary history of species or populations. Though much progress has been made in understanding the parameters that shape the evolution of mitochondria and mtDNA, many questions still remain, and a better understanding of the role this organelle plays in regulating organismal fitness is becoming increasingly critical for accurate phylogeny reconstruction. In population biology, the limited information content of one nonrecombining genetic marker can compromise evolutionary inference, and the effects of nuclear genetic variation—and environmental factors—in mtDNA fitness differences can compound these problems. In systematics, the limited gene set, biased amino acid composition, and problems of compensatory substitutions can cloud phylogenetic signal. Dissecting the functional bases of these biases offers both challenges and opportunities in comparative biology.

Protein database searches using compositionally adjusted substitution matrices

Almost all protein database search methods use amino acid substitution matrices for scoring, optimizing, and assessing the statistical significance of sequence alignments. Much care and effort has therefore gone into constructing substitution matrices, and the quality of search results can depend strongly upon the choice of the proper matrix. A long-standing problem has been the comparison of sequences with biased amino acid compositions, for which standard substitution matrices are not optimal. To address this problem, we have recently developed a general procedure for transforming a standard matrix into one appropriate for the comparison of two sequences with arbitrary, and possibly differing compositions. Such adjusted matrices yield, on average, improved alignments and alignment scores when applied to the comparison of proteins with markedly biased compositions. Here we review the application of compositionally adjusted matrices and consider whether they may also be applied fruitfully to general purpose protein sequence database searches, in which related sequence pairs do not necessarily have strong compositional biases. Although it is not advisable to apply compositional adjustment indiscriminately, we describe several simple criteria under which invoking such adjustment is on average beneficial. In a typical database search, at least one of these criteria is satisfied by over half the related sequence pairs. Compositional substitution matrix adjustment is now available in NCBI's protein-protein version of blast.

Improvement of domain linker prediction by incorporating loop‐length‐dependent characteristics

Protein dissection into structural domains that can fold in isolation is an important issue in both functional and structural proteomics. Here, we analyzed inter- and intradomain loop sequences (respectively named domain linker and nonlinker loops) and computed a domain linker likelihood score, which was used for developing a domain boundary prediction protocol. The analysis confirmed our previous results indicating that the amino acid composition in terms of glycine, proline, aspartic acid, asparagine, lysine, and histidine significantly differs between linker and nonlinker loops. However, a detailed examination revealed that the amino acid composition bias actually depends on the loop length. Indeed, significant frequency deviations were observed for glycine, proline, and aspartic acid in short linker and nonlinker loops, whereas deviations were observed for aspartic acid, proline, asparagine, and lysine in long linker and nonlinker loops. Finally, we incorporated this loop-length-dependent amino acid composition bias in a simple linker prediction protocol, which predicted linkers with a 40.6% specificity and a 36.1% sensitivity. These figures are 4.4 and 2.4% higher than those obtained with our former prediction protocol that does not incorporate loop-length-dependent characteristics. This result should have practical significance for experimental protein dissection, since the probability of obtaining a stably folding structural domain by randomly dissecting a protein sequence is estimated to be 12.6%.

The Pairwise Energy Content Estimated from Amino Acid Composition Discriminates between Folded and Intrinsically Unstructured Proteins

The structural stability of a protein requires a large number of interresidue interactions. The energetic contribution of these can be approximated by low-resolution force fields extracted from known structures, based on observed amino acid pairing frequencies. The summation of such energies, however, cannot be carried out for proteins whose structure is not known or for intrinsically unstructured proteins. To overcome these limitations, we present a novel method for estimating the total pairwise interaction energy, based on a quadratic form in the amino acid composition of the protein. This approach is validated by the good correlation of the estimated and actual energies of proteins of known structure and by a clear separation of folded and disordered proteins in the energy space it defines. As the novel algorithm has not been trained on unstructured proteins, it substantiates the concept of protein disorder, i.e. that the inability to form a well-defined 3D structure is an intrinsic property of many proteins and protein domains. This property is encoded in their sequence, because their biased amino acid composition does not allow sufficient stabilizing interactions to form. By limiting the calculation to a predefined sequential neighborhood, the algorithm was turned into a position-specific scoring scheme that characterizes the tendency of a given amino acid to fall into an ordered or disordered region. This application we term IUPred and compare its performance with three generally accepted predictors, PONDR VL3H, DISOPRED2 and GlobPlot on a database of disordered proteins.

Bioinformatic identification of polymerizing and transmembrane mucins in the puffer fish Fugu rubripes.

Mucins are large glycoproteins characterized by mucin domains that show little sequence conservation and are rich in the amino acids Ser, Thr, and Pro. To effectively predict mucins from genomic and protein sequences obtained from genome projects, we developed a strategy based on the amino acid compositional bias characteristic of the mucin domains. This strategy is combined with an analysis of other features commonly found in mucins. Our method has now been used to predict mucins in the puffer fish Fugu rubripes that were previously not identified or annotated. At least three gel-forming mucins were found with the same general domain structure as the human MUC2 mucin. In addition one transmembrane mucin was identified with SEA and EGF domains as found in the mammalian transmembrane mucins. These results suggest that the number of gel-forming mucins has been conserved during evolution of the vertebrates, whereas the family of transmembrane mucins has been markedly expanded in the higher vertebrates.

Analysis of known bacterial protein vaccine antigens reveals biased physical properties and amino acid composition

Many vaccines have been developed from live attenuated forms of bacterial pathogens or from killed bacterial cells. However, an increased awareness of the potential for transient side-effects following vaccination has prompted an increased emphasis on the use of sub-unit vaccines, rather than those based on whole bacterial cells. The identification of vaccine sub-units is often a lengthy process and bioinformatics approaches have recently been used to identify candidate protein vaccine antigens. Such methods ultimately offer the promise of a more rapid advance towards preclinical studies with vaccines. We have compared the properties of known bacterial vaccine antigens against randomly selected proteins and identified differences in the make-up of these two groups. A computer algorithm that exploits these differences allows the identification of potential vaccine antigen candidates from pathogenic bacteria on the basis of their amino acid composition, a property inherently associated with sub-cellular location.

Probabilistic analysis indicates discordant gene trees in chloroplast evolution.

Analyses of whole-genome data often reveal that some genes have evolutionary histories that diverge from the majority phylogeny estimated for the entire genome. We present a probabilistic model that deals with heterogeneity among gene trees, implement it via the Gibbs sampler, and apply it to the plastid genome. Plastids and their genomes are transmitted as a single block without recombination, hence homogeneity among gene trees within this genome is expected. Nevertheless, previous work has revealed clear heterogeneity among plastid genes (e.g., Delwiche and Palmer 1996). Other studies, using whole plastid genomes of various algae and land plants, found little additional heterogeneity (Martin et al. 1998; Adachi et al. 2000). We augment the earlier studies by using a data set of 14 taxa: 6 land plants, 2 green algae, a diatom, 2 red algae and a cryptophyte, the cyanelle of the glaucocystophyte Cyanophora, and the blue-green alga Synechocystis as an outgroup. Contrary to the earlier analyses, we cannot find even a single, dominant consensus tree. Therefore, we formulate a probabilistic model that divides the genes into two sets: those that follow the consensus tree and those that have independent gene trees. No particular tree is supported by more than three-fourths of the genes. But the set of genes that follows a certain tree is fairly independent of data processing and the method of analysis. With one possible exception, we find no evidence for collinear or functionally related genes to follow similar trees. The phylogenetic pattern also seems independent of bias in amino acid composition. Among possible explanations for the observed phenomenon, the hypothesis that different genes have different covarion structures is difficult to assess. But gene duplication may be possible through the inverted or direct repeat regions, while horizontal gene transfer seems less likely. In contrast to green algae and land plants, inverted repeat regions in red algae and in Cyanophora show abundant differences among the copies. Thus, genes may get duplicated when they are recruited into the inverted repeat region and one of the two copies may be lost after leaving the inverted repeat region.

Unique Amino Acid Composition of Proteins in Halophilic Bacteria

The amino acid compositions of proteins from halophilic archaea were compared with those from non-halophilic mesophiles and thermophiles, in terms of the protein surface and interior, on a genome-wide scale. As we previously reported for proteins from thermophiles, a biased amino acid composition also exists in halophiles, in which an abundance of acidic residues was found on the protein surface as compared to the interior. This general feature did not seem to depend on the individual protein structures, but was applicable to all proteins encoded within the entire genome. Unique protein surface compositions are common in both halophiles and thermophiles. Statistical tests have shown that significant surface compositional differences exist among halophiles, non-halophiles, and thermophiles, while the interior composition within each of the three types of organisms does not significantly differ. Although thermophilic proteins have an almost equal abundance of both acidic and basic residues, a large excess of acidic residues in halophilic proteins seems to be compensated by fewer basic residues. Aspartic acid, lysine, asparagine, alanine, and threonine significantly contributed to the compositional differences of halophiles from meso- and thermophiles. Among them, however, only aspartic acid deviated largely from the expected amount estimated from the dinucleotide composition of the genomic DNA sequence of the halophile, which has an extremely high G+C content (68%). Thus, the other residues with large deviations (Lys, Ala, etc.) from their non-halophilic frequencies could have arisen merely as “dragging effects” caused by the compositional shift of the DNA, which would have changed to increase principally the fraction of aspartic acid alone.

Evolutionary analysis ofArabidopsis, cyanobacterial, and chloroplast genomes reveals plastid phylogeny and thousands of cyanobacterial genes in the nucleus

Chloroplasts were once free-living cyanobacteria that became endosymbionts, but the genomes of contemporary plastids encode only approximately 5-10% as many genes as those of their free-living cousins, indicating that many genes were either lost from plastids or transferred to the nucleus during the course of plant evolution. Previous estimates have suggested that between 800 and perhaps as many as 2,000 genes in the Arabidopsis genome might come from cyanobacteria, but genome-wide phylogenetic surveys that could provide direct estimates of this number are lacking. We compared 24,990 proteins encoded in the Arabidopsis genome to the proteins from three cyanobacterial genomes, 16 other prokaryotic reference genomes, and yeast. Of 9,368 Arabidopsis proteins sufficiently conserved for primary sequence comparison, 866 detected homologues only among cyanobacteria and 834 other branched with cyanobacterial homologues in phylogenetic trees. Extrapolating from these conserved proteins to the whole genome, the data suggest that approximately 4,500 of Arabidopsis protein-coding genes ( approximately 18% of the total) were acquired from the cyanobacterial ancestor of plastids. These proteins encompass all functional classes, and the majority of them are targeted to cell compartments other than the chloroplast. Analysis of 15 sequenced chloroplast genomes revealed 117 nuclear-encoded proteins that are also still present in at least one chloroplast genome. A phylogeny of chloroplast genomes inferred from 41 proteins and 8,303 amino acids sites indicates that at least two independent secondary endosymbiotic events have occurred involving red algae and that amino acid composition bias in chloroplast proteins strongly affects plastid genome phylogeny.

Protein surface amino acid compositions distinctively differ between thermophilic and mesophilic bacteria

One of the well-known observations of proteins from thermophilic bacteria is the bias of the amino acid composition in which charged residues are present in large numbers, and polar residues are scarce. On the other hand, it has been reported that the molecular surfaces of proteins are adapted to their subcellular locations, in terms of the amino acid composition. Thus, it would be reasonable to expect that the differences in the amino acid compositions between proteins of thermophilic and mesophilic bacteria would be much greater on the protein surface than in the interior. We performed systematic comparisons between proteins from thermophilic bacteria and mesophilic bacteria, in terms of the amino acid composition of the protein surface and the interior, as well as the entire amino acid chains, by using sequence information from the genome projects. The biased amino acid composition of thermophilic proteins was confirmed, and the differences from those of mesophilic proteins were most obvious in the compositions of the protein surface. In contrast to the surface composition, the interior composition was not distinctive between the thermophilic and mesophilic proteins. The frequency of the amino acid pairs that are closely located in the space was also analyzed to show the same trend of the single amino acid compositions. Interestingly, extracellular proteins from mesophilic bacteria showed an inverse trend against thermophilic proteins (i.e. a reduced number of charged residues and rich in polar residues). Nuclear proteins from eukaryotes, which are known to be abundant in positive charges, showed different compositions as a whole from the thermophiles. These results suggest that the bias of the amino acid composition of thermophilic proteins is due to the residues on the protein surfaces, which may be constrained by the extreme environment.

Low-Complexity Regions in Plasmodium falciparum Proteins

Full-sequence data available for Plasmodium falciparumchromosomes 2 and 3 are exploited to perform a statistical analysis of the long tracts of biased amino acid composition that characterize the vast majority of P. falciparum proteins and to make a comparison with similarly defined tracts from other simple eukaryotes. When the relatively minor subset of prevalently hydrophobic segments is discarded from the set of low-complexity segments identified by current segmentation methods in P. falciparum proteins, a good correspondence is found between prevalently hydrophilic low-complexity segments and the species-specific, rapidly diverging insertions detected by multiple-alignment procedures when sequences of bona fide homologs are available. Amino acid preferences are fairly uniform in the set of hydrophilic low-complexity segments identified in the twoP. falciparum chromosomes sequenced, as well as in sequenced genes from Plasmodium berghei, but differ from those observed in Saccharomyces cerevisiae and Dictyostelium discoideum. In the two plasmodial species, amino acid frequencies do not correlate with properties such as hydrophilicity, small volume, or flexibility, which might be expected to characterize residues involved in nonglobular domains but do correlate with A-richness in codons. An effect of phenotypic selection versus neutral drift, however, is suggested by the predominance of asparagine over lysine.

Low-complexity regions in Plasmodium falciparum proteins.

Full-sequence data available for Plasmodium falciparum chromosomes 2 and 3 are exploited to perform a statistical analysis of the long tracts of biased amino acid composition that characterize the vast majority of P. falciparum proteins and to make a comparison with similarly defined tracts from other simple eukaryotes. When the relatively minor subset of prevalently hydrophobic segments is discarded from the set of low-complexity segments identified by current segmentation methods in P. falciparum proteins, a good correspondence is found between prevalently hydrophilic low-complexity segments and the species-specific, rapidly diverging insertions detected by multiple-alignment procedures when sequences of bona fide homologs are available. Amino acid preferences are fairly uniform in the set of hydrophilic low-complexity segments identified in the two P. falciparum chromosomes sequenced, as well as in sequenced genes from Plasmodium berghei, but differ from those observed in Saccharomyces cerevisiae and Dictyostelium discoideum. In the two plasmodial species, amino acid frequencies do not correlate with properties such as hydrophilicity, small volume, or flexibility, which might be expected to characterize residues involved in nonglobular domains but do correlate with A-richness in codons. An effect of phenotypic selection versus neutral drift, however, is suggested by the predominance of asparagine over lysine.

Issues in predicting protein function from sequence.

Identifying homologues, defined as genes that arose from a common evolutionary ancestor, is often a relatively straightforward task, thanks to recent advances made in estimating the statistical significance of sequence similarities found from database searches. The extent by which homologues possess similarities in function, however, is less amenable to statistical analysis. Consequently, predicting function by homology is a qualitative, rather than quantitative, process and requires particular care to be taken. This review focuses on the various approaches that have been developed to predict function from the scale of the atom to that of the organism. Similarities in homologues' functions differ considerably at each of these different scales and also vary for different domain families. It is argued that due attention should be paid to all available clues to function, including orthologue identification, conservation of particular residue types, and the co-occurrence of domains in proteins. Pitfalls in database searching methods arising from amino acid compositional bias and database size effects are also discussed.

Gum arabic glycoprotein contains glycomodules of both extensin and arabinogalactan-glycoproteins

Gum arabic glycoprotein (GAGP) is a large molecular weight, hydroxyproline-rich arabinogalactan-protein (AGP) component of gum arabic. GAGP has a simple, highly biased amino acid composition indicating a repetitive polypeptide backbone. Previous work (Qi, W., Fong, C., Lamport, D.T.A., 1991. Plant Physiology 96, 848), suggested small (∼11 residue) repetitive peptide motifs each with three Hyp-arabinoside attachment sites and a single Hyp-arabinogalactan polysaccharide attachment site. We tested that hypothesis by sequence analysis of the GAGP polypeptide after HF-deglycosylation. A family of closely related peptides confirmed the presence of a repetitive 19-residue consensus motif. However, the motif: Ser-Hyp-Hyp-Hyp-Thr-Leu-Ser-Hyp-Ser-Hyp-Thr-Hyp-Thr-Hyp-Hyp-Leu-Gly-Pro-His, was about twice the size anticipated. Thus, judging by Hyp-glycoside profiles of GAGP, the consensus motif contained six Hyp-arabinosides rather than three and two Hyp-polysaccharides rather than one. We inferred the glycosylation sites based on the Hyp contiguity hypothesis which predicts arabinosides on contiguous Hyp residues and arabinogalactan polysaccharides on clustered non-contiguous Hyp residues, i.e. the GAGP motif would consist of arabinosylated contiguous Hyp blocks flanking two central Hyp-polysaccharides. We predict this rigidifies the glycoprotein, enhances the overall symmetry of the glycopeptide motif, and may explain some of the remarkable properties of gum arabic.

A Method for Predicting Transcription Factors Using Local Regions of Biased Amino Acid Composition

A Method for Predicting Transcription Factors Using Local Regions of Biased Amino Acid Composition

Structure and expression of the highly repetitive histone H1-related sperm chromatin proteins from winter flounder.

In the late stages of spermatogenesis, winter flounder produce a family of high molecular mass (80-200 kDa) basic nuclear proteins (HMrBNPs) that combine with the normal complement of histones to produce condensed sperm chromatin with an increased nucleosomal repeat length. The HMrBNPs have a biased amino-acid composition in which Arg, Ser, Lys and Pro are abundant because of their presence in many simple peptide repeats. The organization of these repeats was deduced by cDNA cloning. The predominant repeating units are related 26- and 30-amino-acid sequences that in turn are linked by 6-amino-acid spacers to form 58- and 62-amino-acid repeats. Subsets of these repeats are also present, such as a dispersed 20-amino-acid repeat and a tandem array of nine heptapeptides at the C-terminus. The HMrBNPs appear to have evolved from an extreme H1 variant that has an N-terminal tail of HMrBNP-like sequence linked to an H1 globular region. Based on sequences of the most abundant HMrBNP cDNAs, and the lack of hybridization between HMrBNP mRNAs and a DNA probe for the H1 globular region, the latter domain appears to have been lost during expansion and amplification of the HMrBNP-like repeats. Transcripts of the HMrBNP and H1 variant genes are present in testis RNAs only during the mid-spermatid stage of spermatogenesis, at the same time that HMrBNPs in their highly phosphorylated form first appear in the nucleus. Judging by the lack of a lag between HMrBNP mRNA synthesis and translation, the mRNAs for these highly basic proteins are not stored for any length of time. Instead, the deposition of HMrBNPs onto DNA, which coincides with the major reorganization and silencing of the chromatin, may be controlled by dephosphorylation.

The wheat LEA protein Em functions as an osmoprotective molecule in Saccharomyces cerevisiae.

The biased amino acid composition and aperiodic (random coil) configuration of Group 1 late embryogenesis-abundant (LEA) proteins imply that these proteins are capable of binding large amounts of water. While Group 1 LEAs have been predicted to contribute to osmotic stress protection in both embryonic and vegetative tissues, biochemical support has been lacking. We have used Saccharomyces cerevisiae as a model system to test the putative osmoprotective function of a wheat Group 1 LEA protein, Em. We demonstrate that expression of Em protein in yeast cells is not deleterious to growth in media of normal osmolarity and attenuates the growth inhibition normally observed in media of high osmolarity. Enhanced growth is observed in the presence of a variety of osmotically active compounds indicating that Em protein is capable of mitigating the detrimental effect of low water potential in a relatively non-specific manner. These results are the first biochemical demonstration of an osmoprotective function for a Group 1 LEA and suggest that the yeast expression system will be useful in dissecting the mechanism of protection through structure-function studies.

SIRE-1, a copia/Ty1-like retroelement from soybean, encodes a retroviral envelope-like protein.

The soybean genome hosts a family of several hundred, relatively homogeneous copies of a large, copia/Ty1-like retroelement designated SIRE-1. A copy of this element has been recovered from a Glycine max genomic library. DNA sequence analysis of two SIRE-1 subclones revealed that SIRE-1 contains a long, uninterrupted, ORF between the 3' end of the pol ORF and the 3' long terminal repeat (LTR), a region that harbors the env gene in retroviral genomes. Conceptual translation of this second ORF produces a 70-kDa protein. Computer analyses of the amino acid sequence predicted patterns of transmembrane domains, alpha-helices, and coiled coils strikingly similar to those found in mammalian retroviral envelope proteins. In addition, a 65-residue, proline-rich domain is characterized by a strong amino acid compositional bias virtually identical to that of the 60-amino acid, proline-rich neutralization domain of the feline leukemia virus surface protein. The assignment of SIRE-1 to the copia/Ty1 family was confirmed by comparison of the conceptual translation of its reverse transcriptase-like domain with those of other retroelements. This finding suggests the presence of a proretrovirus in a plant genome and is the strongest evidence to date for the existence of a retrovirus-like genome closely related to copia/Ty1 retrotransposons.

Synonymous nucleotide divergence: what is "saturation"?

The nucleotide divergence at synonymous third sites between two lineages will increase with time since the latest common ancestor, up to some saturation level. The "null-hypothesis divergence" is defined as the percentage of difference predicted at synonymous third sites, allowing for amino acid composition and codon bias, but assuming that codon bias is the same at all sites occupied by a given amino acid, when equilibrium has been reached between forward and backward substitutions. For two highly expressed genes, gapA and ompA, in the enterobacteria, the estimated values of the null-hypothesis divergence are 39.3 and 38.15%, respectively, compared to estimated values of saturation divergence of 19.0 and 25.4%. A possible explanation for this discrepancy is that different codons for a given amino acid are favored at different sites in the same gene.