Testing a biosynthetic theory of the genetic code: fact or artifact?

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.

What Would an Alien Amino Acid Alphabet Look Like and Why?

Life on Earth builds genetically encoded proteins by using a standard alphabet of just 20 L-α-amino acids, although many others were available to life's origins and early evolution. To better understand the causes of this foundational evolutionary outcome, we extend previous analyses which have identified a highly unusual distribution of biophysical properties within the set used by life. Specifically, we use a heuristic search algorithm to identify other sets of amino acids, from a library of plausible alternatives, that emulate life's signature. We find that a subset of amino acids seems predisposed to forming such sets. We present other examples of such alphabets under various assumptions, along with analysis and reasoning about why each might be simplistic. We do so to introduce the central, open question that remains: while fundamental biophysics related to protein folding can potentially reduce a library of 10⁵⁴ possible amino acid alphabets by 7 orders of magnitude, the framework of assumptions that does so leaves a further 10⁴⁵ possibilities. It is therefore tempting to ask what additional assumptions can further reduce these 45 orders of magnitude? We thus conclude with a focus on library and alphabet construction as a useful target for subsequent research that may help future science speak with more confidence about what an alien amino acid alphabet would look like and why.

The Origin of Genetic Code and Translation in the Framework of Current Concepts on the Origin of Life

The origin of genetic code and translation system is probably the central and most difficult problem in the investigations on the origin of life and one of the most complex problems in the evolutionary biology in general. There are multiple hypotheses on the emergence and development of existing genetic systems that propose the mechanisms for the origin and early evolution of genetic code, as well as for the emergence of replication and translation. Here, we discuss the most well-known of these hypotheses, although none of them provides a description of the early evolution of genetic systems without gaps and assumptions. The RNA world hypothesis is a currently prevailing scientific idea on the early evolution of biological and pre-biological structures, the main advantage of which is the assumption that RNAs as the first living systems were self-sufficient, i.e., capable of functioning as both catalysts and templates. However, this hypothesis has also significant limitations. In particular, no ribozymes with processive polymerase activity have been yet discovered or synthesized. Taking into account the mutual need of proteins and nucleic acids in each other in the current world, many authors propose the early evolution scenarios based on the co-evolution of these two classes of organic molecules. They postulate that the emergence of translation was necessary for the replication of nucleic acids, in contrast to the RNA world hypothesis, according to which the emergence of translation was preceded by the era of self-replicating RNAs. Although such scenarios are less parsimonious from the evolutionary point of view, since they require simultaneous emergence and evolution of two classes of organic molecules, as well as the emergence of synchronized replication and translation, their major advantage is that they explain the development of processive and much more accurate protein-dependent replication.

Fitting the standard genetic code into its triplet table

Minimally evolved codes are constructed here; these have randomly chosen standard genetic code (SGC) triplets, completed with completely random triplet assignments. Such "genetic codes" have not evolved, but retain SGC qualities. Retained qualities are basic, part of the underpinning of coding. For example, the sensitivity of coding to arbitrary assignments, which must be < ∼10%, is intrinsic. Such sensitivity comes from the elementary combinatorial properties of coding and constrains any SGC evolution hypothesis. Similarly, assignment of last-evolved functions is difficult because of late kinetic phenomena, likely common across codes. Census of minimally evolved code assignments shows that shape and size of wobble domains controls the code's fit into a coding table, strongly shifting accuracy of codon assignments. Access to the SGC therefore requires a plausible pathway to limited randomness, avoiding difficult completion while fitting a highly ordered, degenerate code into a preset three-dimensional space. Three-dimensional late Crick wobble in a genetic code assembled by lateral transfer between early partial codes satisfies these varied, simultaneous requirements. By allowing parallel evolution of SGC domains, this origin can yield shortened evolution to SGC-level order and allow the code to arise in smaller populations. It effectively yields full codes. Less obviously, it unifies previously studied chemical, biochemical, and wobble order in amino acid assignment, including a stereochemical minority of triplet-amino acid associations. Finally, fusion of intermediates into the final SGC is credible, mirroring broadly accepted later cellular evolution.

Optimal Evolution of the Standard Genetic Code

The Standard Genetic Code (SGC) exists in every known organism on Earth. SGC evolution via early unique codon assignment, then later wobble, yields coding resembling the near-universal code. Below, later wobble is shown to also create an optimal route to accurate codon assignment. Time of optimal codon assignment matches the previously defined mean time for ordered coding, exhibiting ≥ 90% of SGC order. Accurate evolution is also accessible, sufficiently frequent to appear in populations of 10³ to 10⁴ codes. SGC-like coding capacity, code order, and accurate assignments therefore arise together, in one attainable evolutionary intermediate. Examples, which plausibly resemble coding at evolutionary domain separation, are characterized.

Evolution of the Standard Genetic Code

