Evolution of the aminoacyl-tRNA synthetases and the origin of the genetic code

The polyphyletic origins of glycyl-tRNA synthetase and lysyl-tRNA synthetase and their implications

I analyzed the polyphyletic origin of glycyl-tRNA synthetase (GlyRS) and lysyl-tRNA synthetase (LysRS), making plausible the following implications. The fact that the genetic code needed to evolve aminoacyl-tRNA synthetases (ARSs) only very late would be in perfect agreement with a late origin, in the main phyletic lineages, of both GlyRS and LysRS. Indeed, as suggested by the coevolution theory, since the genetic code was structured by biosynthetic relationships between amino acids and as these occurred on tRNA-like molecules which were evidently already loaded with amino acids during its structuring, this made possible a late origin of ARSs. All this corroborates the coevolution theory of the origin of the genetic code to the detriment of theories which would instead predict an early intervention of the action of ARSs in organizing the genetic code. Furthermore, the assembly of the GlyRS and LysRS protein domains in main phyletic lineages is itself at least evidence of the possibility that ancestral genes were assembled using pieces of genetic material that coded these protein domains. This is in accordance with the exon theory of genes which postulates that ancestral exons coded for protein domains or modules that were assembled to form the first genes. This theory is exemplified precisely in the evolution of both GlyRS and LysRS which occurred through the assembly of protein domains in the main phyletic lineages, as analyzed here. Furthermore, this late assembly of protein domains of these proteins into the two main phyletic lineages, i.e. a polyphyletic origin of both GlyRS and LysRS, appears to corroborate the progenote evolutionary stage for both LUCA and at least the first part of the evolutionary stages of the ancestor of bacteria and that of archaea. Indeed, this polyphyletic origin would imply that the genetic code was still evolving because at least two ARSs, i.e. proteins that make the genetic code possible today, were still evolving. This would imply that the evolutionary stages involved were characterized not by cells but by protocells, that is, by progenotes because this is precisely the definition of a progenote. This conclusion would be strengthened by the observation that both GlyRS and LysRS originating in the phyletic lineages leading to bacteria and archaea, would demonstrate that, more generally, proteins were most likely still in rapid and progressive evolution. Namely, a polyphyletic origin of proteins which would qualify at least the initial phase of the evolutionary stage of the ancestor of bacteria and that of archaea as stages belonging to the progenote.

Theories of the origin of the genetic code: Strong corroboration for the coevolution theory

I analyzed all the theories and models of the origin of the genetic code, and over the years, I have considered the main suggestions that could explain this origin. The conclusion of this analysis is that the coevolution theory of the origin of the genetic code is the theory that best captures the majority of observations concerning the organization of the genetic code. In other words, the biosynthetic relationships between amino acids would have heavily influenced the origin of the organization of the genetic code, as supported by the coevolution theory. Instead, the presence in the genetic code of physicochemical properties of amino acids, which have also been linked to the physicochemical properties of anticodons or codons or bases by stereochemical and physicochemical theories, would simply be the result of natural selection. More explicitly, I maintain that these correlations between codons, anticodons or bases and amino acids are in fact the result not of a real correlation between amino acids and codons, for example, but are only the effect of the intervention of natural selection. Specifically, in the genetic code table we expect, for example, that the most similar codons - that is, those that differ by only one base - will have more similar physicochemical properties. Therefore, the 64 codons of the genetic code table ordered in a certain way would also represent an ordering of some of their physicochemical properties. Now, a study aimed at clarifying which physicochemical property of amino acids has influenced the allocation of amino acids in the genetic code has established that the partition energy of amino acids has played a role decisive in this. Indeed, under some conditions, the genetic code was found to be approximately 98% optimized on its columns. In this same work, it was shown that this was most likely the result of the action of natural selection. If natural selection had truly allocated the amino acids in the genetic code in such a way that similar amino acids also have similar codons - this, not through a mechanism of physicochemical interaction between, for example, codons and amino acids - then it might turn out that even different physicochemical properties of codons (or anticodons or bases) show some correlation with the physicochemical properties of amino acids, simply because the partition energy of amino acids is correlated with other physicochemical properties of amino acids. It is very likely that this would inevitably lead to a correlation between codons (or anticodons or bases) and amino acids. In other words, since the codons (anticodons or bases) are ordered in the genetic code, that is to say, some of their physicochemical properties should also be ordered by a similar order, and given that the amino acids would also appear to have been ordered in the genetic code by selection natural, then it should inevitably turn out that there is a correlation between, for example, the hydrophobicity of anticodons and that of amino acids. Instead, the intervention of natural selection in organizing the genetic code would appear to be highly compatible with the main mechanism of structuring the genetic code as supported by the coevolution theory. This would make the coevolution theory the only plausible explanation for the origin of the genetic code.

