Selection on Synthesis Cost Affects Interprotein Amino Acid Usage in All Three Domains of Life

Most investigations of the forces shaping protein evolution have focussed on protein function. However, cells are typically 50%–75% protein by dry weight, with protein expression levels distributed over five orders of magnitude. Cells may, therefore, be under considerable selection pressure to incorporate amino acids that are cheap to synthesize into proteins that are highly expressed. Such selection pressure has been demonstrated to alter amino acid usage in a few organisms, but whether “cost selection” is a general phenomenon remains unknown. One reason for this is that reliable protein expression level data is not available for most organisms. Accordingly, I have developed a new method for detecting cost selection. This method depends solely on interprotein gradients in amino acid usage. Applying it to an analysis of 43 whole genomes from all three domains of life, I show that selection on the synthesis cost of amino acids is a pervasive force in shaping the composition of proteins. Moreover, some amino acids have different price tags for different organisms—the cost of amino acids is changed for organisms living in hydrothermal vents compared with those living at the sea surface or for organisms that have difficulty acquiring elements such as nitrogen compared with those that do not—so I also investigated whether differences between organisms in amino acid usage might reflect differences in synthesis or acquisition costs. The results suggest that organisms evolve to alter amino acid usage in response to environmental conditions. Electronic supplementary material The online version of this article (doi:) contains supplementary material, which is available to authorized users. [Reviewing Editor: Hector Musto]     

A search for structural factors responsible for the stability of proteins from thermophilic organisms

A database was designed to include 392 pairs of homologous proteins from thermophilic and mesophilic organisms. Proteins from thermophilic organisms proved to contain more atom-atom contacts per residue as compared with their mesophilic homologs. Solvent-accessible exterior amino acid residues contribute to the increase in the number of contacts. The amino acid composition was analyzed for internal (solvent-inaccessible) and exterior amino acid residues of thermophilic and mesophilic proteins. The exterior residues of thermophils have higher contents of Lys, Arg, and Glu and lower contents of Ala, Asp, Asn, Gln, Ser, and Thr as compared with mesophilic proteins. Interior protein regions did not differ in amino acid composition.     

Similarities in protein amino acid composition in connection with principles of protein evolution

The amino acid composition of human alcohol dehydrogenase (ADH) was compared with alcohol dehydrogenases from different organisms and with other proteins. Similar amino acid sequences in human ADH (template protein) and in other proteins were determined by means of an original computer program. Analysis of amino acid motifs reveals that the ADHs from evolutionary more close organisms have more common amino acid sequences. The quantity measure of amino acid similarity was the number of similar motifs in analyzed protein per protein length. This value was measured for ADHs and for different proteins. For ADHs, this quotient was higher than for proteins with different functions; for vertebrates it correlated with evolutionary closeness. The similar operation of motif comparison was made with the help of program complex “MEME”. The analysis of ADHs revealed 4 motifs common to 6 of 10 tested organisms and no such motifs for proteins of different function. The conclusion is that general amino composition is more important for protein function than amino acid order and for enzymes of similar function it better correlates with evolutionary distance between organisms.     

Beta-subunit of ATP-synthase: a useful marker for studying the phylogenetic relationship of eubacteria

The genes encoding the beta-subunits of ATP-synthases (ATPases) from Bacteroides fragilis DSM 2151, Cytophaga lytica DSM 2039 and 'Taxeobacter ocellatus' were cloned. The nucleotide sequences were determined completely for the genes of the first two organisms and to a major part for that of 'T. ocellatus'. The predicted amino acid sequences were compared with previously published amino acid sequences of beta-subunits. Two characteristic insertions were found in genes from organisms belonging to the so-called bacteroides-cytophaga-flavo-bacterium group. The remaining structure shows a high degree of sequence similarity within this group. These data support the conclusions drawn from comparative 16S rRNA sequence analyses that organisms in this phenotypically heterogeneous group are phylogenetically related. A phylogenetic tree was constructed based on a distance matrix of optimally aligned amino acid sequences of beta-subunits of ATPases of various eubacteria, chloroplasts and mitochondria. It is in good agreement with a tree derived from 16S rRNA sequence analyses.     

Analysis of differences in amino acid substitution patterns, using multilevel G-tests