A near-universal Standard Genetic Code (SGC) implies a single origin for present Earth life. To study this unique event, I compute paths to the SGC, comparing different plausible histories. Notably, SGC-like coding emerges from traditional evolutionary mechanisms, and a superior route can be identified. To objectively measure evolution, progress values from 0 (random coding) to 1 (SGC-like) are defined: these measure fractions of random-code-to-SGC distance. Progress types are spacing/distance/delta Polar Requirement, detecting space between identical assignments/mutational distance to the SGC/chemical order, respectively. The coding system is based on selected RNAs performing aminoacyl-RNA synthetase reactions. Acceptor RNAs exhibit SGC-like Crick wobble; alternatively, non-wobbling triplets uniquely encode 20 amino acids/start/stop. Triplets acquire 22 functions by stereochemistry, selection, coevolution, or at random. Assignments also propagate to an assigned triplet’s neighborhood via single mutations, but can also decay. A vast code universe makes futile evolutionary paths plentiful. Thus, SGC evolution is critically sensitive to disorder from random assignments. Evolution also inevitably slows near coding completion. The SGC likely avoided these difficulties, and two suitable paths are compared. In late wobble, a majority of non-wobble assignments are made before wobble is adopted. In continuous wobble, a uniquely advantageous early intermediate yields an ordered SGC. Revised coding evolution (limited randomness, late wobble, concentration on amino acid encoding, chemically conservative coevolution with a chemically ordered elite) produces varied full codes with excellent joint progress values. A population of only 600 independent coding tables includes SGC-like members; a Bayesian path toward more accurate SGC evolution is available.

Nonenzymatic Metabolic Reactions and Life's Origins

Prebiotic chemistry aims to explain how the biochemistry of life as we know it came to be. Most efforts in this area have focused on provisioning compounds of importance to life by multistep synthetic routes that do not resemble biochemistry. However, gaining insight into why core metabolism uses the molecules, reactions, pathways, and overall organization that it does requires us to consider molecules not only as synthetic end goals. Equally important are the dynamic processes that build them up and break them down. This perspective has led many researchers to the hypothesis that the first stage of the origin of life began with the onset of a primitive nonenzymatic version of metabolism, initially catalyzed by naturally occurring minerals and metal ions. This view of life's origins has come to be known as "metabolism first". Continuity with modern metabolism would require a primitive version of metabolism to build and break down ketoacids, sugars, amino acids, and ribonucleotides in much the same way as the pathways that do it today. This review discusses metabolic pathways of relevance to the origin of life in a manner accessible to chemists, and summarizes experiments suggesting several pathways might have their roots in prebiotic chemistry. Finally, key remaining milestones for the protometabolic hypothesis are highlighted.

Genetic code and metabolism: The perpetual waltz

Since the origin of life, metabolism and protein synthesis have evolved together to balance the vast amounts of ATP and amino acids required for genetic translation with the rest of the cell's energy needs. A new study offers satisfying insights into a long-standing evolutionary mystery surrounding a fused, bifunctional aminoacyl-tRNA synthetase. To avoid depleting cells from an essential amino acid generated by the Krebs cycle, harvesting for Glu and Pro by the translation machinery was unified in animals, thus preventing a Pro-hungry translational apparatus from depleting the cell of essential Glu reserves.

Genetic Code Evolution Investigated through the Synthesis and Characterisation of Proteins from Reduced-Alphabet Libraries

The universal genetic code of 20 amino acids is the product of evolution. It is believed that earlier versions of the code had fewer residues. Many theories for the order in which amino acids were integrated into the code have been proposed, considering factors ranging from prebiotic chemistry to codon capture. Several meta-analyses combined these theories to yield a feasible consensus chronology of the genetic code's evolution, but there is a dearth of experimental data to test the hypothesised order. We used combinatorial chemistry to synthesise libraries of random polypeptides that were based on different subsets of the 20 standard amino acids, thus representing different stages of a plausible history of the alphabet. Four libraries were comprised of the five, nine, and 16 most ancient amino acids, and all 20 extant residues for a direct side-by-side comparison. We characterised numerous variants from each library for their solubility and propensity to form secondary, tertiary or quaternary structures. Proteins from the two most ancient libraries were more likely to be soluble than those from the extant library. Several individual protein variants exhibited inducible protein folding and other traits typical of intrinsically disordered proteins. From these libraries, we can infer how primordial protein structure and function might have evolved with the genetic code.

Optimization of the standard genetic code according to three codon positions using an evolutionary algorithm

Many biological systems are typically examined from the point of view of adaptation to certain conditions or requirements. One such system is the standard genetic code (SGC), which generally minimizes the cost of amino acid replacements resulting from mutations or mistranslations. However, no full consensus has been reached on the factors that caused the evolution of this feature. One of the hypotheses suggests that code optimality was directly selected as an advantage to preserve information about encoded proteins. An important feature that should be considered when studying the SGC is the different roles of the three codon positions. Therefore, we investigated the robustness of this code regarding the cost of amino acid replacements resulting from substitutions in these positions separately and the sum of these costs. We applied a modified evolutionary algorithm and included four models of the genetic code assuming various restrictions on its structure. The SGC was compared both with the codes that minimize the objective function and those that maximize it. This approach allowed us to place the SGC in the global space of possible codes, which is a more appropriate and unbiased comparison than that with randomly generated codes because they are characterized by relatively uniform amino acid assignments to codons. The SGC appeared to be well optimized at the global scale, but its individual positions were not fully optimized because there were codes that were optimized for only one codon position and simultaneously outperformed the SGC at the other positions. We also found that different code structures may lead to the same optimality and that random codes can show a tendency to minimize costs under some of the genetic code models. Our results suggest that the optimality of SGC could be a by-product of other processes.