Crick Wobble and Superwobble in Standard Genetic Code Evolution

Wobble coding is inevitable during evolution of the Standard Genetic Code (SGC). It ultimately splits half of NN U/C/A/G coding boxes with different assignments. Further, it contributes to pervasive SGC order by reinforcing close spacing for identical SGC assignments. But wobble cannot appear too soon, or it will inhibit encoding and more decisively, obstruct evolution of full coding tables. However, these prior results assumed Crick wobble, NN U/C and NN A/G, read by a single adaptor RNA. Superwobble translates NN U/C/A/G codons, using one adaptor RNA with an unmodified 5' anticodon U (appropriate to earliest coding) in modern mitochondria, plastids, and mycoplasma. Assuming the SGC was selected when evolving codes most resembled it, characteristics of the critical selection events can be calculated. For example, continuous superwobble infrequently evolves SGC-like coding tables. So, continuous superwobble is a very improbable origin hypothesis. In contrast, late-arising superwobble shares late Crick wobble's frequent resemblance to SGC order. Thus late superwobble is possible, but yields SGC-like assignments less frequently than late Crick wobble. Ancient coding ambiguity, most simply, arose from Crick wobble alone. This is consistent with SGC assignments to NAN codons.

Evolution of the genetic code

Diverse models have been advanced for the evolution of the genetic code. Here, models for tRNA, aminoacyl-tRNA synthetase (aaRS) and genetic code evolution were combined with an understanding of EF-Tu suppression of tRNA 3rd anticodon position wobbling. The result is a highly detailed scheme that describes the placements of all amino acids in the standard genetic code. The model describes evolution of 6-, 4-, 3-, 2- and 1-codon sectors. Innovation in column 3 of the code is explained. Wobbling and code degeneracy are explained. Separate distribution of serine sectors between columns 2 and 4 of the code is described. We conclude that very little chaos contributed to evolution of the genetic code and that the pattern of evolution of aaRS enzymes describes a history of the evolution of the code. A model is proposed to describe the biological selection for the earliest evolution of the code and for protocell evolution.

The phylogenetic distribution of the glutaminyl-tRNA synthetase and Glu-tRNAGln amidotransferase in the fundamental lineages would imply that the ancestor of archaea, that of eukaryotes and LUCA were progenotes

The function of the glutaminyl-tRNA synthetase and Glu-tRNA^Gln amidotransferase might be related to the origin of the genetic code because, for example, glutaminyl-tRNA synthetase catalyses the fundamental reaction that makes the genetic code. If the evolutionary stage of the origin of these two enzymes could be unambiguously identified, then the genetic code should still have been originating at that particular evolutionary stage because the fundamental reaction that makes the code itself was still evidently evolving. This would result in that particular evolutionary moment being attributed to the evolutionary stage of the progenote because it would have a relationship between the genotype and the phenotype not yet fully realized because the genetic code was precisely still originating. I then analyzed the distribution of the glutaminyl-tRNA synthetase and Glu-tRNA^Gln aminodotrasferase in the main phyletic lineages. Since in some cases the origin of these two enzymes can be related to the evolutionary stages of ancestors of archaea and eukaryotes, this would indicate these ancestors as progenotes because at that evolutionary moment the genetic code was evidently still evolving, thus realizing the definition of progenote. The conclusion that the ancestor of archaea and that of eukaryotes were progenotes would imply that even the last universal common ancestor (LUCA) was a progenote because it appeared, on the tree of life, temporally before these ancestors.

The Cipher of the Genetic Code

A new approach to understanding of the genetic code is developed. In order to overcome the key paradox (and Darwinian selection problem) that the highly complex amino acid Phe is encoded by the simplest codons (UUY), and the simplest Gly encoded by the most complex codons (GGN); as well as the paradox of the duplication of some amino acids in the encoding process (Leu, Ser, Arg), we proposed an extension of the notion (and concept) of genetic code. For a better (and lighter) understanding of genetic coding, we proposed a hypothesis after that (under the conditions of allowed metaphoricity and modeling in biology) genetic code has to be understood, analogously in cryptology, as the unity of three entities: the code, the cipher of the code and the key of the cipher. In this hierarchy the term (and notion) "genetic code" remains what has been from the beginning: a connection between four-letter alphabet (four Py-Pu nucleotides, in form of codons) and a twenty-letter alphabet (twenty amino acids); the cipher is a specific chemical complementarity in chemical properties of molecules in the form: similarity in dissimilarity versus dissimilarity in similarity ("Sim in Diss vs Diss in Sim") and the key of cipher: the complementarity on the binary tree of the genetic code in the form: 0-15, 1-14, 2-13, …, 6-9, 7-8. These concepts improve understanding that within the two main Genetic Code Tables (of the nucleotide doublets and nucleotide Triplets) exists a sophisticated nuancing and balancing in the properties of the constituents of GC, including the balance of the number of molecules, atoms, and nucleons.