In this paper, a new algorithm is presented, which makes possible multilevel comparison of BLOSUM protein substitution matrices based on data from different groups of organisms. As an example, a comparison between substitution matrices based on data from two groups of bacterial genomes with different GC content is presented. Our approach includes evaluating the number of amino acid pairs in BLOCKS databases created separately for the two groups of bacteria using protein sequences deposited in the COG database. Differences of distributions of amino acid pair counts are tested using the chi-squared based G-test. Different analysis levels make it possible to distinguish different patterns of amino acid substitution. Application of the algorithm reveals statistically significant differences in amino acid substitution patterns between AT-rich and GC-rich groups of bacterial organisms. The differences are particularly visible in the overall substitution pattern, amino acid conservation pattern and in comparison of substitution patterns for single amino acids. The algorithm presented in this paper can be considered a novel method for multi-level comparison of amino acid substitution patterns. The presented approach is not limited to bacterial organisms and BLOSUM substitution matrices. Statistically significant differences between substitution patterns in the two groups of bacterial organisms with respect to amino acid conservation pattern can be the evidence of different rate of evolutionary change between AT-rich and GC-rich bacterial organisms.     

Search for structural factors responsible for the stability of proteins from thermophilic organisms

Database including 392 homologous pairs of proteins from thermophilic and mesophilic organisms was created. Using this database we have found that proteins from termophilic organisms contain more atom-atom contacts per residue in comparison with mesophilic homologues. Contribution to increase of the number of contacts gives exterior amino acid residues, accessible for the solvent. Amino acid composition of interior, inaccessible for the solvent, and exterior amino acid residues of proteins from thermophilic and mesophilic organisms were analyzed. We have obtained that exterior residues of proteins from thermophilic organisms contain more such amino acid residues as Lys, Arg and Glu and smaller such amino acid residues as Ala, Asp, Asn. Gln, Ser, and Thr in comparison with proteins from mesophilic organisms. Amino acid compositions of interior residues of considered proteins are not different.     

Relative Amino Acid Composition Signatures of Organisms and Environments

Background

Identifying organism-environment interactions at the molecular level is crucial to understanding how organisms adapt to and change the chemical and molecular landscape of their habitats. In this work we investigated whether relative amino acid compositions could be used as a molecular signature of an environment and whether such a signature could also be observed at the level of the cellular amino acid composition of the microorganisms that inhabit that environment.

Methodologies/Principal Findings

To address these questions we collected and analyzed environmental amino acid determinations from the literature, and estimated from complete genomic sequences the global relative amino acid abundances of organisms that are cognate to the different types of environment. Environmental relative amino acid abundances clustered into broad groups (ocean waters, host-associated environments, grass land environments, sandy soils and sediments, and forest soils), indicating the presence of amino acid signatures specific for each environment. These signatures correlate to those found in organisms. Nevertheless, relative amino acid abundance of organisms was more influenced by GC content than habitat or phylogeny.

Conclusions

Our results suggest that relative amino acid composition can be used as a signature of an environment. In addition, we observed that the relative amino acid composition of organisms is not highly determined by environment, reinforcing previous studies that find GC content to be the major factor correlating to amino acid composition in living organisms.     

Cold-adapted tubulins in the glacier ice worm, Mesenchytraeus solifugus

Glacier ice worms, Mesenchytraeus solifugus and related species, are the only known annelids that survive obligately in glacier ice and snow. One fundamental component of cold temperature adaptation is the ability to polymerize tubulin, which typically depolymerizes at low physiological temperatures (e.g., <10 degrees C) in most temperate species. In this study, we isolated two alpha-tubulin (Msalpha1, Msalpha2) and two beta-tubulin (Msbeta1, Msbeta2) subunits from an ice worm cDNA library, and compared their predicted amino acid sequences with homologues from other cold-adapted organisms (e.g., Antarctic fish, ciliate) in an effort to identify species-specific amino acid substitutions that contribute to cold temperature-dependent tubulin polymerization. Our comparisons and predicted protein structures suggest that ice worm-specific amino acid substitutions stabilize lateral contact associations, particularly between beta-tubulin protofilaments, but these substitutions occur at different positions in comparison with other cold-adapted tubulins. The ice worm tubulin gene family appears relatively small, comprising one primary alpha- and one primary beta-tubulin monomers, though minor isoforms and psuedogenes were identified. Our analyses suggest that variation occurs in the strategies (i.e., species-specific amino acid substitutions, gene number) by which cold-adapted taxa have evolved the ability to polymerize tubulin at low physiological temperatures.     

Using First Passage Statistics to Extract Environmentally Dependent Amino Acid Correlations