From the primordial soup to self-driving cars: standards and their role in natural and technological innovation

Standards are specifications to which the elements of a technology must conform. Here, we apply this notion to the biochemical 'technologies' of nature, where objects like DNA and proteins, as well as processes like the regulation of gene activity are highly standardized. We introduce the concept of standards with multiple examples, ranging from the ancient genetic material RNA, to Palaeolithic stone axes, and digital electronics, and we discuss common ways in which standards emerge in nature and technology. We then focus on the question of how standards can facilitate technological and biological innovation. Innovation-enhancing standards include those of proteins and digital electronics. They share common features, such as that few standardized building blocks can be combined through standard interfaces to create myriad useful objects or processes. We argue that such features will also characterize the most innovation-enhancing standards of future technologies.

An integrated, structure- and energy-based view of the genetic code

The principles of mRNA decoding are conserved among all extant life forms. We present an integrative view of all the interaction networks between mRNA, tRNA and rRNA: the intrinsic stability of codon-anticodon duplex, the conformation of the anticodon hairpin, the presence of modified nucleotides, the occurrence of non-Watson-Crick pairs in the codon-anticodon helix and the interactions with bases of rRNA at the A-site decoding site. We derive a more information-rich, alternative representation of the genetic code, that is circular with an unsymmetrical distribution of codons leading to a clear segregation between GC-rich 4-codon boxes and AU-rich 2:2-codon and 3:1-codon boxes. All tRNA sequence variations can be visualized, within an internal structural and energy framework, for each organism, and each anticodon of the sense codons. The multiplicity and complexity of nucleotide modifications at positions 34 and 37 of the anticodon loop segregate meaningfully, and correlate well with the necessity to stabilize AU-rich codon-anticodon pairs and to avoid miscoding in split codon boxes. The evolution and expansion of the genetic code is viewed as being originally based on GC content with progressive introduction of A/U together with tRNA modifications. The representation we present should help the engineering of the genetic code to include non-natural amino acids.

Origins of tmRNA: the missing link in the birth of protein synthesis?

The RNA world hypothesis refers to the early period on earth in which RNA was central in assuring both genetic continuity and catalysis. The end of this era coincided with the development of the genetic code and protein synthesis, symbolized by the apparition of the first non-random messenger RNA (mRNA). Modern transfer-messenger RNA (tmRNA) is a unique hybrid molecule which has the properties of both mRNA and transfer RNA (tRNA). It acts as a key molecule during trans-translation, a major quality control pathway of modern bacterial protein synthesis. tmRNA shares many common characteristics with ancestral RNA. Here, we present a model in which proto-tmRNAs were the first molecules on earth to support non-random protein synthesis, explaining the emergence of early genetic code. In this way, proto-tmRNA could be the missing link between the first mRNA and tRNA molecules and modern ribosome-mediated protein synthesis.

Analysis of the optimality of the standard genetic code

Many theories have been proposed attempting to explain the origin of the genetic code. In this work, we compare performance of the standard genetic code against millions of randomly generated codes. On left, ability of genetic codes to encode additional information and their robustness to frameshift mutations.

Ancestral AlaX editing enzymes for control of genetic code fidelity are not tRNA-specific

Accurate protein synthesis requires the hydrolytic editing of tRNAs incorrectly aminoacylated by aminoacyl-tRNA synthetases (ARSs). Recognition of cognate tRNAs by ARS is less error-prone than amino acid recognition, and, consequently, editing domains are generally believed to act only on the tRNAs cognate to their related ARSs. For example, the AlaX family of editing domains, including the editing domain of alanyl-tRNA synthetase and the related free-standing trans-editing AlaX enzymes, are thought to specifically act on tRNA(Ala), whereas the editing domains of threonyl-tRNA synthetases are specific for tRNA(Thr). Here we show that, contrary to this belief, AlaX-S, the smallest of the extant AlaX enzymes, deacylates Ser-tRNA(Thr) in addition to Ser-tRNA(Ala) and that a single residue is important to determine this behavior. Our data indicate that promiscuous forms of AlaX are ancestral to tRNA-specific AlaXs. We propose that former AlaX domains were used to maintain translational fidelity in earlier stages of genetic code evolution when mis-serylation of several tRNAs was possible.

Genetic code evolution started with the incorporation of glycine, followed by other small hydrophilic amino acids

