Goldilocks and RNA: where Mg2+ concentration is just right.

Magnesium, the most abundant divalent cation in cells, catalyzes RNA cleavage but also promotes RNA folding. Because folding can protect RNA from cleavage, we predicted a 'Goldilocks landscape', with local maximum in RNA lifetime at Mg2+ concentrations required for folding. Here, we use simulation and experiment to discover an innate and sophisticated mechanism of control of RNA lifetime. By simulation we characterized RNA Goldilocks landscapes and their dependence on cleavage and folding parameters. Experiments with yeast tRNAPhe and the Tetrahymena ribozyme P4-P6 domain show that structured RNAs can inhabit Goldilocks peaks. The Goldilocks peaks are tunable by differences in folded and unfolded cleavage rate constants, Mg2+ binding cooperativity, and Mg2+ affinity. Different folding and cleavage parameters produce Goldilocks landscapes with a variety of features. Goldilocks behavior allows ultrafine control of RNA chemical lifetime, whereas non-folding RNAs do not display Goldilocks peaks of protection. In sum, the effects of Mg2+ on RNA persistence are expected to be pleomorphic, both protecting and degrading RNA. In evolutionary context, Goldilocks behavior may have been a selectable trait of RNA in an early Earth environment containing Mg2+ and other metals.

Creative destruction: New protein folds from old.

Mechanisms of emergence and divergence of protein folds pose central questions in biological sciences. Incremental mutation and stepwise adaptation explain relationships between topologically similar protein folds. However, the universe of folds is diverse and riotous, suggesting more potent and creative forces are at play. Sequence and structure similarity are observed between distinct folds, indicating that proteins with distinct folds may share common ancestry. We found evidence of common ancestry between three distinct ß-barrel folds: Scr kinase family homology (SH3), oligonucleotide/oligosaccharide-binding (OB), and cradle loop barrel (CLB). The data suggest a mechanism of fold evolution that interconverts SH3, OB, and CLB. This mechanism, which we call creative destruction, can be generalized to explain many examples of fold evolution including circular permutation. In creative destruction, an open reading frame duplicates or otherwise merges with another to produce a fused polypeptide. A merger forces two ancestral domains into a new sequence and spatial context. The fused polypeptide can explore folding landscapes that are inaccessible to either of the independent ancestral domains. However, the folding landscapes of the fused polypeptide are not fully independent of those of the ancestral domains. Creative destruction is thus partially conservative; a daughter fold inherits some motifs from ancestral folds. After merger and refolding, adaptive processes such as mutation and loss of extraneous segments optimize the new daughter fold. This model has application in disease states characterized by genetic instability. Fused proteins observed in cancer cells are likely to experience remodeled folding landscapes and realize altered folds, conferring new or altered functions.

rRNA expansion segment 7 in eukaryotes: from Signature Fold to tentacles.

The ribosomal core is universally conserved across the tree of life. However, eukaryotic ribosomes contain diverse rRNA expansion segments (ESs) on their surfaces. Sites of ES insertions are predicted from sites of insertion of micro-ESs in archaea. Expansion segment 7 (ES7) is one of the most diverse regions of the ribosome, emanating from a short stem loop and ranging to over 750 nucleotides in mammals. We present secondary and full-atom 3D structures of ES7 from species spanning eukaryotic diversity. Our results are based on experimental 3D structures, the accretion model of ribosomal evolution, phylogenetic relationships, multiple sequence alignments, RNA folding algorithms and 3D modeling by RNAComposer. ES7 contains a distinct motif, the 'ES7 Signature Fold', which is generally invariant in 2D topology and 3D structure in all eukaryotic ribosomes. We establish a model in which ES7 developed over evolution through a series of elementary and recursive growth events. The data are sufficient to support an atomic-level accretion path for rRNA growth. The non-monophyletic distribution of some ES7 features across the phylogeny suggests acquisition via convergent processes. And finally, illustrating the power of our approach, we constructed the 2D and 3D structure of the entire LSU rRNA of Mus musculus.

ProteoVision: web server for advanced visualization of ribosomal proteins.