In this work, we study the first passage statistics of amino acid primary sequences, that is the probability of observing an amino acid for the first time at a certain number of residues away from a fixed amino acid. By using this rich mathematical framework, we are able to capture the background distribution for an organism, and infer lengths at which the first passage has a probability that differs from what is expected. While many features of an organism''s genome are due to natural selection, others are related to amino acid chemistry and the environment in which an organism lives, constraining the randomness of genomes upon which selection can further act. We therefore use this approach to infer amino acid correlations, and then study how these correlations vary across a wide range of organisms under a wide range of optimal growth temperatures. We find a nearly universal exponential background distribution, consistent with the idea that most amino acids are globally uncorrelated from other amino acids in genomes. When we are able to extract significant correlations, these correlations are reliably dependent on optimal growth temperature, across phylogenetic boundaries. Some of the correlations we extract, such as the enhanced probability of finding, for the first time, a cysteine three residues away from a cysteine or glutamic acid two residues away from an arginine, likely relate to thermal stability. However, other correlations, likely appearing on alpha helical surfaces, have a less clear physiochemical interpretation and may relate to thermal stability or unusual metabolic properties of organisms that live in a high temperature environment.     

Production of amino acids by free-living heterotrophic nitrogen-fixing bacteria

Summary Interest of microbial production of amino acids has been increased greatly since development of biotechnological methods. These methods represent a perspective way applied in a future large-scale manufacture of inexpensive amino acids. In this context, the isolation of producing organisms that may be exploited in the desing of alternative methods for the production of amino acids could be of primary importance.In this review we will describe the liberation of amino acids (methionine, lysine, arginine, tryptophane and glutamic acid) byAzotobacter andAzospirillum during growth in culture media with different carbon sources under diazotrophic and adiazotrophic conditions. These organisms may be useful in developing new methods for the industrial production of amino acids.     

Amino acid cost and codon-usage biases in 6 prokaryotic genomes: a whole-genome analysis

For most prokaryotic organisms, amino acid biosynthesis represents a significant portion of their overall energy budget. The difference in the cost of synthesis between amino acids can be striking, differing by as much as 7-fold. Two prokaryotic organisms, Escherichia coli and Bacillus subtilis, have been shown to preferentially utilize less costly amino acids in highly expressed genes, indicating that parsimony in amino acid selection may confer a selective advantage for prokaryotes. This study confirms those findings and extends them to 4 additional prokaryotic organisms: Chlamydia trachomatis, Chlamydophila pneumoniae AR39, Synechocystis sp. PCC 6803, and Thermus thermophilus HB27. Adherence to codon-usage biases for each of these 6 organisms is inversely correlated with a coding region's average amino acid biosynthetic cost in a fashion that is independent of chemoheterotrophic, photoautotrophic, or thermophilic lifestyle. The obligate parasites C. trachomatis and C. pneumoniae AR39 are incapable of synthesizing many of the 20 common amino acids. Removing auxotrophic amino acids from consideration in these organisms does not alter the overall trend of preferential use of energetically inexpensive amino acids in highly expressed genes.     

A concept of amino acid archaeorelation: Origin of life and the genetic code

Using certain assumptions, an attempt has been made to obtain a quantitative evaluation of the similarity of the bulk amino acid compositions of various organisms and of the quantitative ratios of amino acids in the products of the presumable prebiological ructions and in the natural proteins. The results of comparison between the molar ratios of amino acids in synthesized products and in natural proteins and the quantitative distribution of triplets among amino acids in the existing genetic code are given. Proceeding from the correlations found and the data reported in literature, a concept has been suggested amounting, in essence, to the assertion that the structure of the genetic code (the quantitative distribution of the triplets) is in a certain way determined by the pre-existing, primary ratio (archaeorelation) of amino acids.     

Cloning and nucleotide sequencing of genes for small, acid-soluble spore proteins of Bacillus cereus, Bacillus stearothermophilus, and "Thermoactinomyces thalpophilus"