We propose that glycine was the first amino acid to be incorporated into the genetic code, followed by serine, aspartic and/or glutamic acid-small hydrophilic amino acids that all have codons in the bottom right-hand corner of the standard genetic code table. Because primordial ribosomal synthesis is presumed to have been rudimentary, this stage would have been characterized by the synthesis of short, water-soluble peptides, the first of which would have comprised polyglycine. Evolution of the code is proposed to have occurred by the duplication and mutation of tRNA sequences, which produced a radiation of codon assignment outwards from the bottom right-hand corner. As a result of this expansion, we propose a trend from small hydrophilic to hydrophobic amino acids, with selection for longer polypeptides requiring a hydrophobic core for folding and stability driving the incorporation of hydrophobic amino acids into the code.

Coding nucleic acids are chaperons for protein folding: a novel theory of protein folding

The arguments for nucleic acid chaperons are reviewed and three new lines of evidence are added. (1) It was found that amino acids encoded by codons in short nucleic acid loops frequently form turns and helices in the corresponding protein structures. (2) The amino acids encoded by partially complementary (1st and 3rd nucleotides) codons are more frequently co-located in the encoded proteins than expected by chance. (3) There are significant correlations between thermodynamic changes (ddG) caused by codon mutations in nucleic acids and the thermodynamic changes caused by the corresponding amino acid mutations in the encoded proteins. We conclude that the concept of the Proteomic Code and nucleic acid chaperons seems correct from the bioinformatics point of view, and we expect to see direct biochemical experiments and evidence in the near future.

Imprints of the genetic code in the ribosome

The establishment of the genetic code remains elusive nearly five decades after the code was elucidated. The stereochemical hypothesis postulates that the code developed from interactions between nucleotides and amino acids, yet supporting evidence in a biological context is lacking. We show here that anticodons are selectively enriched near their respective amino acids in the ribosome, and that such enrichment is significantly correlated with the canonical code over random codes. Ribosomal anticodon-amino acid enrichment further reveals that specific codons were reassigned during code evolution, and that the code evolved through a two-stage transition from ancient amino acids without anticodon interaction to newer additions with anticodon interaction. The ribosome thus serves as a molecular fossil, preserving biological evidence that anticodon-amino acid interactions shaped the evolution of the genetic code.

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.

The close relationship between the biosynthetic families of amino acids and the organisation of the genetic code

By generating random codes and applying Fisher's exact test, we confirm that the biosynthetic families of amino acids are intimately involved in the organisation of the genetic code. This observation corroborates the coevolution theory of genetic code origin. As the amino acids belonging to the single biosynthetic families have codons that are contiguous in the genetic code, they must have entered the code itself by means of a clustering mechanism, which must clearly have been compatible with the mechanism on which this theory is based because this too envisages the clustering of biosynthetically correlated amino acids within the code.

An extension of the coevolution theory of the origin of the genetic code

BACKGROUND

The coevolution theory of the origin of the genetic code suggests that the genetic code is an imprint of the biosynthetic relationships between amino acids. However, this theory does not seem to attribute a role to the biosynthetic relationships between the earliest amino acids that evolved along the pathways of energetic metabolism. As a result, the coevolution theory is unable to clearly define the very earliest phases of genetic code origin. In order to remove this difficulty, I here suggest an extension of the coevolution theory that attributes a crucial role to the first amino acids that evolved along these biosynthetic pathways and to their biosynthetic relationships, even when defined by the non-amino acid molecules that are their precursors.

RESULTS

It is re-observed that the first amino acids to evolve along these biosynthetic pathways are predominantly those codified by codons of the type GNN, and this observation is found to be statistically significant. Furthermore, the close biosynthetic relationships between the sibling amino acids Ala-Ser, Ser-Gly, Asp-Glu, and Ala-Val are not random in the genetic code table and reinforce the hypothesis that the biosynthetic relationships between these six amino acids played a crucial role in defining the very earliest phases of genetic code origin.

CONCLUSION

All this leads to the hypothesis that there existed a code, GNS, reflecting the biosynthetic relationships between these six amino acids which, as it defines the very earliest phases of genetic code origin, removes the main difficulty of the coevolution theory. Furthermore, it is here discussed how this code might have naturally led to the code codifying only for the domains of the codons of precursor amino acids, as predicted by the coevolution theory. Finally, the hypothesis here suggested also removes other problems of the coevolution theory, such as the existence for certain pairs of amino acids with an unclear biosynthetic relationship between the precursor and product amino acids and the collocation of Ala between the amino acids Val and Leu belonging to the pyruvate biosynthetic family, which the coevolution theory considered as belonging to different biosyntheses.

REVIEWERS

This article was reviewed by Rob Knight, Paul Higgs (nominated by Laura Landweber), and Eugene Koonin.

Evolutionary patterns in the sequence and structure of transfer RNA: a window into early translation and the genetic code