ProteoVision is a web server designed to explore protein structure and evolution through simultaneous visualization of multiple sequence alignments, topology diagrams and 3D structures. Starting with a multiple sequence alignment, ProteoVision computes conservation scores and a variety of physicochemical properties and simultaneously maps and visualizes alignments and other data on multiple levels of representation. The web server calculates and displays frequencies of amino acids. ProteoVision is optimized for ribosomal proteins but is applicable to analysis of any protein. ProteoVision handles internally generated and user uploaded alignments and connects them with a selected structure, found in the PDB or uploaded by the user. It can generate de novo topology diagrams from three-dimensional structures. All displayed data is interactive and can be saved in various formats as publication quality images or external datasets or PyMol Scripts. ProteoVision enables detailed study of protein fragments defined by Evolutionary Classification of protein Domains (ECOD) classification. ProteoVision is available at http://proteovision.chemistry.gatech.edu/.

Fold Evolution before LUCA: Common Ancestry of SH3 Domains and OB Domains.

SH3 and OB are the simplest, oldest, and most common protein domains within the translation system. SH3 and OB domains are ß-barrels that are structurally similar but are topologically distinct. To transform an OB domain to a SH3 domain, ß-strands must be permuted in a multistep and evolutionarily implausible mechanism. Here, we explored relationships between SH3 and OB domains of ribosomal proteins, initiation, and elongation factors using a combined sequence- and structure-based approach. We detect a common core of SH3 and OB domains, as a region of significant structure and sequence similarity. The common core contains four ß-strands and a loop, but omits the fifth ß-strand, which is variable and is absent from some OB and SH3 domain proteins. The structure of the common core immediately suggests a simple permutation mechanism for interconversion between SH3 and OB domains, which appear to share an ancestor. The OB domain was formed by duplication and adaptation of the SH3 domain core, or vice versa, in a simple and probable transformation. By employing the folding algorithm AlphaFold2, we demonstrated that an ancestral reconstruction of a permuted SH3 sequence folds into an OB structure, and an ancestral reconstruction of a permuted OB sequence folds into a SH3 structure. The tandem SH3 and OB domains in the universal ribosomal protein uL2 share a common ancestor, suggesting that the divergence of these two domains occurred before the last universal common ancestor.

Adaptation and Exaptation: From Small Molecules to Feathers.

Evolution works by adaptation and exaptation. At an organismal level, exaptation and adaptation are seen in the formation of organelles and the advent of multicellularity. At the sub-organismal level, molecular systems such as proteins and RNAs readily undergo adaptation and exaptation. Here we suggest that the concepts of adaptation and exaptation are universal, synergistic, and recursive and apply to small molecules such as metabolites, cofactors, and the building blocks of extant polymers. For example, adenosine has been extensively adapted and exapted throughout biological evolution. Chemical variants of adenosine that are products of adaptation include 2' deoxyadenosine in DNA and a wide array of modified forms in mRNAs, tRNAs, rRNAs, and viral RNAs. Adenosine and its variants have been extensively exapted for various functions, including informational polymers (RNA, DNA), energy storage (ATP), metabolism (e.g., coenzyme A), and signaling (cyclic AMP). According to Gould, Vrba, and Darwin, exaptation imposes a general constraint on interpretation of history and origins; because of exaptation, extant function should not be used to explain evolutionary history. While this notion is accepted in evolutionary biology, it can also guide the study of the chemical origins of life. We propose that (i) evolutionary theory is broadly applicable from the dawn of life to the present time from molecules to organisms, (ii) exaptation and adaptation were important and simultaneous processes, and (iii) robust origin of life models can be constructed without conflating extant utility with historical basis of origins.

TwinCons: Conservation score for uncovering deep sequence similarity and divergence.

We have developed the program TwinCons, to detect noisy signals of deep ancestry of proteins or nucleic acids. As input, the program uses a composite alignment containing pre-defined groups, and mathematically determines a 'cost' of transforming one group to the other at each position of the alignment. The output distinguishes conserved, variable and signature positions. A signature is conserved within groups but differs between groups. The method automatically detects continuous characteristic stretches (segments) within alignments. TwinCons provides a convenient representation of conserved, variable and signature positions as a single score, enabling the structural mapping and visualization of these characteristics. Structure is more conserved than sequence. TwinCons highlights alternative sequences of conserved structures. Using TwinCons, we detected highly similar segments between proteins from the translation and transcription systems. TwinCons detects conserved residues within regions of high functional importance for the ribosomal RNA (rRNA) and demonstrates that signatures are not confined to specific regions but are distributed across the rRNA structure. The ability to evaluate both nucleic acid and protein alignments allows TwinCons to be used in combined sequence and structural analysis of signatures and conservation in rRNA and in ribosomal proteins (rProteins). TwinCons detects a strong sequence conservation signal between bacterial and archaeal rProteins related by circular permutation. This conserved sequence is structurally colocalized with conserved rRNA, indicated by TwinCons scores of rRNA alignments of bacterial and archaeal groups. This combined analysis revealed deep co-evolution of rRNA and rProtein buried within the deepest branching points in the tree of life.

