Preferred codons and amino acid couples in hyperthermophiles

Background

Most organisms grow at temperatures from 20 to 50°C but some prokaryotes, including Archaea and Bacteria, are capable of withstanding higher temperatures, from 60 to >100°C. What makes these cells so resistant to heat? Their biomolecules must be sufficiently stable, especially proteins, to work under these extreme conditions, but the bases for thermostability remains elusive.

Results

The preferential usage of certain couples of amino acids and codons in thermal adaptation was investigated, by comparative proteome analysis, using 28 complete genomes from 18 mesophiles, 4 thermophiles, and 6 hyperthermophiles. In the hyperthermophiles proteomes, whenever the percent of Glu (E) and Lys (K) Increased, the percent of Gln (Q) and His (H) decreased, so that the E+K/Q+H ratio was > 4,5; in the mesophiles proteomes, it was < 2,5 and in the thermophiles an intermediary value was observed. The E+K/Q+H ratios for chaperonins, potentially thermostable proteins, were higher than their proteome ratios whereas, for DNA ligases, not necessarily thermostable, they followed the proteome ones. Analysis of codon usage revealed that hyperthermophiles preferred AGR codons for Arg in detriment of CGN codons, which were preferred by mesophiles.

Conclusions

The results suggested that the E+K/Q+H ratio may provide a useful mark for distinguishing hyperthermophilic, thermophilic and mesophilic prokaryotes and that the high percent of the amino acid couple E+K, consistently associated to the low percent of the pair Q+H, could contribute to protein thermostability. Second, the preference for AGR codons for Arg was a signature of all hyperthermophilics so far analyzed.

     

Thermo-search: lifestyle and thermostability analysis

Thermo-search is an online web tool for the analysis of proteomes and individual proteins according to the ratio of two couplets of preferred and avoided amino acids in hyperthermophiles, thermophiles and mesophiles. It displays the ratio between glutamic acid plus lysine (E+K) and glutamine plus histidine (Q+H), which is higher in thermophilic proteomes and thermostable proteins than in mesophilic proteomes and thermo labile proteins. Thermo-search allows a rapid screen of the CRM database for thermostable proteins in their functional categories and a visualization of the (E+K)/(Q+H) average ratio between organisms, allowing a comparison of their lifestyles.     

Prediction of potential thermostable proteins in Xylella fastidiosa

The average protein (E+K)/(Q+H) ratio in organisms has already been demonstrated to have a strong correlation with their optimal growth temperature. Employing the Thermo-Search web tool, we used this ratio as a basis to look for thermostable proteins in a mesophile, Xylella fastidiosa. Nine proteins were chosen to have their three-dimensional structures modeled by homology, using mainly proteins from mesophiles as templates. Resulting models featured a high number of hydrophobic interactions, a property that has been previously associated with thermostability. These results demonstrate the interesting possibility of using the (E+K)/(Q+H) ratio to find individual thermostable proteins in mesophilic organisms.     

Global Patterns of Protein Domain Gain and Loss in Superkingdoms

Domains are modules within proteins that can fold and function independently and are evolutionarily conserved. Here we compared the usage and distribution of protein domain families in the free-living proteomes of Archaea, Bacteria and Eukarya and reconstructed species phylogenies while tracing the history of domain emergence and loss in proteomes. We show that both gains and losses of domains occurred frequently during proteome evolution. The rate of domain discovery increased approximately linearly in evolutionary time. Remarkably, gains generally outnumbered losses and the gain-to-loss ratios were much higher in akaryotes compared to eukaryotes. Functional annotations of domain families revealed that both Archaea and Bacteria gained and lost metabolic capabilities during the course of evolution while Eukarya acquired a number of diverse molecular functions including those involved in extracellular processes, immunological mechanisms, and cell regulation. Results also highlighted significant contemporary sharing of informational enzymes between Archaea and Eukarya and metabolic enzymes between Bacteria and Eukarya. Finally, the analysis provided useful insights into the evolution of species. The archaeal superkingdom appeared first in evolution by gradual loss of ancestral domains, bacterial lineages were the first to gain superkingdom-specific domains, and eukaryotes (likely) originated when an expanding proto-eukaryotic stem lineage gained organelles through endosymbiosis of already diversified bacterial lineages. The evolutionary dynamics of domain families in proteomes and the increasing number of domain gains is predicted to redefine the persistence strategies of organisms in superkingdoms, influence the make up of molecular functions, and enhance organismal complexity by the generation of new domain architectures. This dynamics highlights ongoing secondary evolutionary adaptations in akaryotic microbes, especially Archaea.     