Transfer RNA (tRNA) molecules play vital roles during protein synthesis. Their acceptor arms are aminoacylated with specific amino acid residues while their anticodons delimit codon specificity. The history of these two functions has been generally linked in evolutionary studies of the genetic code. However, these functions could have been differentially recruited as evolutionary signatures were left embedded in tRNA molecules. Here we built phylogenies derived from the sequence and structure of tRNA, we forced taxa into monophyletic groups using constraint analyses, tested competing evolutionary hypotheses, and generated timelines of amino acid charging and codon discovery. Charging of Sec, Tyr, Ser and Leu appeared ancient, while specificities related to Asn, Met, and Arg were derived. The timelines also uncovered an early role of the second and then first codon bases, identified codons for Ala and Pro as the most ancient, and revealed important evolutionary take-overs related to the loss of the long variable arm in tRNA. The lack of correlation between ancestries of amino acid charging and encoding indicated that the separate discoveries of these functions reflected independent histories of recruitment. These histories were probably curbed by co-options and important take-overs during early diversification of the living world.

Transfer RNA and the origins of diversified life

The evolution of the transfer RNA (tRNA) molecule is controversial but embeds the history of protein biosynthesis, the genetic code, and the origins of diversified life. A new phylogenetic method based on RNA structure that we developed provides new lines of evidence to support the genome tag hypothesis and confirms that the 'top half' of tRNA is more ancient than the 'bottom half'. Timelines of amino acid charging function generated from constraint analyses showed that selenocysteine, tyrosine, serine, and leucine specificities were ancient, while those related to asparagine, methionine, and arginine were more recent. The timelines also uncovered an early role of the second and then first codon bases, identified codons for alanine and proline as the most ancient, and revealed important evolutionary take-overs related to the loss of the long variable arm of tRNA. Furthermore, organismal timelines showed Archaea was the oldest superkingdom, followed by viruses, and superkingdoms Eukarya and Bacteria in that order supporting conclusions from recent phylogenomic studies of protein architecture. Strikingly, results showed that the origin of viruses was not only ancient but was linked to Archaea, supporting the notion that the archaeal lineage is the most ancient on earth and its origin predated diversification of tRNA function and specificity.

A content-centric organization of the genetic code

The codon table for the canonical genetic code can be rearranged in such a way that the code is divided into four quarters and two halves according to the variability of their GC and purine contents, respectively. For prokaryotic genomes, when the genomic GC content increases, their amino acid contents tend to be restricted to the GC-rich quarter and the purine-content insensitive half, where all codons are fourfold degenerate and relatively mutation-tolerant. Conversely, when the genomic GC content decreases, most of the codons retract to the AU-rich quarter and the purine-content sensitive half; most of the codons not only remain encoding physicochemically diversified amino acids but also vary when transversion (between purine and pyrimidine) happens. Amino acids with sixfold-degenerate codons are distributed into all four quarters and across the two halves; their fourfold-degenerate codons are all partitioned into the purine-insensitive half in favorite of robustness against mutations. The features manifested in the rearranged codon table explain most of the intrinsic relationship between protein coding sequences (the informational content) and amino acid compositions (the functional content). The renovated codon table is useful in predicting abundant amino acids and positioning the amino acids with related or distinct physicochemical properties.

Error Minimization and Coding Triplet/Binding Site Associations Are Independent Features of the Canonical Genetic Code

The canonical genetic code has been reported both to be error minimizing and to show stereochemical associations between coding triplets and binding sites. In order to test whether these two properties are unexpectedly overlapping, we generated 200,000 randomized genetic codes using each of five randomization schemes, with and without randomization of stop codons. Comparison of the code error (difference in polar requirement for single-nucleotide codon interchanges) with the coding triplet concentrations in RNA binding sites for eight amino acids shows that these properties are independent and uncorrelated. Thus, one is not the result of the other, and error minimization and triplet associations probably arose independently during the history of the genetic code. We explicitly show that prior fixation of a stereochemical core is consistent with an effective later minimization of error.

Origins of the genetic code: the escaped triplet theory

There is very significant evidence that cognate codons and/or anticodons are unexpectedly frequent in RNA-binding sites for seven of eight biological amino acids that have been tested. This suggests that a substantial fraction of the genetic code has a stereochemical basis, the triplets having escaped from their original function in amino acid-binding sites to become modern codons and anticodons. We explicitly show that this stereochemical basis is consistent with subsequent optimization of the code to minimize the effect of coding mistakes on protein structure. These data also strengthen the argument for invention of the genetic code in an RNA world and for the RNA world itself.

Designing a neural network for the constraint optimization of the fitness functions devised based on the load minimization of the genetic code

Nonrandom patterns in codon assignments are supported by many statistical and biochemical studies in the last two decades. The canonical genetic code is known to be highly efficient in minimizing the effects of mistranslational errors and point mutations, an ability, which in term is designated "load minimization". Prior studies have included many attempts at quantitative estimation of the fraction of randomly generated codes, which in terms of load minimization, score higher than the canonical genetic code. In this study, a neural network, which estimates a highly optimized genetic code in a relatively short period of time has been devised. Several fitness functions were used throughout this text. Meanwhile, we have made use of two cost measure matrices, PAM74-100 and mutation matrix.