Root of the Tree: The Significance, Evolution, and Origins of the Ribosome.

The ribosome is an ancient molecular fossil that provides a telescope to the origins of life. Made from RNA and protein, the ribosome translates mRNA to coded protein in all living systems. Universality, economy, centrality and antiquity are ingrained in translation. The translation machinery dominates the set of genes that are shared as orthologues across the tree of life. The lineage of the translation system defines the universal tree of life. The function of a ribosome is to build ribosomes; to accomplish this task, ribosomes make ribosomal proteins, polymerases, enzymes, and signaling proteins. Every coded protein ever produced by life on Earth has passed through the exit tunnel, which is the birth canal of biology. During the root phase of the tree of life, before the last common ancestor of life (LUCA), exit tunnel evolution is dominant and unremitting. Protein folding coevolved with evolution of the exit tunnel. The ribosome shows that protein folding initiated with intrinsic disorder, supported through a short, primitive exit tunnel. Folding progressed to thermodynamically stable ß-structures and then to kinetically trapped α-structures. The latter were enabled by a long, mature exit tunnel that partially offset the general thermodynamic tendency of all polypeptides to form ß-sheets. RNA chaperoned the evolution of protein folding from the very beginning. The universal common core of the ribosome, with a mass of nearly 2 million Daltons, was finalized by LUCA. The ribosome entered stasis after LUCA and remained in that state for billions of years. Bacterial ribosomes never left stasis. Archaeal ribosomes have remained near stasis, except for the superphylum Asgard, which has accreted rRNA post LUCA. Eukaryotic ribosomes in some lineages appear to be logarithmically accreting rRNA over the last billion years. Ribosomal expansion in Asgard and Eukarya has been incremental and iterative, without substantial remodeling of pre-existing basal structures. The ribosome preserves information on its history.

Proteotoxic stress promotes entrapment of ribosomes and misfolded proteins in a shared cytosolic compartment.

Cells continuously monitor protein synthesis to prevent accumulation of aberrant polypeptides. Insufficient capacity of cellular degradative systems, chaperone shortage or high levels of mistranslation by ribosomes can result in proteotoxic stress and endanger proteostasis. One of the least explored reasons for mistranslation is the incorrect functioning of the ribosome itself. To understand how cells deal with ribosome malfunction, we introduced mutations in the Expansion Segment 7 (ES7L) of 25S rRNA that allowed the formation of mature, translationally active ribosomes but induced proteotoxic stress and compromised cell viability. The ES7L-mutated ribosomes escaped nonfunctional rRNA Decay (NRD) and remained stable. Remarkably, ES7L-mutated ribosomes showed increased segregation into cytoplasmic foci containing soluble misfolded proteins. This ribosome entrapment pathway, termed TRAP (Translational Relocalization with Aberrant Polypeptides), was generalizable beyond the ES7L mutation, as wild-type ribosomes also showed increased relocalization into the same compartments in cells exposed to proteotoxic stressors. We propose that during TRAP, assembled ribosomes associated with misfolded nascent chains move into cytoplasmic compartments enriched in factors that facilitate protein quality control. In addition, TRAP may help to keep translation at its peak efficiency by preventing malfunctioning ribosomes from active duty in translation.

Cutting in-line with iron: ribosomal function and non-oxidative RNA cleavage.