Correlation between codon usage and thermostability

Variations of arginine codon usage between organisms may have important implications to thermostability. The preferential usage of AGR codons for arginine in thermophiles and hyperthermophiles implies positive error minimization, contributing to avoid mutations that could harm protein thermostability. This bias is not a mere consequence of increased G + C content, as it has been previously suggested, and may represent a new mechanism of adaptation to protein thermostability.     

A multi-factors rational design strategy for enhancing the thermostability of Escherichia coli AppA phytase

Despite recent advances in our understanding of the importance of protein surface properties for protein thermostability,there are seldom studies on multi-factors rational design strategy, so a more scientific, simple and effective rational strategy is urgent for protein engineering. Here, we first attempted to use a three-factors rational design strategy combining three common structural features, protein flexibility, protein surface, and salt bridges. Escherichia coli AppA phytase was used as a model enzyme to improve its thermostability. Moreover, the structure and enzyme features of the thermostable mutants designed by our strategy were analyzed roundly. For the single mutants, two (Q206E and Y311K), in five exhibited thermostable property with a higher success rate of prediction (40 %). For the multiple mutants, the themostable sites were combined with another site, I427L, we obtained by directed evolution, Q206E/I427L, Y311K/I427L, and Q206E/Y311K/I427L, all exhibited thermostable property. The Y311K/I427L doubled thermostability (61.7 %, and was compared to 30.97 % after being heated at 80 °C for 10 min) and catalytic efficiency (4.46 was compared to 2.37) improved more than the wild-type AppA phytase almost without hampering catalytic activity. These multi-factors of rational design strategy can be applied practically as a thermostabilization strategy instead of the conventional single-factor approach.     

Paradoxical thermostable enzymes from psychrophile: molecular characterization and potentiality for biotechnological application

NAD(P)⁺-dependent aldehyde dehydrogenase (EC 1.2.1.5) and aspartase (EC 4.3.1.1) in the cells of an atypical psychrophile from Antarctic seawater, Cytophaga sp. KUC-1, were paradoxically thermostable, although they derived from a psychrophile. Both enzymes showed the highest activity at about 55 °C, and also active even under cold conditions. The enzymes contained more Ile residues than the enzymes from mesophiles. The Ile/Ile + Val + Leu ratio of the Cytophaga thermostable enzymes was much higher than that of the enzymes from mesophiles. As compared with the enzymes from other microorganisms, the Cytophaga thermostable enzymes have the structural differences in the C-terminal region of the enzymes. Therefore, the C-terminal region might be important for the paradoxical thermostability of the enzymes. The psychrophilic microorganism produces not only psychrophilic enzyme, but thermostable enzyme with psychrophilicity. Therefore, the psychrophilic microorganism is one of the candidates for isolation of novel biocatalysts, which have potential for various industrial applications.     

Impact of translational selection on codon usage bias in the archaeon Methanococcus maripaludis

Patterns of codon usage have been extensively studied among Bacteria and Eukaryotes, but there has been little investigation of species from the third domain of life, the Archaea. Here, we examine the nature of codon usage bias in a methanogenic archaeon, Methanococcus maripaludis. Genome-wide patterns of codon usage are dominated by a strong A + T bias, presumably largely reflecting mutation patterns. Nevertheless, there is variation among genes in the use of a subset of putatively translationally optimal codons, which is strongly correlated with gene expression level. In comparison with Bacteria such as Escherichia coli, the strength of selected codon usage bias in highly expressed genes in M. maripaludis seems surprisingly high given its moderate growth rate. However, the pattern of selected codon usage differs between M. maripaludis and E. coli: in the archaeon, strongly selected codon usage bias is largely restricted to twofold degenerate amino acids (AAs). Weaker bias among the codons for fourfold degenerate AAs is consistent with the small number of tRNA genes in the M. maripaludis genome.     