As found previously with other Bacillus species, spores of B. stearothermophilus and "Thermoactinomyces thalpophilus" contained significant levels of small, acid-soluble spore proteins (SASP) which were rapidly degraded during spore germination and which reacted with antibodies raised against B. megaterium SASP. Genes coding for a B. stearothermophilus and a "T. thalpophilus" SASP as well as for two B. cereus SASP were cloned, their nucleotide sequences were determined, and the amino acid sequences of the SASP coded for were compared. Strikingly, all of the amino acid residues previously found to be conserved in this group of SASP both within and between two other Bacillus species (B. megaterium and B. subtilis) were also conserved in the SASP coded for by the B. cereus genes as well as those coded for by the genes from the more distantly related organisms B. stearothermophilus and "T. thalpophilus." This finding strongly suggests that there is significant selective pressure to conserve SASP primary sequence and thus that these proteins serve some function other than simply amino acid storage.     

Biosynthesis and control of beta-lactam antibiotics: the early steps in the "classical" tripeptide pathway

The interaction between growth and secondary metabolism develops from physiological responses of the producer organism to its environment. Nutrients are channelled into primary growth processes or into secondary processes such as antibiotic biosynthesis by a variety of metabolic controls, the nature of which has been extensively studied in organisms producing beta-lactam antibiotics via the tripeptide, delta-(L-alpha-aminoadipyl)-L-cysteinyl-D-valine. In the following article we review the early stages of beta-lactam biosynthesis in fungi and actinomycetes, keeping in mind the regulation of primary pathways that provide the amino acid precursors of this group of antibiotics, as well as the regulation of the secondary pathway itself. Of special importance to organisms engaging in secondary metabolism are the control mechanisms that suppress the nonessential process during rapid growth but allow secondary metabolic genes to be expressed and resources to be diverted when environmental factors generate the appropriate biochemical signals.     

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. Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

Porcine thyrotropin. The amino-acid sequence of the alpha and beta subunits.

The amino acid sequences of both the alpha and beta subunits of porcine thyrotropin have been studied. Bovine thyrotropin primary structure was taken as a model for ordering the tryptic peptides of porcine thyrotropin. The amino acid sequence of the alpha subunit is identical to that of porcine luteinizing hormone, while oligosaccharide side-chains differ in composition. The primary structure of the beta subunit differs from that of bovine thyrotropin by six amino acid replacements, in positions 22, 24, 26, 36, 62 and 69, and by the absence of a methionyl residue at the carboxy terminus. Chemical evolutions of thyrotropin and luteinizing hormone are compared.     

The amino acid sequence of cytochrome c553 from Microcystis aeruginosa

Cytochrome c553 is an electron donor to P700 in the photosynthetic electron transfer chain of cyanobacteria and eukaryotic algae. We have purified this cytochrome from the cyanobacterium Microcystis aeruginosa and determined its amino acid sequence. When the amino acid sequence of this protein is compared to sequences of cytochromes c553 from other organisms, one sees that the evolution of net charge is more pronounced than the evolution of overall structure, further documenting a pronounced shift in the isoelectric point of this protein during the evolution of cyanobacteria. Cyanobacteria and algae also contain cytochrome c550 (Mr 15,500) which is quite different from cytochrome c553 (Mr 10,500). When the amino acid sequence of cytochrome c553 is compared to that of cytochrome c550, two regions of similar sequence are recognized.     

Post-translational isoprenylation of tryptophan

Bacillus subtilis and related bacilli produce a post-translationally modified oligopeptide, ComX pheromone, that stimulates natural genetic competence controlled by quorum sensing. The ComX pheromones are formed by geranylation or farnesylation on a tryptophan residue at the 3 position of its indole ring. This results in the formation of a tricyclic structure including, a newly formed five-membered ring, similar to proline. Isoprenylation of ComX to form ComX pheromones is essential for pheromonal activity, and is functionally more crucial than its amino acid sequence. The ComX pheromone is the first example of isoprenoidal modifiations of tryptophan residues in living organisms and post-translational isoprenylation of any amino acid in prokaryotes. Because the presence of geranylated compounds is unusual in primary and secondary metabolites outside the plant kingdom, post-translational geranylation in bacilli is unprecedented in nature.     

Optimal growth temperature of prokaryotes correlates with class II amino acid composition

Partitioning of aminoacyl-tRNA synthetases and their associated amino acids into two classes allows us to distinguish between thermophilic and mesophilic species based only on amino acids composition. The CLASSDB program has been developed for amino acid content analysis in organisms treated individually or pooled together to form a pattern of characteristic properties. A strong correlation has been observed between optimal growth temperature (OGT) of organisms and class II amino acids content. Amino acid composition in organisms closely related phylogenetically but dissimilar in their OGT testifies that thermo-adaptation happens rather rapidly on the time scale of evolution.