Horizontal transfer of code fragments between protocells can explain the origins of the genetic code without vertical descent

Theories of the origin of the genetic code typically appeal to natural selection and/or mutation of hereditable traits to explain its regularities and error robustness, yet the present translation system presupposes high-fidelity replication. Woese's solution to this bootstrapping problem was to assume that code optimization had played a key role in reducing the effect of errors caused by the early translation system. He further conjectured that initially evolution was dominated by horizontal exchange of cellular components among loosely organized protocells ("progenotes"), rather than by vertical transmission of genes. Here we simulated such communal evolution based on horizontal transfer of code fragments, possibly involving pairs of tRNAs and their cognate aminoacyl tRNA synthetases or a precursor tRNA ribozyme capable of catalysing its own aminoacylation, by using an iterated learning model. This is the first model to confirm Woese's conjecture that regularity, optimality, and (near) universality could have emerged via horizontal interactions alone.

Evolutionary Limitation and Opportunities for Developing tRNA Synthetase Inhibitors with 5-Binding-Mode Classification

Aminoacyl-tRNA synthetases (aaRSs) are enzymes that catalyze the transfer of amino acids to their cognate tRNAs as building blocks for translation. Each of the aaRS families plays a pivotal role in protein biosynthesis and is indispensable for cell growth and survival. In addition, aaRSs in higher species have evolved important non-translational functions. These translational and non-translational functions of aaRS are attractive for developing antibacterial, antifungal, and antiparasitic agents and for treating other human diseases. The interplay between amino acids, tRNA, ATP, EF-Tu and non-canonical binding partners, had shaped each family with distinct pattern of key sites for regulation, with characters varying among species across the path of evolution. These sporadic variations in the aaRSs offer great opportunity to target these essential enzymes for therapy. Up to this day, growing numbers of aaRS inhibitors have been discovered and developed. Here, we summarize the latest developments and structural studies of aaRS inhibitors, and classify them with distinct binding modes into five categories.

Toward major evolutionary transitions theory 2.0

The impressive body of work on the major evolutionary transitions in the last 20 y calls for a reconstruction of the theory although a 2D account (evolution of informational systems and transitions in individuality) remains. Significant advances include the concept of fraternal and egalitarian transitions (lower-level units like and unlike, respectively). Multilevel selection, first without, then with, the collectives in focus is an important explanatory mechanism. Transitions are decomposed into phases of origin, maintenance, and transformation (i.e., further evolution) of the higher level units, which helps reduce the number of transitions in the revised list by two so that it is less top-heavy. After the transition, units show strong cooperation and very limited realized conflict. The origins of cells, the emergence of the genetic code and translation, the evolution of the eukaryotic cell, multicellularity, and the origin of human groups with language are reconsidered in some detail in the light of new data and considerations. Arguments are given why sex is not in the revised list as a separate transition. Some of the transitions can be recursive (e.g., plastids, multicellularity) or limited (transitions that share the usual features of major transitions without a massive phylogenetic impact, such as the micro- and macronuclei in ciliates). During transitions, new units of reproduction emerge, and establishment of such units requires high fidelity of reproduction (as opposed to mere replication).

Ancestral Reconstruction of a Pre-LUCA Aminoacyl-tRNA Synthetase Ancestor Supports the Late Addition of Trp to the Genetic Code