Presence of tRNA-dependent pathways correlates with high cysteine content in methanogenic Archaea

Archeal proteomes can be clustered into two groups based on their cysteine content. One group of proteomes displays a low cysteine content ( approximately 0.7% of the entire proteome), whereas the second group contains twice as many cysteines as the first ( approximately 1.3%). All cysteine-rich organisms belong to the methanogenic Archaea, which generates special cysteine clusters associated with primitive metabolic reactions. Our findings suggest that cysteine plays an important role in early forms of life.     

A key role in catalysis for His89 of adenylosuccinate lyase of Bacillus subtilis

Adenylosuccinate lyase of Bacillus subtilis is a tetrameric enzyme which catalyzes the cleavage of adenylosuccinate to AMP and fumarate. We have mutated His(89), one of three conserved histidines, to Gln, Ala, Glu, and Arg. The enzymes were expressed in Escherichia coli and purified to homogeneity. As compared to a specific activity of 1. 56 micromol of adenylosuccinate converted/min/mg protein for wild-type enzyme, the mutant enzymes exhibit specific activities of 0.0225, 0.0036, 0.0036, and 0.0009 for H89Q, H89A, H89E, and H89R, respectively. Circular dichroism and FPLC gel filtration reveal that mutant enzymes have a similar conformation and oligomeric state to that of wild-type enzyme. In H89Q, the K(M) for adenylosuccinate increases slightly to 2.5-fold that of wild-type, the K(M) for fumarate is elevated 3.3-fold, and the K(M) for AMP is 13 times higher than that observed in wild-type enzyme. The catalytic efficiency of the H89Q enzyme is compromised, with k(cat)/K(M) reduced 174-fold in the direction of AMP formation. These data suggest that His(89) plays a role in both the binding of the AMP portion of the substrate and in correctly orienting the substrate for catalysis. Incubation of H89Q with inactive H141Q enzyme [Lee, T. T., Worby, C., Bao, Z.-Q., Dixon, J. E., and Colman, R. F. (1999) Biochemistry 38, 22-32] leads to a 30-fold increase in activity. This intersubunit complementation indicates that His(89) and His(141) from different subunits participate in the active site and that both are required for catalysis.     

Bivalent cations and amino-acid composition contribute to the thermostability of Bacillus licheniformis xylose isomerase.

Comparative analysis of genome sequence data from mesophilic and hyperthermophilic micro-organisms has revealed a strong bias against specific thermolabile amino-acid residues (i.e. N and Q) in hyperthermophilic proteins. The N + Q content of class II xylose isomerases (XIs) from mesophiles, moderate thermophiles, and hyperthermophiles was examined. It was found to correlate inversely with the growth temperature of the source organism in all cases examined, except for the previously uncharacterized XI from Bacillus licheniformis DSM13 (BLXI), which had an N + Q content comparable to that of homologs from much more thermophilic sources. To determine whether BLXI behaves as a thermostable enzyme, it was expressed in Escherichia coli, and the thermostability and activity properties of the recombinant enzyme were studied. Indeed, it was optimally active at 70-72 degrees C, which is significantly higher than the optimal growth temperature (37 degrees C) of B. licheniformis. The kinetic properties of BLXI, determined at 60 degrees C with glucose and xylose as substrates, were comparable to those of other class II XIs. The stability of BLXI was dependent on the metallic cation present in its two metal-binding sites. The enzyme thermostability increased in the order apoenzyme < Mg2+-enzyme < Co2+-enzyme approximately Mn2+-enzyme, with melting temperatures of 50.3 degrees C, 53.3 degrees C, 73.4 degrees C, and 73.6 degrees C. BLXI inactivation was first-order in all conditions examined. The energy of activation for irreversible inactivation was also strongly influenced by the metal present, ranging from 342 kJ x mol(-1) (apoenzyme) to 604 kJ x mol(-1) (Mg2+-enzyme) to 1166 kJ x mol(-1) (Co2+-enzyme). These results suggest that the first irreversible event in BLXI unfolding is the release of one or both of its metals from the active site. Although N + Q content was an indicator of thermostability for class II XIs, this pattern may not hold for other sets of homologous enzymes. In fact, the extremely thermostable alpha-amylase from B. licheniformis was found to have an average N + Q content compared with homologous enzymes from a variety of mesophilic and thermophilic sources. Thus, it would appear that protein thermostability is a function of more complex molecular determinants than amino-acid content alone.     