Divalent metal cations are essential to the structure and function of the ribosome. Previous characterizations of the ribosome performed under standard laboratory conditions have implicated Mg2+ as a primary mediator of ribosomal structure and function. Possible contributions of Fe2+ as a ribosomal cofactor have been largely overlooked, despite the ribosome's early evolution in a high Fe2+ environment, and the continued use of Fe2+ by obligate anaerobes inhabiting high Fe2+ niches. Here, we show that (i) Fe2+ cleaves RNA by in-line cleavage, a non-oxidative mechanism that has not previously been shown experimentally for this metal, (ii) the first-order in-line rate constant with respect to divalent cations is >200 times greater with Fe2+ than with Mg2+, (iii) functional ribosomes are associated with Fe2+ after purification from cells grown under low O2 and high Fe2+ and (iv) a small fraction of Fe2+ that is associated with the ribosome is not exchangeable with surrounding divalent cations, presumably because those ions are tightly coordinated by rRNA and deeply buried in the ribosome. In total, these results expand the ancient role of iron in biochemistry and highlight a possible new mechanism of iron toxicity.

Frontiers in Prebiotic Chemistry and Early Earth Environments.

The Prebiotic Chemistry and Early Earth Environments (PCE3) Consortium is a community of researchers seeking to understand the origins of life on Earth and in the universe. PCE3 is one of five Research Coordination Networks (RCNs) within NASA's Astrobiology Program. Here we report on the inaugural PCE3 workshop, intended to cross-pollinate, transfer information, promote cooperation, break down disciplinary barriers, identify new directions, and foster collaborations. This workshop, entitled, "Building a New Foundation", was designed to propagate current knowledge, identify possibilities for multidisciplinary collaboration, and ultimately define paths for future collaborations. Presentations addressed the likely conditions on early Earth in ways that could be incorporated into prebiotic chemistry experiments and conceptual models to improve their plausibility and accuracy. Additionally, the discussions that followed among workshop participants helped to identify within each subdiscipline particularly impactful new research directions. At its core, the foundational knowledge base presented in this workshop should underpin future workshops and enable collaborations that bridge the many disciplines that are part of PCE3.

Selective incorporation of proteinaceous over nonproteinaceous cationic amino acids in model prebiotic oligomerization reactions.

Numerous long-standing questions in origins-of-life research center on the history of biopolymers. For example, how and why did nature select the polypeptide backbone and proteinaceous side chains? Depsipeptides, containing both ester and amide linkages, have been proposed as ancestors of polypeptides. In this paper, we investigate cationic depsipeptides that form under mild dry-down reactions. We compare the oligomerization of various cationic amino acids, including the cationic proteinaceous amino acids (lysine, Lys; arginine, Arg; and histidine, His), along with nonproteinaceous analogs of Lys harboring fewer methylene groups in their side chains. These analogs, which have been discussed as potential prebiotic alternatives to Lys, are ornithine, 2,4-diaminobutyric acid, and 2,3-diaminopropionic acid (Orn, Dab, and Dpr). We observe that the proteinaceous amino acids condense more extensively than these nonproteinaceous amino acids. Orn and Dab readily cyclize into lactams, while Dab and Dpr condense less efficiently. Furthermore, the proteinaceous amino acids exhibit more selective oligomerization through their α-amines relative to their side-chain groups. This selectivity results in predominantly linear depsipeptides in which the amino acids are α-amine-linked, analogous to today's proteins. These results suggest a chemical basis for the selection of Lys, Arg, and His over other cationic amino acids for incorporation into proto-proteins on the early Earth. Given that electrostatics are key elements of protein-RNA and protein-DNA interactions in extant life, we hypothesize that cationic side chains incorporated into proto-peptides, as reported in this study, served in a variety of functions with ancestral nucleic acid polymers in the early stages of life.

Human ribosomal G-quadruplexes regulate heme bioavailability.

The in vitro formation of stable G-quadruplexes (G4s) in human rRNA was recently reported. However, their formation in cells and their cellular roles were not resolved. Here, by taking a chemical biology approach that integrates results from immunofluorescence, G4 ligands, heme-affinity reagents, and a genetically encoded fluorescent heme sensor, we report that human ribosomes can form G4s in vivo that regulate heme bioavailability. Immunofluorescence experiments indicate that the vast majority of extra-nuclear G4s are associated with rRNA. Moreover, titrating human cells with a G4 ligand alters the ability of ribosomes to bind heme and disrupts cellular heme bioavailability as measured by a genetically encoded fluorescent heme sensor. Overall, these results suggest that ribosomes play a role in regulating heme homeostasis.

