Ancestral Insertions and Expansions of rRNA do not Support an Origin of the Ribosome in Its Peptidyl Transferase Center

Origin of life: Drawing the big picture

Trying to provide a broad overview about the origin of life in Earth, the most significant transitions of life before cells are listed and discussed. The current approach emphasizes the symbiotic relationships that emerged with life. We propose a rational, stepwise scenario for the origin of life that starts with the origin of the first biomolecules and steps forward until the origins of the first cells. Along this path, we aim to provide a brief, though comprehensive theoretical model that will consider the following steps: (i) how nucleotides and other biomolecules could be made prebiotically in specific prebiotic refuges; (ii) how the first molecules of RNAs were formed; (iii) how the proto-peptidyl transferase center was built by the concatenation of proto-tRNAs; (iv) how the ribosome and the genetic code could be structured; (v) how progenotes could live and reproduce as "naked" ribonucleoprotein molecules; (vi) how peptides started to bind molecules in the prebiotic soup allowing biochemical pathways to evolve from those bindings; (vii) how genomes got bigger by the symbiotic relationship of progenotes and lateral transference of genetic material; (viii) how the progenote LUCA has been formed by assembling most biochemical routes; (ix) how the first virion capsids probably emerged and evolved; (x) how phospholipid membranes emerged probably twice by the evolution of lipid-binding proteins; (xi) how DNA synthesis have been formed in parallel in Bacteria and Archaea; and, finally, (xii) how DNA-based cells of Bacteria and Archaeabacteria have been constituted. The picture provided is conjectural and present epistemological gaps. Future research will help to advance into the elucidation of gaps and confirmation/refutation of current statements.

The story of rRNA expansion segments: Finding functionality amidst diversity

Expansion segments (ESs) are multinucleotide insertions present across phyla at specific conserved positions in eukaryotic rRNAs. ESs are generally absent in bacterial rRNAs with some exceptions, while the archaeal rRNAs have microexpansions at regions that coincide with those of eukaryotic ESs. Although there is an increasing prominence of ribosomes, especially the ribosomal proteins, in fine-tuning gene expression through translation regulation, the role of rRNA ESs is relatively underexplored. While rRNAs have been established as the major catalytic hub in ribosome function, the presence of ESs widens their scope as a species-specific regulatory hub of protein synthesis. In this comprehensive review, we have elaborately discussed the current understanding of the functional aspects of rRNA ESs of cytoplasmic eukaryotic ribosomes and discuss their past, present, and future. This article is categorized under: RNA Structure and Dynamics > Influence of RNA Structure in Biological Systems Translation > Ribosome Structure/Function Translation > Regulation.

A Short Tale of the Origin of Proteins and Ribosome Evolution

Proteins are the workhorses of the cell and have been key players throughout the evolution of all organisms, from the origin of life to the present era. How might life have originated from the prebiotic chemistry of early Earth? This is one of the most intriguing unsolved questions in biology. Currently, however, it is generally accepted that amino acids, the building blocks of proteins, were abiotically available on primitive Earth, which would have made the formation of early peptides in a similar fashion possible. Peptides are likely to have coevolved with ancestral forms of RNA. The ribosome is the most evident product of this coevolution process, a sophisticated nanomachine that performs the synthesis of proteins codified in genomes. In this general review, we explore the evolution of proteins from their peptide origins to their folding and regulation based on the example of superoxide dismutase (SOD1), a key enzyme in oxygen metabolism on modern Earth.

Codon assignment evolvability in theoretical minimal RNA rings

Genetic code codon-amino acid assignments evolve for 15 (AAA, AGA, AGG, ATA, CGG, CTA, CTG. CTC, CTT, TAA, TAG, TCA, TCG, TGA and TTA (GNN codons notably absent)) among 64 codons (23.4%) across the 31 genetic codes (NCBI list completed with recently suggested green algal mitochondrial genetic codes). Their usage in 25 theoretical minimal RNA rings is examined. RNA rings are designed in silico to code once over the shortest length for all 22 coding signals (start and stop codons and each amino acid according to the standard genetic code). Though designed along coding constraints, RNA rings resemble ancestral tRNA loops, assigning to each RNA ring a putative anticodon, a cognate amino acid and an evolutionary genetic code integration rank for that cognate amino acid. Analyses here show 1. biases against/for evolvable codons in the two first vs last thirds of RNA ring coding sequences, 2. RNA rings with evolvable codons have recent cognates, and 3. evolvable codon and cytosine numbers in RNA ring compositions are positively correlated. Applying alternative genetic codes to RNA rings designed for nonredundant coding according to the standard genetic code reveals unsuspected properties of the standard genetic code and of RNA rings, notably on codon assignment evolvability and the special role of cytosine in relation to codon assignment evolvability and of the genetic code's coding structure.