The genetic code was likely complete in its current form by the time of the last universal common ancestor (LUCA). Several scenarios have been proposed for explaining the code's pre-LUCA emergence and expansion, and the relative order of the appearance of amino acids used in translation. One co-evolutionary model of genetic code expansion proposes that at least some amino acids were added to the code by the ancient divergence of aminoacyl-tRNA synthetase (aaRS) families. Of all the amino acids used within the genetic code, Trp is most frequently claimed as a relatively recent addition. We observe that, since TrpRS and TyrRS are paralogous protein families retaining significant sequence similarity, the inferred sequence composition of their ancestor can be used to evaluate this co-evolutionary model of genetic code expansion. We show that ancestral sequence reconstructions of the pre-LUCA paralog ancestor of TyrRS and TrpRS have several sites containing Tyr, yet a complete absence of sites containing Trp. This is consistent with the paralog ancestor being specific for the utilization of Tyr, with Trp being a subsequent addition to the genetic code facilitated by a process of aaRS divergence and neofunctionalization. Only after this divergence could Trp be specifically encoded and incorporated into proteins, including the TyrRS and TrpRS descendant lineages themselves. This early absence of Trp is observed under both homogeneous and non-homogeneous models of ancestral sequence reconstruction. Simulations support that this observed absence of Trp is unlikely to be due to chance or model bias. These results support that the final stages of genetic code evolution occurred well within the "protein world," and that the presence-absence of Trp within conserved sites of ancient protein domains is a likely measure of their relative antiquity, permitting the relative timing of extremely early events within protein evolution before LUCA.

Probable relationship between partitions of the set of codons and the origin of the genetic code

Here we study the distribution of randomly generated partitions of the set of amino acid-coding codons. Some results are an application from a previous work, about the Stirling numbers of the second kind and triplet codes, both to the cases of triplet codes having four stop codons, as in mammalian mitochondrial genetic code, and hypothetical doublet codes. Extending previous results, in this work it is found that the most probable number of blocks of synonymous codons, in a genetic code, is similar to the number of amino acids when there are four stop codons, as well as it could be for a primigenious doublet code. Also it is studied the integer partitions associated to patterns of synonymous codons and it is shown, for the canonical code, that the standard deviation inside an integer partition is one of the most probable. We think that, in some early epoch, the genetic code might have had a maximum of the disorder or entropy, independent of the assignment between codons and amino acids, reaching a state similar to "code freeze" proposed by Francis Crick. In later stages, maybe deterministic rules have reassigned codons to amino acids, forming the natural codes, such as the canonical code, but keeping the numerical features describing the set partitions and the integer partitions, like a "fossil numbers"; both kinds of partitions about the set of amino acid-coding codons.

The protein invasion: a broad review on the origin of the translational system

Translation, coded peptide synthesis, arguably exists at the heart of modern cellular life. By orchestrating an incredibly complex interaction between tRNAs, mRNAs, aaRSs, the ribosome, and numerous other small molecules, the translational system allows the interpretation of data in the form of DNA to create massively complex proteins which control and enact almost every cellular function. A natural question then, is how did this system evolve? Here we present a broad review of the existing theories of the last two decades on the origin of the translational system. We attempt to synthesize the wide variety of ideas as well as organize them into modular components, addressing the evolution of the peptide-RNA interaction, tRNA, mRNA, the ribosome, and the first proteins separately. We hope to provide both a comprehensive overview of the literature as well as a framework for future discussions and novel theories.

On origin of genetic code and tRNA before translation

BACKGROUND

Synthesis of proteins is based on the genetic code - a nearly universal assignment of codons to amino acids (aas). A major challenge to the understanding of the origins of this assignment is the archetypal "key-lock vs. frozen accident" dilemma. Here we re-examine this dilemma in light of 1) the fundamental veto on "foresight evolution", 2) modular structures of tRNAs and aminoacyl-tRNA synthetases, and 3) the updated library of aa-binding sites in RNA aptamers successfully selected in vitro for eight amino acids.

RESULTS

The aa-binding sites of arginine, isoleucine and tyrosine contain both their cognate triplets, anticodons and codons. We have noticed that these cases might be associated with palindrome-dinucleotides. For example, one-base shift to the left brings arginine codons CGN, with CG at 1-2 positions, to the respective anticodons NCG, with CG at 2-3 positions. Formally, the concomitant presence of codons and anticodons is also expected in the reverse situation, with codons containing palindrome-dinucleotides at their 2-3 positions, and anticodons exhibiting them at 1-2 positions. A closer analysis reveals that, surprisingly, RNA binding sites for Arg, Ile and Tyr "prefer" (exactly as in the actual genetic code) the anticodon(2-3)/codon(1-2) tetramers to their anticodon(1-2)/codon(2-3) counterparts, despite the seemingly perfect symmetry of the latter. However, since in vitro selection of aa-specific RNA aptamers apparently had nothing to do with translation, this striking preference provides a new strong support to the notion of the genetic code emerging before translation, in response to catalytic (and possibly other) needs of ancient RNA life. Consistently with the pre-translation origin of the code, we propose here a new model of tRNA origin by the gradual, Fibonacci process-like, elongation of a tRNA molecule from a primordial coding triplet and 5'DCCA3' quadruplet (D is a base-determinator) to the eventual 76 base-long cloverleaf-shaped molecule.