Zinc through the three domains of life

Zinc is one of the metal ions essential for life, as it is required for the proper functioning of a large number of proteins. Despite its importance, the annotation of zinc-binding proteins in gene banks or protein domain databases still has significant room for improvement. In the present work, we compiled a list of known zinc-binding protein domains and of known zinc-binding sequence motifs (zinc-binding patterns), and then used them jointly to analyze the proteome of 57 different organisms to obtain an overview of zinc usage by archaeal, bacterial, and eukaryotic organisms. Zinc-binding proteins are an abundant fraction of these proteomes, ranging between 4% and 10%. The number of zinc-binding proteins correlates linearly with the total number of proteins encoded by the genome of an organism, but the proportionality constant of Eukaryota (8.8%) is significantly higher than that observed in Bacteria and Archaea (from 5% to 6%). Most of this enrichment is due to the larger portfolio of regulatory proteins in Eukaryota.     

Evolution of the mitochondrial genetic code I. Origin of AGR serine and stop codons in metazoan mitochondria

Summary AGA and AGG (AGR) are arginine codons in the universal genetic code. These codons are read as serine or are used as stop codons in metazoan mitochondria. The arginine residues coded by AGR in yeast orTrypanosoma are coded by arginine CGN throughout metazoan mitochondria. AGR serine sites in metazoan mitochondria are occupied mainly in corresponding sites in yeast orTrypanosoma mitochondria by UCN serine, AGY serine, or codons for amino acids other than serine or arginine. Based on these observations, we propose the following evolutionary events. AGR codons became unassigned because of deletion of tRNA Arg (UCU) and elimination of AGR codons by conversion to CGN arginine codons. Upon acquisition by serine tRNA of pairing ability with AGR codons, some codons for amino acids other than arginine mutated to AGR, and were caputed by anticodon GCU in serine tRNA. During vertebrate mitochondrial evolution, AGR stop codons presumably were created from UAG stop by deletion of the first nucleotide U and by use of R as the third nucleotide that had existed next to the ancestral UAG stop.     

Designing thermostable proteins: ancestral mutants of 3-isopropylmalate dehydrogenase designed by using a phylogenetic tree

We have recently developed a new method for designing thermostable proteins using phylogenetic trees of enzymes. In this study, we investigated a method for designing proteins with improved stability using 3-isopropylmalate dehydrogenase (IPMDH) from Thermus thermophilus as a model enzyme. We designed 12 mutant enzymes, each having an ancestral amino acid residue that was present in the common ancestor of Bacteria and Archaea. At least six of the 12 ancestral mutants tested showed thermal stability higher than that of the original enzyme. The results supported the hyperthermophilic universal ancestor hypothesis. The effect of ancestral residues on IPMDHs of several organisms and on the related enzyme isocitrate dehydrogenase was summarised and analysed. The effect of an ancestral residue on thermostability did not depend on the degree of conservation of the residue at the site, suggesting that the stabilisation of these mutant proteins is not related to sequence conservation but to the antiquity of the introduced residues. The results suggest also that this method could be an efficient way of designing mutant enzymes with higher thermostability based only on the primary structure and a phylogenetic tree.     

Occurrence of disordered patterns and homorepeats in eukaryotic and bacterial proteomes