A blueprint for academic laboratories to produce SARS-CoV-2 quantitative RT-PCR test kits.

Widespread testing for the presence of the novel coronavirus severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) in individuals remains vital for controlling the COVID-19 pandemic prior to the advent of an effective treatment. Challenges in testing can be traced to an initial shortage of supplies, expertise, and/or instrumentation necessary to detect the virus by quantitative RT-PCR (RT-qPCR), the most robust, sensitive, and specific assay currently available. Here we show that academic biochemistry and molecular biology laboratories equipped with appropriate expertise and infrastructure can replicate commercially available SARS-CoV-2 RT-qPCR test kits and backfill pipeline shortages. The Georgia Tech COVID-19 Test Kit Support Group, composed of faculty, staff, and trainees across the biotechnology quad at Georgia Institute of Technology, synthesized multiplexed primers and probes and formulated a master mix composed of enzymes and proteins produced in-house. Our in-house kit compares favorably with a commercial product used for diagnostic testing. We also developed an environmental testing protocol to readily monitor surfaces for the presence of SARS-CoV-2. Our blueprint should be readily reproducible by research teams at other institutions, and our protocols may be modified and adapted to enable SARS-CoV-2 detection in more resource-limited settings.

Water and Life: The Medium is the Message.

Water, the most abundant compound on the surface of the Earth and probably in the universe, is the medium of biology, but is much more than that. Water is the most frequent actor in the chemistry of metabolism. Our quantitation here reveals that water accounts for 99.4% of metabolites in Escherichia coli by molar concentration. Between a third and a half of known biochemical reactions involve consumption or production of water. We calculated the chemical flux of water and observed that in the life of a cell, a given water molecule frequently and repeatedly serves as a reaction substrate, intermediate, cofactor, and product. Our results show that as an E. coli cell replicates in the presence of molecular oxygen, an average in vivo water molecule is chemically transformed or is mechanistically involved in catalysis ~ 3.7 times. We conclude that, for biological water, there is no distinction between medium and chemical participant. Chemical transformations of water provide a basis for understanding not only extant biochemistry, but the origins of life. Because the chemistry of water dominates metabolism and also drives biological synthesis and degradation, it seems likely that metabolism co-evolved with biopolymers, which helps to reconcile polymer-first versus metabolism-first theories for the origins of life.

Multiple prebiotic metals mediate translation.

Today, Mg2+ is an essential cofactor with diverse structural and functional roles in life's oldest macromolecular machine, the translation system. We tested whether ancient Earth conditions (low O2, high Fe2+, and high Mn2+) can revert the ribosome to a functional ancestral state. First, SHAPE (selective 2'-hydroxyl acylation analyzed by primer extension) was used to compare the effect of Mg2+, Fe2+, and Mn2+ on the tertiary structure of rRNA. Then, we used in vitro translation reactions to test whether Fe2+ or Mn2+ could mediate protein production, and quantified ribosomal metal content. We found that (i) Mg2+, Fe2+, and Mn2+ had strikingly similar effects on rRNA folding; (ii) Fe2+ and Mn2+ can replace Mg2+ as the dominant divalent cation during translation of mRNA to functional protein; and (iii) Fe and Mn associate extensively with the ribosome. Given that the translation system originated and matured when Fe2+ and Mn2+ were abundant, these findings suggest that Fe2+ and Mn2+ played a role in early ribosomal evolution.

Translation: The Universal Structural Core of Life.

The Universal Gene Set of Life (UGSL) is common to genomes of all extant organisms. The UGSL is small, consisting of <100 genes, and is dominated by genes encoding the translation system. Here we extend the search for biological universality to three dimensions. We characterize and quantitate the universality of structure of macromolecules that are common to all of life. We determine that around 90% of prokaryotic ribosomal RNA (rRNA) forms a common core, which is the structural and functional foundation of rRNAs of all cytoplasmic ribosomes. We have established a database, which we call the Sparse and Efficient Representation of the Extant Biology (the SEREB database). This database contains complete and cross-validated rRNA sequences of species chosen, as far as possible, to sparsely and efficiently sample all known phyla. Atomic-resolution structures of ribosomes provide data for structural comparison and validation of sequence-based models. We developed a similarity statistic called pairing adjusted sequence entropy, which characterizes paired nucleotides by their adherence to covariation and unpaired nucleotides by conventional conservation of identity. For canonically paired nucleotides the unit of structure is the nucleotide pair. For unpaired nucleotides, the unit of structure is the nucleotide. By quantitatively defining the common core of rRNA, we systematize the conservation and divergence of the translational system across the tree of life, and can begin to understand the unique evolutionary pressures that cause its universality. We explore the relationship between ribosomal size and diversity, geological time, and organismal complexity.