CONCLUSION

Taken together, our findings necessarily imply that primordial tRNAs, tRNA aminoacylating ribozymes, and (later) the translation machinery in general have been co-evolving to ''fit'' the (likely already defined) genetic code, rather than the opposite way around. Coding triplets in this primal pre-translational code were likely similar to the anticodons, with second and third nucleotides being more important than the less specific first one. Later, when the code was expanding in co-evolution with the translation apparatus, the importance of 2-3 nucleotides of coding triplets "transferred" to the 1-2 nucleotides of their complements, thus distinguishing anticodons from codons. This evolutionary primacy of anticodons in genetic coding makes the hypothesis of primal stereo-chemical affinity between amino acids and cognate triplets, the hypothesis of coding coenzyme handles for amino acids, the hypothesis of tRNA-like genomic 3' tags suggesting that tRNAs originated in replication, and the hypothesis of ancient ribozymes-mediated operational code of tRNA aminoacylation not mutually contradicting but rather co-existing in harmony.

Origin and evolution of the genetic code: the universal enigma

The genetic code is nearly universal, and the arrangement of the codons in the standard codon table is highly nonrandom. The three main concepts on the origin and evolution of the code are the stereochemical theory, according to which codon assignments are dictated by physicochemical affinity between amino acids and the cognate codons (anticodons); the coevolution theory, which posits that the code structure coevolved with amino acid biosynthesis pathways; and the error minimization theory under which selection to minimize the adverse effect of point mutations and translation errors was the principal factor of the code's evolution. These theories are not mutually exclusive and are also compatible with the frozen accident hypothesis, that is, the notion that the standard code might have no special properties but was fixed simply because all extant life forms share a common ancestor, with subsequent changes to the code, mostly, precluded by the deleterious effect of codon reassignment. Mathematical analysis of the structure and possible evolutionary trajectories of the code shows that it is highly robust to translational misreading but there are numerous more robust codes, so the standard code potentially could evolve from a random code via a short sequence of codon series reassignments. Thus, much of the evolution that led to the standard code could be a combination of frozen accident with selection for error minimization although contributions from coevolution of the code with metabolic pathways and weak affinities between amino acids and nucleotide triplets cannot be ruled out. However, such scenarios for the code evolution are based on formal schemes whose relevance to the actual primordial evolution is uncertain. A real understanding of the code origin and evolution is likely to be attainable only in conjunction with a credible scenario for the evolution of the coding principle itself and the translation system.

In silico detection of tRNA sequence features characteristic to aminoacyl-tRNA synthetase class membership

Aminoacyl tRNA synthetases (aaRS) are grouped into Class I and II based on primary and tertiary structure and enzyme properties suggesting two independent phylogenetic lineages. Analogously, tRNA molecules can also form two respective classes, based on the class membership of their corresponding aaRS. Although some aaRS–tRNA interactions are not extremely specific and require editing mechanisms to avoid misaminoacylation, most aaRS–tRNA interactions are rather stereospecific. Thus, class-specific aaRS features could be mirrored by class-specific tRNA features. However, previous investigations failed to detect conserved class-specific nucleotides. Here we introduce a discrete mathematical approach that evaluates not only class-specific ‘strictly present’, but also ‘strictly absent’ nucleotides. The disjoint subsets of these elements compose a unique partition, named extended consensus partition (ECP). By analyzing the ECP for both Class I and II tDNA sets from 50 (13 archaeal, 30 bacterial and 7 eukaryotic) species, we could demonstrate that class-specific tRNA sequence features do exist, although not in terms of strictly conserved nucleotides as it had previously been anticipated. This finding demonstrates that important information was hidden in tRNA sequences inaccessible for traditional statistical methods. The ECP analysis might contribute to the understanding of tRNA evolution and could enrich the sequence analysis tool repertoire.

Extremely thermophilic translation system in the common ancestor commonote: ancestral mutants of Glycyl-tRNA synthetase from the extreme thermophile Thermus thermophilus