On the optimality of the genetic code, with the consideration of coevolution theory by comparison of prominent cost measure matrices

Statistical and biochemical studies have revealed non-random patterns in codon assignments. The canonical genetic code is known to be highly efficient in minimizing the effects of mistranslation errors and point mutations, since it is known that when an amino acid is converted to another due to error, the biochemical properties of the resulted amino acid are usually very similar to those of the original one. In this study, using altered forms of the fitness functions used in the prior studies, we have optimized the parameters involved in the calculation of the error minimizing property of the genetic code so that the genetic code outscores the random codes as much as possible. This work also compares two prominent matrices, the Mutation Matrix and Point Accepted Mutations 74-100 (PAM(74-100)). It has been resulted that the hypothetical properties of the coevolution theory of the genetic code are already considered in PAM(74-100), giving more evidence on the existence of bias towards the genetic code in this matrix. Furthermore, our results indicate that PAM(74-100) is biased towards the single base mistranslation occurrences in second codon position as well as the frequency of amino acids. Thus PAM(74-100) is not a suitable substitution matrix for the studies conducted on the evolution of the genetic code.

On the structural regularity in nucleobases and amino acids and relationship to the origin and evolution of the genetic code

To explore how chemical structures of both nucleobases and amino acids may have played a role in shaping the genetic code, numbers of sp2 hybrid nitrogen atoms in nucleobases were taken as a determinative measure for empirical stereo-electronic property to analyze the genetic code. Results revealed that amino acid hydropathy correlates strongly with the sp2 nitrogen atom numbers in nucleobases rather than with the overall electronic property such as redox potentials of the bases, reflecting that stereo-electronic property of bases may play a role. In the rearranged code, five simple but stereo-structurally distinctive amino acids (Gly, Pro, Val, Thr and Ala) and their codon quartets form a crossed intersection "core". Secondly, a re-categorization of the amino acids according to their beta-carbon stereochemistry, verified by charge density (at beta-carbon) calculation, results in five groups of stereo-structurally distinctive amino acids, the group leaders of which are Gly, Pro, Val, Thr and Ala, remarkably overlapping the above "core". These two lines of independent observations provide empirical arguments for a contention that a seemingly "frozen" "core" could have formed at a certain evolutionary stage. The possible existence of this codon "core" is in conformity with a previous evolutionary model whereby stereochemical interactions may have shaped the code. Moreover, the genetic code listed in UCGA succession together with this codon "core" has recently facilitated an identification of the unprecedented icosikaioctagon symmetry and bi-pyramidal nature of the genetic code.

A new classification scheme of the genetic code

Since the early days of the discovery of the genetic code nonrandom patterns have been searched for in the code in the hope of providing information about its origin and early evolution. Here we present a new classification scheme of the genetic code that is based on a binary representation of the purines and pyrimidines. This scheme reveals known patterns more clearly than the common one, for instance, the classification of strong, mixed, and weak codons as well as the ordering of codon families. Furthermore, new patterns have been found that have not been described before: Nearly all quantitative amino acid properties, such as Woese's polarity and the specific volume, show a perfect correlation to Lagerkvist's codon-anticodon binding strength. Our new scheme leads to new ideas about the evolution of the genetic code. It is hypothesized that it started with a binary doublet code and developed via a quaternary doublet code into the contemporary triplet code. Furthermore, arguments are presented against suggestions that a "simpler" code, where only the midbase was informational, was at the origin of the genetic code.

Expanding the Genetic Code

Although chemists can synthesize virtually any small organic molecule, our ability to rationally manipulate the structures of proteins is quite limited, despite their involvement in virtually every life process. For most proteins, modifications are largely restricted to substitutions among the common 20 amino acids. Herein we describe recent advances that make it possible to add new building blocks to the genetic codes of both prokaryotic and eukaryotic organisms. Over 30 novel amino acids have been genetically encoded in response to unique triplet and quadruplet codons including fluorescent, photoreactive, and redox-active amino acids, glycosylated amino acids, and amino acids with keto, azido, acetylenic, and heavy-atom-containing side chains. By removing the limitations imposed by the existing 20 amino acid code, it should be possible to generate proteins and perhaps entire organisms with new or enhanced properties.

Coevolution theory of the genetic code at age thirty

The coevolution theory of the genetic code, which postulates that prebiotic synthesis was an inadequate source of all twenty protein amino acids, and therefore some of them had to be derived from the coevolving pathways of amino acid biosynthesis, has been assessed in the light of the discoveries of the past three decades. Its four fundamental tenets regarding the essentiality of amino acid biosynthesis, role of pretran synthesis, biosynthetic imprint on codon allocations and mutability of the encoded amino acids are proven by the new knowledge. Of the factors that guided the evolutionary selection of the universal code, the relative contributions of Amino Acid Biosynthesis: Error Minimization: Stereochemical Interaction are estimated to first approximation as 40,000,000:400:1, which suggests that amino acid biosynthesis represents the dominant factor shaping the code. The utility of the coevolution theory is demonstrated by its opening up experimental expansions of the code and providing a basis for locating the root of life.