Combining the motif discovery and disorder protein segment identification in PDB allows us to create the first and largest library of disordered patterns. At present the library includes 109 disordered patterns. Here we offer a comprehensive analysis of the occurrence of selected disordered patterns and 20 homorepeats of 6 residues long in 123 proteomes. 27 disordered patterns occur sparsely in all considered proteomes, but the patterns of low-complexity-homorepeats-appear more often in eukaryotic than in bacterial proteomes. A comparative analysis of the number of proteins containing homorepeats of 6 residues long and the disordered selected patterns in these proteomes has been performed. The matrices of correlation coefficients between numbers of proteins where at least once a homorepeat of six residues long for each of 20 types of amino acid residues and 109 disordered patterns from the library appears in 9 kingdoms of eukaryota and 5 phyla of bacteria have been calculated. As a rule, the correlation coefficients are higher inside the considered kingdom than between them. The largest fraction of homorepeats of 6 residues belongs to Amoebozoa proteomes (D. discoideum), 46%. Moreover, the longest uninterrupted repeats belong to S306 from D. discoideum (Amoebozoa). Homorepeats of some amino acids occur more frequently than others and the type of homorepeats varies across different proteomes, . For example, E6 appears most frequent for all considered proteomes for Chordata, Q6 for Arthropoda, S6 for Nematoda. The averaged occurrence of multiple long runs of 6 amino acids in a decreasing order for 97 eukaryotic proteomes is as follows: Q6, S6, A6, G6, N6, E6, P6, T6, D6, K6, L6, H6, R6, F6, V6, I6, Y6, C6, M6, W6, and for 26 bacterial proteomes it is A6, G6, P6, and the others occur seldom. This suggests that such short similar motifs are responsible for common functions for nonhomologous, unrelated proteins from different organisms.     

Consensus temporal order of amino acids and evolution of the triplet code

Forty different single-factor criteria and multi-factor hypotheses about chronological order of appearance of amino acids in the early evolution are summarized in consensus ranking. All available knowledge and thoughts about origin and evolution of the genetic code are thus combined in a single list where the amino acids are ranked chronologically. Due to consensus nature of the chronology it has several important properties not visible in individual rankings by any of the initial criteria. Nine amino acids of the Miller's imitation of primordial environment are all ranked as topmost (G, A, V, D, E, P, S, L, T). This result does not change even after several criteria related to Miller's data are excluded from calculations. The consensus order of appearance of the 20 amino acids on the evolutionary scene also reveals a unique and strikingly simple chronological organization of 64 codons, that could not be figured out from individual criteria: New codons appear in descending order of their thermostability, as complementary pairs, with the complements recruited sequentially from the codon repertoires of the earlier or simultaneously appearing amino acids. These three rules (Thermostability, Complementarity and Processivity) hold strictly as well as leading position of the earliest amino acids according to Miller. The consensus chronology of amino acids, G/A, V/D, P, S, E/L, T, R, N, K, Q, I, C, H, F, M, Y, W, and the derived temporal order for codons may serve, thus, as a justified working model of choice for further studies on the origin and evolution of the genetic code.     

A conformation-specific interhelical salt bridge in the K+ binding site of gastric H,K-ATPase