Based on phylogenetic analysis of 16 S and 18 S rRNAs, the common ancestor of all organisms (Commonote) was proposed to be hyperthermophilic. We have previously tested this hypothesis using enzymes with ancestral residues that are inferred by molecular phylogenetic analysis. The ancestral mutant enzymes involved in metabolic systems show higher thermal stability than wild-type enzymes, consistent with the hyperthermophile common ancestor hypothesis. Here, we have extended the experiments to include an enzyme of the translation system, glycyl-tRNA synthetase (GlyRS). The translation system often shows a phylogenetic tree that is similar to the rRNA tree. Thus, it is likely that the tree represents the evolutionary route of the organisms. The maximum-likelihood tree of alpha(2) type GlyRS was constructed. From this analysis the ancestral sequence of GlyRS was deduced and individual or pairs of ancestral residues were introduced into Thermus thermophilus GlyRS. The ancestral mutants were expressed in Escherichia coli, purified and activity measured. The thermostability of eight mutated proteins was evaluated by CD (circular dichroism) measurements. Six mutants showed higher thermostability than wild-type enzyme and seven mutants showed higher activity than wild-type enzyme at 70 degrees C, suggesting an extremely thermophilic translation system in the common ancestor Commonote.

An examination of the carbon isotope effects associated with amino acid biosynthesis

Stable carbon isotope ratios (delta(13)C) were determined for alanine, proline, phenylalanine, valine, leucine, isoleucine, aspartate (aspartic acid and asparagine), glutamate (glutamic acid and glutamine), lysine, serine, glycine, and threonine from metabolically diverse microorganisms. The microorganisms examined included fermenting bacteria, organotrophic, chemolithotrophic, phototrophic, methylotrophic, methanogenic, acetogenic, acetotrophic, and naturally occurring cryptoendolithic communities from the Dry Valleys of Antarctica. Here we demonstrated that reactions involved in amino acid biosynthesis can be used to distinguish amino acids formed by life from those formed by nonbiological processes. The unique patterns of delta(13)C imprinted by life on amino acids produced a biological bias. We also showed that, by applying discriminant function analysis to the delta(13)C value of a pool of amino acids formed by biological activity, it was possible to identify key aspects of intermediary carbon metabolism in the microbial world. In fact, microorganisms examined in this study could be placed within one of three metabolic groups: (1) heterotrophs that grow by oxidizing compounds containing three or more carbon-to-carbon bonds (fermenters and organotrophs), (2) autotrophs that grow by taking up carbon dioxide (chemolitotrophs and phototrophs), and (3) acetoclastic microbes that grow by assimilation of formaldehyde or acetate (methylotrophs, methanogens, acetogens, and acetotrophs). Furthermore, we demonstrated that cryptoendolithic communities from Antarctica grouped most closely with the autotrophs, which indicates that the dominant metabolic pathways in these communities are likely those utilized for CO(2 )fixation. We propose that this technique can be used to determine the dominant metabolic types in a community and reveal the overall flow of carbon in a complex ecosystem.

The triplet genetic code had a doublet predecessor

Information theoretic analysis of genetic languages indicates that the naturally occurring 20 amino acids and the triplet genetic code arose by duplication of 10 amino acids of class-II and a doublet genetic code having codons NNY and anticodons GNN. Evidence for this scenario is presented based on the properties of aminoacyl-tRNA synthetases, amino acids and nucleotide bases.

A harmonic structure of the genetic code

In this paper is presented a new, very harmonic structure of the genetic code (GC) within a system of "4 x 5" (and/or of "5 x 4") of amino acids (AAs) in two variants. In first variant, the five rows within the system start with one polar charged amino acid (AA) each, making first column, consisting from five polar charged AAs (D, R, K, H, E). Five polar non-charged AAs (N, P, Y, W, Q) follow, then five non-polar AAs as last column (A, L, F, V, I) and, finally, five polar or non-polar AAs, in a combination, as first to last column (A as non-polar; S, T as polar, and G, P as ambivalent AAs). A second variant is subsequent to this one-"4 x 5" system with five nitrogen AAs (K, R, P, H, W), five oxygen (D, E, Y, S, T), five solely carbon (A, L, F, V, I) and five "combined" AAs (G with hydrogen as side chain; C and M with carbon and sulfur; N and Q with carbon, oxygen and nitrogen). A strict balance of atom and nucleon number as well as molecule mass follows the classification in both system variants.

Amino acid biogenesis, evolution of the genetic code and aminoacyl-tRNA synthetases

The aminoacyl-tRNA synthetases (aaRSs) ensure the fidelity of the translation of the genetic code, covalently attaching appropriate amino acids to the corresponding nucleic acid adaptor molecules-tRNA. The fundamental role of aminoacylation reaction catalysed by aaRSs implies that representatives of the family are thought to be among the earliest proteins to appear. Based on sequence analysis and catalytic domain structure, aaRSs have been partitioned into two classes of 10 enzymes each. However, based on the structural and sequence data only, it will not be easily understood that the present partitioning is not governed by chance. Our findings suggest that organization of amino acid biosynthetic pathways and clustering of aaRSs into different classes are intimately related to one another. A plausible explanation for such a relationship is dictated by early link between aaRSs and amino acids biosynthetic proteins. The aaRSs catalytic cores are highly relevant to the ancient metabolic reactions, namely, amino acids and cofactors biosynthesis. In particular we show that class II aaRSs mostly associated with the primordial amino acids, while class I aaRSs are usually related to amino acids evolved lately. Reasoning from this we propose a possible chronology of genetic code evolution.