José M. The Ancient History of Peptidyl Transferase Center Formation as Told by Conservation and Information Analyses

The peptidyl transferase center (PTC) is the catalytic center of the ribosome and forms part of the 23S ribosomal RNA. The PTC has been recognized as the earliest ribosomal part and its origins embodied the First Universal Common Ancestor (FUCA). The PTC is frequently assumed to be highly conserved along all living beings. In this work, we posed the following questions: (i) How many 100% conserved bases can be found in the PTC? (ii) Is it possible to identify clusters of informationally linked nucleotides along its sequence? (iii) Can we propose how the PTC was formed? (iv) How does sequence conservation reflect on the secondary and tertiary structures of the PTC? Aiming to answer these questions, all available complete sequences of 23S ribosomal RNA from Bacteria and Archaea deposited on GenBank database were downloaded. Using a sequence bait of 179 bp from the PTC of Thermus termophilus, we performed an optimum pairwise alignment to retrieve the PTC region from 1424 filtered 23S rRNA sequences. These PTC sequences were multiply aligned, and the conserved regions were assigned and observed along the primary, secondary, and tertiary structures. The PTC structure was observed to be more highly conserved close to the adenine located at the catalytical site. Clusters of interrelated, co-evolving nucleotides reinforce previous assumptions that the PTC was formed by the concatenation of proto-tRNAs and important residues responsible for its assembly were identified. The observed sequence variation does not seem to significantly affect the 3D structure of the PTC ribozyme.

Comparisons between small ribosomal RNA and theoretical minimal RNA ring secondary structures confirm phylogenetic and structural accretion histories

Ribosomal RNAs are complex structures that presumably evolved by tRNA accretions. Statistical properties of tRNA secondary structures correlate with genetic code integration orders of their cognate amino acids. Ribosomal RNA secondary structures resemble those of tRNAs with recent cognates. Hence, rRNAs presumably evolved from ancestral tRNAs. Here, analyses compare secondary structure subcomponents of small ribosomal RNA subunits with secondary structures of theoretical minimal RNA rings, presumed proto-tRNAs. Two independent methods determined different accretion orders of rRNA structural subelements: (a) classical comparative homology and phylogenetic reconstruction, and (b) a structural hypothesis assuming an inverted onion ring growth where the three-dimensional ribosome's core is most ancient and peripheral elements most recent. Comparisons between (a) and (b) accretions orders with RNA ring secondary structure scales show that recent rRNA subelements are: 1. more like RNA rings with recent cognates, indicating ongoing coevolution between tRNA and rRNA secondary structures; 2. less similar to theoretical minimal RNA rings with ancient cognates. Our method fits (a) and (b) in all examined organisms, more with (a) than (b). Results stress the need to integrate independent methods. Theoretical minimal RNA rings are potential evolutionary references for any sequence-based evolutionary analyses, independent of the focal data from that study.

Identification of a circular code periodicity in the bacterial ribosome: origin of codon periodicity in genes?