On the optimality of the genetic code, with the consideration of termination codons

The existence of nonrandom patterns in codon assignments is supported by many statistical and biochemical studies. The canonical genetic code is known to be highly efficient in minimizing the effects of mistranslation errors and point mutations. For example, it is known that when an error induces the conversion of an amino acid to another, the biochemical properties of the resulting amino acid are usually very similar to that of the original. Prior studies include many attempts at quantitative estimation of the fraction of randomly generated codes which, based upon load minimization, score higher than the canonical genetic code. In this study, we took into consideration both the relative frequencies of amino acids and nonsense mistranslations, factors which had been previously ignored. Incorporation of these parameters, resulted in a fitness function (phi) which rendered the canonical genetic code to be highly optimized with respect to load minimization. Considering termination codons, we applied a biosynthetic version of the coevolution theory, however, with low significance. We employed a revised cost for the precursor-product pairs of amino acids and showed that the significance of this approach depends on the cost measure matrix used by the researcher. Thus, we have compared the two prominent matrices, point accepted mutations 74-100 (PAM(74-100)) and mutation matrix in our study.

The case for an error minimizing standard genetic code

Since discovering the pattern by which amino acids are assigned to codons within the standard genetic code, investigators have explored the idea that natural selection placed biochemically similar amino acids near to one another in coding space so as to minimize the impact of mutations and/or mistranslations. The analytical evidence to support this theory has grown in sophistication and strength over the years, and counterclaims questioning its plausibility and quantitative support have yet to transcend some significant weaknesses in their approach. These weaknesses are illustrated here by means of a simple simulation model for adaptive genetic code evolution. There remain ill explored facets of the 'error minimizing' code hypothesis, however, including the mechanism and pathway by which an adaptive pattern of codon assignments emerged, the extent to which natural selection created synonym redundancy, its role in shaping the amino acid and nucleotide languages, and even the correct interpretation of the adaptive codon assignment pattern: these represent fertile areas for future research.

On the evolution of primitive genetic codes

The primordial genetic code probably has been a drastically simplified ancestor of the canonical code that is used by contemporary cells. In order to understand how the present-day code came about we first need to explain how the language of the building plan can change without destroying the encoded information. In this work we introduce a minimal organism model that is based on biophysically reasonable descriptions of RNA and protein, namely secondary structure folding and knowledge based potentials. The evolution of a population of such organism under competition for a common resource is simulated explicitly at the level of individual replication events. Starting with very simple codes, and hence greatly reduced amino acid alphabets, we observe a diversification of the codes in most simulation runs. The driving force behind this effect is the possibility to produce fitter proteins when the repertoire of amino acids is enlarged.

A common periodic table of codons and amino acids

A periodic table of codons has been designed where the codons are in regular locations. The table has four fields (16 places in each) one with each of the four nucleotides (A, U, G, C) in the central codon position. Thus, AAA (lysine), UUU (phenylalanine), GGG (glycine), and CCC (proline) were placed into the corners of the fields as the main codons (and amino acids) of the fields. They were connected to each other by six axes. The resulting nucleic acid periodic table showed perfect axial symmetry for codons. The corresponding amino acid table also displaced periodicity regarding the biochemical properties (charge and hydropathy) of the 20 amino acids and the position of the stop signals. The table emphasizes the importance of the central nucleotide in the codons and predicts that purines control the charge while pyrimidines determine the polarity of the amino acids. This prediction was experimentally tested.

Generation of a bacterium with a 21 amino acid genetic code

We have generated a completely autonomous bacterium with a 21 amino acid genetic code. This bacterium can biosynthesize a nonstandard amino acid from basic carbon sources and incorporate this amino acid into proteins in response to the amber nonsense codon. The biosynthetic pathway for the amino acid p-aminophenylalanine (pAF) as well as a unique pAF synthetase and cognate tRNA were added to Escherichia coli. Denaturing gel electrophoresis and mass spectrometric analysis show that pAF is incorporated into myoglobin with fidelity and efficiency rivaling those of the common 20 amino acids. This and other such organisms may provide an opportunity to examine the evolutionary consequences of adding new amino acids to the genetic repertoire, as well as generate proteins with new or enhanced biological functions.

Molecular evolution before the origin of species