On the relative content of G,C bases in codons of amino acids corresponding to class I and II aminoacyl-tRNA synthetases

We have analyzed the relative G,C content from protein coding regions of 530 organisms and found that the ratio of the G,C content of the codons of the amino acids corresponding to Class II and Class I aminoacyl-tRNA synthetases decreases in a statistically significant way from prokaryotes to animals. This can be interpreted assuming that an initial asymmetry between the G,C content of codons of Class I and II amino acids existed and has decreased in the course of evolution.

On the classes of aminoacyl-tRNA synthetases and the error minimization in the genetic code

As a consequence of the existence of two classes of aminoacyl-tRNA synthetases (aaRSs), we defined two types of mutations: g (mutations that do not change the class of the involved amino acids) and u (those which change the class). We have found that the mean chemical distance resulting from g mutations is smaller than that corresponding to u mutations, indicating that g mutations are responsible for most of the known minimization of the genetic code. This supports models for the origin and evolution of the code, in which new amino acids were added after duplications or modification of existing aaRSs.

Early fixation of an optimal genetic code

The evolutionary forces that produced the canonical genetic code before the last universal ancestor remain obscure. One hypothesis is that the arrangement of amino acid/codon assignments results from selection to minimize the effects of errors (e.g., mistranslation and mutation) on resulting proteins. If amino acid similarity is measured as polarity, the canonical code does indeed outperform most theoretical alternatives. However, this finding does not hold for other amino acid properties, ignores plausible restrictions on possible code structure, and does not address the naturally occurring nonstandard genetic codes. Finally, other analyses have shown that significantly better code structures are possible. Here, we show that if theoretically possible code structures are limited to reflect plausible biological constraints, and amino acid similarity is quantified using empirical data of substitution frequencies, the canonical code is at or very close to a global optimum for error minimization across plausible parameter space. This result is robust to variation in the methods and assumptions of the analysis. Although significantly better codes do exist under some assumptions, they are extremely rare and thus consistent with reports of an adaptive code: previous analyses which suggest otherwise derive from a misleading metric. However, all extant, naturally occurring, secondarily derived, nonstandard genetic codes do appear less adaptive. The arrangement of amino acid assignments to the codons of the standard genetic code appears to be a direct product of natural selection for a system that minimizes the phenotypic impact of genetic error. Potential criticisms of previous analyses appear to be without substance. That known variants of the standard genetic code appear less adaptive suggests that different evolutionary factors predominated before and after fixation of the canonical code. While the evidence for an adaptive code is clear, the process by which the code achieved this optimization requires further attention.

A gene fusion event in the evolution of aminoacyl-tRNA synthetases

The genes of glutamyl- and prolyl-tRNA synthetases (GluRS and ProRS) are organized differently in the three kingdoms of the tree of life. In bacteria and archaea, distinct genes encode the two proteins. In several organisms from the eukaryotic phylum of coelomate metazoans, the two polypeptides are carried by a single polypeptide chain to form a bifunctional protein. The linker region is made of imperfectly repeated units also recovered as singular or plural elements connected as N-terminal or C-terminal polypeptide extensions in various eukaryotic aminoacyl-tRNA synthetases. Phylogenetic analysis points to the monophyletic origin of this polypeptide motif appended to six different members of the synthetase family, belonging to either of the two classes of aminoacyl-tRNA synthetases. In particular, the monospecific GluRS and ProRS from Caenorhabditis elegans, an acoelomate metazoan, exhibit this recurrent motif as a C-terminal or N-terminal appendage, respectively. Our analysis of the extant motifs suggests a possible series of events responsible for a gene fusion that gave rise to the bifunctional glutamyl-prolyl-tRNA synthetase through recombination between genomic sequences encoding the repeated units.

Physico-chemical constraints connected with the coding properties of the genetic system