Homology modeling of gastric H,K-ATPase based on the E2 model of sarcoplasmic reticulum Ca2+-ATPase (Toyoshima, C., and Nomura, H. (2002) Nature 392, 835-839) revealed the presence of a single high-affinity binding site for K+ and an E2 form-specific salt bridge between Glu820 (M6) and Lys791 (M5). In the E820Q mutant this salt bridge is no longer possible, and the head group of Lys791, together with a water molecule, fills the position of the K+ ion and apparently mimics the K+-filled cation binding pocket. This gives an explanation for the K+-independent ATPase activity and dephosphorylation step of the E820Q mutant (Swarts, H. G. P., Hermsen, H. P. H., Koenderink, J. B., Schuurmans Stekhoven, F. M. A. H., and De Pont, J. J. H. H. M. (1998) EMBO J. 17, 3029-3035) and, indirectly, for its E1 preference. The model is strongly supported by a series of reported mutagenesis studies on charged and polar amino acid residues in the membrane domain. To further test this model, Lys791 was mutated alone and in combination with other crucial residues. In the K791A mutant, the K+ affinity was markedly reduced without altering the E2 preference of the enzyme. The K791A mutation prevented, in contrast to the K791R mutation, the spontaneous dephosphorylation of the E820Q mutant as well as its conformational equilibrium change toward E1. This indicates that the salt bridge is essential for high-affinity K+ binding and the E2 preference of H,K-ATPase. Moreover, its breakage (E820Q) can generate a K+-insensitive activity and an E1 preference. In addition, the study gives a molecular explanation for the electroneutrality of H,K-ATPases.     

Genes for tRNA(Arg) located in the upstream region of the Shiga toxin II operon in enterohemorrhagic Escherichia coli O157:H7.

We found two genes for tRNA(Arg) in the region upstream of genes for Shiga-like toxin type II (SLT-II) in Escherichia coli O157:H7. The two encoded forms of tRNA(Arg) recognize rare codons in E. coli K12 but these rare codons occur in the toxin genes at high frequency.     

Cooperative effect of two surface amino acid mutations (Q252L and E170K) in glucose dehydrogenase from Bacillus megaterium IWG3 on stabilization of its oligomeric state

A thermostable glucose dehydrogenase (GlcDH) mutant of Bacillus megaterium IWG3 harboring the Q252L substitution (Y. Makino, S. Negoro, I. Urabe, and H. Okada, J. Biol. Chem. 264:6381-6385, 1989) is stable at pH values above 9, but only in the presence of 2 M NaCl. Another GlcDH mutant exhibiting increased stability at an alkaline pH in the absence of NaCl has been isolated previously (S.-H. Baik, T. Ide, H. Yoshida, O. Kagami, and S. Harayama, Appl. Microbiol. Biotechnol. 61:329-335, 2003). This mutant had two amino acid substitutions, Q252L and E170K. In the present study, we characterized three GlcDH mutants harboring the substitutions Q252L, E170K, and Q252L/E170K under low-salt conditions. The GlcDH mutant harboring two substitutions, Q252L/E170K, was stable, but mutants harboring a single substitution, either Q252L or E170K, were unstable at an alkaline pH. Gel filtration chromatography analyses demonstrated that the oligomeric state of the Q252/E170K enzyme was stable, while the tetramers of the enzymes harboring a single substitution (Q252L or E170K) dissociated into dimers at an alkaline pH. These results indicated that the Q252L and E170K substitutions synergistically strengthened the interaction at the dimer-dimer interface. The crystal structure of the E170K/Q252L mutant, determined at 2.0-angstroms resolution, showed that residues 170 and 252 are located in a hydrophobic cavity at the subunit-subunit interface. We concluded that these residues in the wild-type enzyme have thermodynamically unfavorable effects, while the Q252L and E170K substitutions increase the subunit-subunit interactions by stabilizing the hydrophobic cavity.     

Mutational effects on thermostable superoxide dismutase from Aquifex pyrophilus: understanding the molecular basis of protein thermostability

We designed two mutants of superoxide dismutase (SOD), one is thermostable and the other is thermolabile, which provide valuable insight to identify amino acid residues essential for the thermostability of the SOD from Aquifex pyrophilus (ApSOD). The mutant K12A, in which Lys12 was replaced by Ala, had increased thermostability compared to that of the wild type. The T(1/2) value of K12A was 210 min and that of the wild type was 175 min at 95 degrees C. However, the thermostability of the mutant E41A, which has a T(1/2) value of 25 min at 95 degrees C, was significantly decreased compared to the wild type of ApSOD. To explain the enhanced thermostability of K12A and thermolabile E41A on the structural basis, the crystal structures of the two SOD mutants have been determined. The results have clearly shown the general significance of hydrogen bonds and ion-pair network in the thermostability of proteins.