Amino acids at conserved sites in the residue sequence of 10 ancient proteins, from 844 phylogenetically diverse sources, were used to specify their time of origin in the interval before species divergence from the last common ancestor (LCA). The order of amino acid addition to the genetic code, based on biosynthesis path length and other molecular evidence, provided a reference for evaluating the 'code age' of each residue profile examined. Significantly earlier estimates were obtained for conserved amino acid residues in these proteins than non-conserved residues. Evidence from the primary structure of 'fossil' proteins thus corroborated the biosynthetic order of amino acid addition to the code.Low potential ferredoxin (Fdxn) had the earliest residue profile among the proteins in this study. A phylogenetic tree for 82 prokaryote Fdxn sequences was rooted midway between bacteria and archaea branches. LCA Fdxn had a 23-residue antecedent whose residue profile matched mid-expansion phase codon assignments and included an amide residue. It contained a highly acidic N-terminal region and a non-charged C-terminal region, with all four cysteine residues. This small protein apparently anchored a [4Fe-4S] cluster, ligated by C-terminal cysteines, to a positively charged mineral surface, consistent with mediating e(-) transfer in a primordial surface system before cells appeared. Its negatively charged N-terminal 'attachment site' was highly mutable during evolution of ancestral Fdxn for Bacteria and Archaea, consistent with a loss of function after cell formation. An initial glutamate to lysine substitution may link 'attachment site' removal to early post-expansion phase entry of basic amino acids to the code. As proteins evidently anchored non-charged amide residues initially, surface attachment of cofactors and other functional groups emerges as a general function of pre-cell proteins.A phylogenetic tree of 107 proteolipid (PL) helix-1 sequences from H(+)-ATPase of bacteria, archaea and eukaryotes had its root between prokaryote branches. LCA PL h1 residue profile optimally fit a late expansion phase codon array. Sequence repeats in transmembrane PL helices h1 and h2 indicated formation of the archetypal PL hairpin structure involved successive tandem duplications, initiated within the gene for an 11-residue (or 4-residue) hydrophobic peptide. Ancestral PL h1 lacked acidic residues, in a fundamental departure from the prototype pre-cell protein. By this stage, proteins with a hydrophobic domain had evolved. Its non-polar, late expansion phase residue profile point to ancestral PL being a component of an early permeable cell membrane. Other indicators of cell formation about this stage of code evolution include phospholipid biosynthesis path length, FtsZ residue profile, and late entry of basic amino acids into the genetic code. Estimates based on conserved residues in prokaryote cell septation protein, FtsZ, and proteins involved with synthesis, transcription and replication of DNA revealed FtsZ, ribonucleotide reductase, RNA polymerase core subunits and 5'-->3' flap exonuclease, FEN-1, originated soon after cells putatively evolved. While reverse transcriptase and topoisomerase I, Topo I, appeared late in the pre-divergence era, when the genetic code was essentially complete. The transition from RNA genes to a DNA genome seemingly proceeded via formation of a DNA-RNA heteroduplex. These results suggest formation of DNA awaited evolution of a catalyst with a hydrophobic domain, capable of sequestering radical bearing intermediates in its synthesis from ribonucleotide precursors. Late formation of topology altering protein, Topo I, further suggests consolidation of genes into chromosomes followed synthesis of comparatively thermostable DNA strands.

Optimality of the genetic code with respect to protein stability and amino-acid frequencies

BACKGROUND

The genetic code is known to be efficient in limiting the effect of mistranslation errors. A misread codon often codes for the same amino acid or one with similar biochemical properties, so the structure and function of the coded protein remain relatively unaltered. Previous studies have attempted to address this question quantitatively, by estimating the fraction of randomly generated codes that do better than the genetic code in respect of overall robustness. We extended these results by investigating the role of amino-acid frequencies in the optimality of the genetic code.

RESULTS

We found that taking the amino-acid frequency into account decreases the fraction of random codes that beat the natural code. This effect is particularly pronounced when more refined measures of the amino-acid substitution cost are used than hydrophobicity. To show this, we devised a new cost function by evaluating in silico the change in folding free energy caused by all possible point mutations in a set of protein structures. With this function, which measures protein stability while being unrelated to the code's structure, we estimated that around two random codes in a billion (109) are fitter than the natural code. When alternative codes are restricted to those that interchange biosynthetically related amino acids, the genetic code appears even more optimal.

CONCLUSIONS

These results lead us to discuss the role of amino-acid frequencies and other parameters in the genetic code's evolution, in an attempt to propose a tentative picture of primitive life.

Computer-modeling origin of a simple genetic apparatus

This computer simulation is based on a model of the origin of life proposed by H. Kuhn and J. Waser, where the evolution of short molecular strands is assumed to take place in a distinct spatiotemporal structured environment. In their model, the prebiotic situation is strongly simplified to grasp essential features of the evolution of the genetic apparatus without attempts to trace the historic path. With the tool of computer implementation confining to principle aspects and focused on critical features of the model, a deeper understanding of the model's premises is achieved. Each generation consists of three steps: (i) construction of devices (entities exposed to selection) presently available; (ii) selection; and (iii) multiplication of the isolated strands (R oligomers) by complementary copying with occasional variation by copying mismatch. In the beginning, the devices are single strands with random sequences; later, increasingly complex aggregates of strands form devices such as a hairpin-assembler device which develop in favorable cases. A monomers interlink by binding to the hairpin-assembler device, and a translation machinery, called the hairpin-assembler-enzyme device, emerges, which translates the sequence of R(1) and R(2) monomers in the assembler strand to the sequence of A(1) and A(2) monomers in the A oligomer, working as an enzyme.