New insights on the origin of the genetic code, based on the analysis of the physico-chemical properties of its molecular constituents (RNA and amino acids), are reported in this paper. We point out a symmetry in the genetic code table and show that it can be explained by the nature of the anticodon-codon interaction. The importance of the strength of this interaction is examined and a correlation is found between the free-energy change (DeltaG(0)) of anticodon-codon association and the volume of the corresponding amino acids. This correlation is investigated in conjunction with the well-known one linking the hydrophobicity of the anticodons with that of the amino acids. We show that they can be considerated separately and that the energy vs. volume correlation may be explained by the process implicating the peptide bond formation between two successive amino acids during translation. This interpretation is supported by a statistical pattern of bases (purines or pyrimidines), observed in present coding genes, and by considerations involving the availability of the different kinds of amino acids. Finally, we try to explain the hydrophobicity correlation when reconstructing the events at the time of the so-called "RNA World". The whole of our investigation shows that the genetic code might be sufficiently robust to exist without the participation of pre-existing proteins, and that this robustness is a consequence of the physico-chemical properties of the four bases of the genetic system.

A possible mechanism of peptide bond formation on ribosome without mediation of peptidyl transferase

Ribosome, the ubiquitous organelle, is the site for protein synthesis in all types of cells. The consecutive peptide bonds are formed by the transpeptidation reaction between carboxyl group of peptidyl moiety and the amino group of the aminoacyl moiety. Both the moieties are attached to the appropiate tRNAs positioned on the ribosome at P and A sites, respectively, through codon-anticodon recognition directed by messenger RNA. The reaction seems to proceed by the nucleophillic attack of the amino group of the aminoacyl tRNA at the A site and on the carboxyl of the ester group of the tRNA at P-site of ribosome. The configuration of the carbon atom of the tetrahedral intermediate may be R or S depending on the direction of the nucleophillic attack. After selecting the favorable conformation of this tetrahedral intermediate quantum mechanical calculations have been carried out to determine the energy needed for its formation. A cyclic intermediate where 2'-OH of the ribose sugar of the P-site tRNA is a member of the ring can be formed from the tetrahedral intermediate. This cyclic intermediate produces a free tRNA and a tRNA attached to a planar peptide unit. Analysis of the energetics using semiempirical method for the formation of a cyclic intermediate indicates that the peptide bond formation through the tetrahedral intermediate in S configuration may not need assistance from any outside agent like an enzyme

The genetic code as a Golden mean determined system

For the first time it is shown that the genetic code, as a binary code, is determined by Golden mean through the unity of the binary-code tree and the Farey tree, corresponding to the Watson-Crick table and Gray code model of the genetic code at the same time. On the other hand, the Golden mean relations also correspond with the essentiality-nonessentiality of protein amino acids, as well as with their splitting into two classes handled by two classes of aminoacyl-tRNA synthetases.

Gene families: the taxonomy of protein paralogs and chimeras

Ancient duplications and rearrangements of protein-coding segments have resulted in complex gene family relationships. Duplications can be tandem or dispersed and can involve entire coding regions or modules that correspond to folded protein domains. As a result, gene products may acquire new specificities, altered recognition properties, or modified functions. Extreme proliferation of some families within an organism, perhaps at the expense of other families, may correspond to functional innovations during evolution. The underlying processes are still at work, and the large fraction of human and other genomes consisting of transposable elements may be a manifestation of the evolutionary benefits of genomic flexibility.

Partition of aminoacyl-tRNA synthetases in two different structural classes dating back to early metabolism: implications for the origin of the genetic code and the nature of protein sequences

We describe, on the molecular level, a possible fuzzy and primordial translation apparatus capable of synthesizing polypeptides from nucleic acids in a world containing a mixture of coevolving molecules of RNA and proteins already arranged in metabolic cycles (including cofactors). Close attention is paid to template-free systems because they are believed to be the immediate ancestors of this primordial translation apparatus. The two classes of aminoacyl-tRNA synthetases (aaRSs), as seen today, are considered as the remnants of such a simple imprecise translation apparatus and are used as guidelines for the construction of the model. Earlier theoretical work by Bedian on a related system is invoked to show how specificity and stability could have been achieved automatically and rather quickly, starting from such an imprecise system, i.e., how the encoded synthesis of proteins could have appeared. Because of the binary nature of the underlying proto-code, the first genetically encoded proteins would then have been alternating copolymers with a high degree of degeneracy, but not random. Indeed, a clear signal for alternating hydrophobic and hydrophilic residues in present-day protein sequences can be detected. Later evolution of the genetic code would have proceeded along lines already discussed by Crick. However, in the initial stages, the translation apparatus proposed here is in fact very similar to the one postulated by Woese, only here it is given a molecular framework. This hypothesis departs from the paradigm of the RNA world in that it supposes that the origin of the genetic code occurred after the apparition of some functional (statistical) proteins first. Implications for protein design are also discussed.