Three-base periodicity (TBP), where nucleotides and higher order n-tuples are preferentially spaced by 3, 6, 9, etc. bases, is a well-known intrinsic property of protein-coding DNA sequences. However, its origins are still not fully understood. One hypothesis is that the periodicity reflects a primordial coding system that was used before the emergence of the modern standard genetic code (SGC). Recent evidence suggests that the X circular code, a set of 20 trinucleotides allowing the reading frames in genes to be retrieved locally, represents a possible ancestor of the SGC. Motifs from the X circular code have been found in the reading frame of protein-coding regions in extant organisms from bacteria to eukaryotes, in many transfer RNA (tRNA) genes and in important functional regions of the ribosomal RNA (rRNA), notably in the peptidyl transferase centre and the decoding centre. Here, we have used a powerful correlation function to search for periodicity patterns involving the 20 trinucleotides of the X circular code in a large set of bacterial protein-coding genes, as well as in the translation machinery, including rRNA and tRNA sequences. As might be expected, we found a strong circular code periodicity 0 modulo 3 in the protein-coding genes. More surprisingly, we also identified a similar circular code periodicity in a large region of the 16S rRNA. This region includes the 3' major domain corresponding to the primordial proto-ribosome decoding centre and containing numerous sites that interact with the tRNA and messenger RNA (mRNA) during translation. Furthermore, 3D structural analysis shows that the periodicity region surrounds the mRNA channel that lies between the head and the body of the SSU. Our results support the hypothesis that the X circular code may constitute an ancestral translation code involved in reading frame retrieval and maintenance, traces of which persist in modern mRNA, tRNA and rRNA despite their long evolution and adaptation to the SGC.

Deamination gradients within codons after 1<->2 position swap predict amino acid hydrophobicity and parallel β-sheet conformational preference

Deaminations C->T and A->G are frequent mutations producing nucleotide content gradients across genomes proportional to singlestrandedness during replication/transcription. Hence, within single codons, deamination risks increase from first to third codon positions, while second codon positions are functionally most crucial. Here genetic codes are analyzed assuming that after anticodons protected codons from deaminations, first and second codon positions swapped (N2N1N3->N1N2N3), with lowest deamination risks for N2 in presumed primitive N2N1N3 codons. N2N1N3, not standard N1N2N3, codon structure minimizes deaminations inversely proportionally to cognate amino acid hydrophobicity and parallel betasheet conformational preference. For N1N2N3, deamination minimization increases with genetic code integration order of cognate amino acids: during the presumed N2N1N3->N1N2N3 codon structure transition, protein synthesis combined direct codon-amino acid interactions for late amino acids and tRNA-based translation for early amino acids. Hence N2N1N3 codons would correspond to tRNA-free translation by spontaneous codon-amino acid affinities, and tRNA-mediated translation presumably caused N2N1N3->N1N2N3 swaps. Results show that rational, not arbitrary rules link codon and amino acid structures. Some analyses detect mitochondrial RNAs and peptides in public data corresponding to systematic position swaps, suggesting occasional swapping polymerase activity.

Accretion history of large ribosomal subunits deduced from theoretical minimal RNA rings is congruent with histories derived from phylogenetic and structural methods

Accretions of tRNAs presumably formed the large complex ribosomal RNA structures. Similarities of tRNA secondary structures with rRNA secondary structures increase with the integration order of their cognate amino acid in the genetic code, indicating tRNA evolution towards rRNA-like structures. Here analyses rank secondary structure subelements of three large ribosomal RNAs (Prokaryota: Archaea: Thermus thermophilus; Bacteria: Escherichia coli; Eukaryota: Saccharomyces cerevisiae) in relation to their similarities with secondary structures formed by presumed proto-tRNAs, represented by 25 theoretical minimal RNA rings. These ranks are compared to those derived from two independent methods (ranks provide a relative evolutionary age to the rRNA substructure), (a) cladistic phylogenetic analyses and (b) 3D-crystallography where core subelements are presumed ancient and peripheral ones recent. Comparisons of rRNA secondary structure subelements with RNA ring secondary structures show congruence between ranks deduced by this method and both (a) and (b) (more with (a) than (b)), especially for RNA rings with predicted ancient cognate amino acid. Reconstruction of accretion histories of large rRNAs will gain from adequately integrating information from independent methods. Theoretical minimal RNA rings, sequences deterministically designed in silico according to specific coding constraints, might produce adequate scales for prebiotic and early life molecular evolution.

Possible Emergence of Sequence Specific RNA Aminoacylation via Peptide Intermediary to Initiate Darwinian Evolution and Code Through Origin of Life