Iron mediates catalysis of nucleic acid processing enzymes: support for Fe(II) as a cofactor before the great oxidation event.

Life originated in an anoxic, Fe2+-rich environment. We hypothesize that on early Earth, Fe2+ was a ubiquitous cofactor for nucleic acids, with roles in RNA folding and catalysis as well as in processing of nucleic acids by protein enzymes. In this model, Mg2+ replaced Fe2+ as the primary cofactor for nucleic acids in parallel with known metal substitutions of metalloproteins, driven by the Great Oxidation Event. To test predictions of this model, we assay the ability of nucleic acid processing enzymes, including a DNA polymerase, an RNA polymerase and a DNA ligase, to use Fe2+ in place of Mg2+ as a cofactor during catalysis. Results show that Fe2+ can indeed substitute for Mg2+ in catalytic function of these enzymes. Additionally, we use calculations to unravel differences in energetics, structures and reactivities of relevant Mg2+ and Fe2+ complexes. Computation explains why Fe2+ can be a more potent cofactor than Mg2+ in a variety of folding and catalytic functions. We propose that the rise of O2 on Earth drove a Fe2+ to Mg2+ substitution in proteins and nucleic acids, a hypothesis consistent with a general model in which some modern biochemical systems retain latent abilities to revert to primordial Fe2+-based states when exposed to pre-GOE conditions.

Functional RNAs: combined assembly and packaging in VLPs.

We describe here a one pot RNA production, packaging and delivery system based on bacteriophage Qß. We demonstrate a method for production of a novel RNAi scaffold, packaged within Qß virus-like particles (VLPs). The RNAi scaffold is a general utility chimera that contains a functional RNA duplex with paired silencing and carrier sequences stabilized by a miR-30 stem-loop. The Qß hairpin on the 5΄ end confers affinity for the Qß coat protein (CP). Silencing sequences can include mature miRNAs and siRNAs, and can target essentially any desired mRNA. The VLP-RNAi assembles upon co-expression of CP and the RNAi scaffold in E. coli. The annealing of the scaffold to form functional RNAs is intramolecular and is therefore robust and concentration independent. We demonstrate dose- and time-dependent inhibition of GFP expression in human cells with VLP-RNAi. In addition, we target the 3΄UTR of oncogenic Ras mRNA and suppress Pan-Ras expression, which attenuates cell proliferation and promotes mortality of brain tumor cells. This combination of RNAi scaffold design with Qß VLP packaging is demonstrated to be target-specific and efficient.

Frozen in Time: The History of Proteins.

The ribosome is imprinted with a detailed molecular chronology of the origins and early evolution of proteins. Here we show that when arranged by evolutionary phase of ribosomal evolution, ribosomal protein (rProtein) segments reveal an atomic level history of protein folding. The data support a model in which aboriginal oligomers evolved into globular proteins in a hierarchical step-wise process. Complexity of assembly and folding of polypeptide increased incrementally in concert with expansion of rRNA. (i) Short random coil proto-peptides bound to rRNA, and (ii) lengthened over time and coalesced into ß-ß secondary elements. These secondary elements (iii) accreted and collapsed, primarily into ß-domains. Domains (iv) accumulated and gained complex super-secondary structures composed of mixtures of α-helices and ß-strands. Early protein evolution was guided and accelerated by interactions with rRNA. rRNA and proto-peptide provided mutual protection from chemical degradation and disassembly. rRNA stabilized polypeptide assemblies, which evolved in a stepwise process into globular domains, bypassing the immense space of random unproductive sequences. Coded proteins originated as oligomers and polymers created by the ribosome, on the ribosome and for the ribosome. Synthesis of increasingly longer products was iteratively coupled with lengthening and maturation of the ribosomal exit tunnel. Protein catalysis appears to be a late byproduct of selection for sophisticated and finely controlled assembly.