One of the most intriguing questions in biological science is how life originated on Earth. A large number of hypotheses have been proposed to explain it, each putting an emphasis on different events leading to functional translation and self-sustained system. Here, we propose a set of interactions that could have taken place in the prebiotic environment. According to our hypothesis, hybridization-induced proximity of short aminoacylated RNAs led to the synthesis of peptides of random sequence. We postulate that among these emerged a type of peptide(s) capable of stimulating the interaction between specific RNAs and specific amino acids, which we call "bridge peptide" (BP). We conclude that translation should have emerged at the same time when the standard genetic code begun to evolve due to the stabilizing effect on RNA-peptide complexes with the help of BPs. Ribosomes, ribozymes, and the enzyme-directed RNA replication could co-evolve within the same period, as logical outcome of RNA-peptide world without the need of RNA only self-sustained step.

Coding of Class I and II Aminoacyl-tRNA Synthetases

The aminoacyl-tRNA synthetases and their cognate transfer RNAs translate the universal genetic code. The twenty canonical amino acids are sufficiently diverse to create a selective advantage for dividing amino acid activation between two distinct, apparently unrelated superfamilies of synthetases, Class I amino acids being generally larger and less polar, Class II amino acids smaller and more polar. Biochemical, bioinformatic, and protein engineering experiments support the hypothesis that the two Classes descended from opposite strands of the same ancestral gene. Parallel experimental deconstructions of Class I and II synthetases reveal parallel losses in catalytic proficiency at two novel modular levels-protozymes and Urzymes-associated with the evolution of catalytic activity. Bi-directional coding supports an important unification of the proteome; affords a genetic relatedness metric-middle base-pairing frequencies in sense/antisense alignments-that probes more deeply into the evolutionary history of translation than do single multiple sequence alignments; and has facilitated the analysis of hitherto unknown coding relationships in tRNA sequences. Reconstruction of native synthetases by modular thermodynamic cycles facilitated by domain engineering emphasizes the subtlety associated with achieving high specificity, shedding new light on allosteric relationships in contemporary synthetases. Synthetase Urzyme structural biology suggests that they are catalytically-active molten globules, broadening the potential manifold of polypeptide catalysts accessible to primitive genetic coding and motivating revisions of the origins of catalysis. Finally, bi-directional genetic coding of some of the oldest genes in the proteome places major limitations on the likelihood that any RNA World preceded the origins of coded proteins.

The progene hypothesis: the nucleoprotein world and how life began

In this article, I review the results of studies on the origin of life distinct from the popular RNA world hypothesis. The alternate scenario postulates the origin of the first bimolecular genetic system (a polynucleotide gene and a polypeptide processive polymerase) with simultaneous replication and translation and includes the following key features: 1. The bimolecular genetic system emerges not from mononucleotides and monoamino acids, but from progenes, namely, trinucleotides aminoacylated on 3'-end by a non-random amino acid (NpNpNp ~ pX ~ Aa, where N--deoxyribo- or ribonucleoside, p--phosphate, X--a bifunctional agent, for example ribose, Aa--amino acid, ~ macroerge bond). Progenes are used as substrates for simultaneous synthesis of a polynucleotide and a polypeptide. Growth of the system is controlled by the growing polypeptide, and the bimolecular genetic system emerges as an extremely rare event. The first living being (virus-like organism protoviroid, Protoviroidum primum) arises and reproduces in prebiotic liposome-like structures using progenes. A population of protoviroids possessing the genetic system evolves in accordance with the Darwinian principle. Early evolution from protoviroid world to protocell world is shortly described. 2. The progene forming mechanism (NpNp + Np ~ pX ~ Aa) makes it possible to explain the emergence of the prebiotic physicochemical group genetic code, as well as the selection of organic compounds for the future genetic system from the racemic environment. 3. The protoviroid is reproduced on a progene basis via replicative transcription-translation (RTT, the first molecular genetic process) that is similar to its modern counterparts. Nothing is required for the emergence and reproduction of the protoviroid except for progenes and conditions for their formation. 4. The general scheme of early evolution is as follows: prebiotic world → protoviroid (nucleoprotein) world → protocell (DNA-RNA-protein) world → LUCA (Last Universal Common Ancestor) → modern cell world. This scheme exclude the existence of an independent RNA world as predecessor of the cellular world.

Computing the origin and evolution of the ribosome from its structure - Uncovering processes of macromolecular accretion benefiting synthetic biology

Accretion occurs pervasively in nature at widely different timeframes. The process also manifests in the evolution of macromolecules. Here we review recent computational and structural biology studies of evolutionary accretion that make use of the ideographic (historical, retrodictive) and nomothetic (universal, predictive) scientific frameworks. Computational studies uncover explicit timelines of accretion of structural parts in molecular repertoires and molecules. Phylogenetic trees of protein structural domains and proteomes and their molecular functions were built from a genomic census of millions of encoded proteins and associated terminal Gene Ontology terms. Trees reveal a ‘metabolic-first’ origin of proteins, the late development of translation, and a patchwork distribution of proteins in biological networks mediated by molecular recruitment. Similarly, the natural history of ancient RNA molecules inferred from trees of molecular substructures built from a census of molecular features shows patchwork-like accretion patterns. Ideographic analyses of ribosomal history uncover the early appearance of structures supporting mRNA decoding and tRNA translocation, the coevolution of ribosomal proteins and RNA, and a first evolutionary transition that brings ribosomal subunits together into a processive protein biosynthetic complex. Nomothetic structural biology studies of tertiary interactions and ancient insertions in rRNA complement these findings, once concentric layering assumptions are removed. Patterns of coaxial helical stacking reveal a frustrated dynamics of outward and inward ribosomal growth possibly mediated by structural grafting. The early rise of the ribosomal ‘turnstile’ suggests an evolutionary transition in natural biological computation. Results make explicit the need to understand processes of molecular growth and information transfer of macromolecules.

Ribosomal accretion, apriorism and the phylogenetic method: a response to Petrov and Williams

Historical (ideographic) and non-historical (nomothetic) studies of ribosomal accretion appear to arrive at diametrically opposite conclusions. Phylogenetic analysis of thousands of RNA molecules and protein structures in hundreds of genomes supports the structural origin of the ribosome in RNA decoding and ribosomal mechanics. Predictions from extant features in a handful of rRNA structural models of the large ribosomal subunit support its origin in protein biosynthesis. In recent correspondence, one of us reported that correcting dismissals of conflicting data and avoiding unwarranted assumptions of the nomothetic method reconciled conclusions. In response, Petrov and Williams dismissed our arguments claiming we did not understand their algorithmic model of ribosomal apical growth. Instead, they controverted the historical approach. Here we show that their objections to the phylogenetic method are unjustified, that their algorithm subjectively guarantees back-in-time molecular deconstructions toward the protein biosynthetic core, and that processes of ribosomal growth are much more complex. We prompt abandoning apriorism, decreasing ad hoc hypotheses and integrating historical and non-historical scientific methods.

The ancient heart of the ribosomal large subunit: a response to Caetano-Anolles

Our recent Accretion Model of ribosomal evolution uses insertion fingerprints and a "trunk-branch" formalism to recapitulate the building up of common core rRNA of the Large Ribosomal Subunit. The Accretion Model is a conservative and natural extension of a method developed by Bokov and Steinberg (Nature 457:977-80, 2009), which confirms the correctness of lower resolution models by Fox and others. In each of these models, the LSU originates with the peptidyl transferase center (PTC), consistent with expectations that the ribosome is the source of defined-sequence functional proteins. In an adjacent note, Caetano-Anolles (J Mol Evol 80:162-165, 2015) disparages the Accretion Model, because it controverts the 'Growth Inferred by Genothermal Ordering' (GIGO) model. GIGO analyzes secondary structures, assigns the origin of the ribosome to a region outside of the PTC, and assumes or deduces that (i) large protein enzymes of defined amino acid sequence predate ribosomal synthesis of proteins, (ii) proteins directly replicate by non-ribosomal mechanisms, (iii) rRNA unfailingly increases in thermodynamic stability over time, and (iv) the Woese and Fox canonical tree of life is mis-rooted. Much of the specific GIGO critique of the Accretion Model is based on confusion about the three-dimensional nature of RNA and trunk-branch polymorphism; the Accretion Model incorporates several types of trunk-branch relationships.