Peak Cenozoic warmth enabled deep-sea sand deposition.

The early Eocene (~ 56-48 million years ago) was marked by peak Cenozoic warmth and sea levels, high CO2, and largely ice-free conditions. This time has been described as a period of increased continental erosion and silicate weathering. However, these conclusions are based largely on geochemical investigation of marine mudstones and carbonates or study of intermontane Laramide basin settings. Here, we evaluate the marine coarse siliciclastic response to early Paleogene hothouse climatic and oceanographic conditions. We compile an inventory of documented sand-rich (turbidite) deep-marine depositional systems, recording 59 instances of early Eocene turbidite systems along nearly all continental margins despite globally-elevated sea levels. Sand-rich systems were widespread on active margins (42 instances), but also on passive margins (17 instances). Along passive margins, 13 of 17 early Eocene systems are associated with known Eocene-age fluvial systems, consistent with a fluvial clastic response to Paleogene warming. We suggest that deep-marine sedimentary basins preserve clastic records of early Eocene climatic extremes. We also suggest that in addition to control by eustasy and tectonism, climate-driven increases in sediment supply (e.g., drainage integration, global rainfall, denudation) may significantly contribute to the global distribution and volume of coarse-grained deep-marine deposition despite high sea level.

Navigating mental health challenges in graduate school.

Many graduate students experience mental health struggles that lead them to question their place in academia. Two scientists who experienced extreme lows in graduate school reflect on what helped them during their low points, and suggest strategies for everyone to contribute to mentally healthier workplaces in academia.

Early Evolution of Transcription Systems and Divergence of Archaea and Bacteria.

DNA template-dependent multi-subunit RNA polymerases (RNAPs) found in all three domains of life and some viruses are of the two-double-Ψ-ß-barrel (DPBB) type. The 2-DPBB protein format is also found in some RNA template-dependent RNAPs and a major replicative DNA template-dependent DNA polymerase (DNAP) from Archaea (PolD). The 2-DPBB family of RNAPs and DNAPs probably evolved prior to the last universal common cellular ancestor (LUCA). Archaeal Transcription Factor B (TFB) and bacterial σ factors include homologous strings of helix-turn-helix units. The consequences of TFB-σ homology are discussed in terms of the evolution of archaeal and bacterial core promoters. Domain-specific DPBB loop inserts functionally connect general transcription factors to the RNAP active site. Archaea appear to be more similar to LUCA than Bacteria. Evolution of bacterial σ factors from TFB appears to have driven divergence of Bacteria from Archaea, splitting the prokaryotic domains.

Evolution of Life on Earth: tRNA, Aminoacyl-tRNA Synthetases and the Genetic Code.

Life on Earth and the genetic code evolved around tRNA and the tRNA anticodon. We posit that the genetic code initially evolved to synthesize polyglycine as a cross-linking agent to stabilize protocells. We posit that the initial amino acids to enter the code occupied larger sectors of the code that were then invaded by incoming amino acids. Displacements of amino acids follow selection rules. The code sectored from a glycine code to a four amino acid code to an eight amino acid code to an ~16 amino acid code to the standard 20 amino acid code with stops. The proposed patterns of code sectoring are now most apparent from patterns of aminoacyl-tRNA synthetase evolution. The Elongation Factor-Tu GTPase anticodon-codon latch that checks the accuracy of translation appears to have evolved at about the eight amino acid to ~16 amino acid stage. Before evolution of the EF-Tu latch, we posit that both the 1st and 3rd anticodon positions were wobble positions. The genetic code evolved via tRNA charging errors and via enzymatic modifications of amino acids joined to tRNAs, followed by tRNA and aminoacyl-tRNA synthetase differentiation. Fidelity mechanisms froze the code by inhibiting further innovation.

The 3-Minihelix tRNA Evolution Theorem.

Transfer RNA (tRNA) is the central intellectual property in the evolution of life on Earth. tRNA evolved from repeats and inverted repeats of known sequence. The anticodon and the T stem-loop-stems are homologs with significant conserved sequence identity. A number of models have been advanced to explain tRNA evolution. No 2-minihelix model or accretion model (built a stem at a time) can be correct, in part because of anticodon and T stem-loop-stem identity. Only a 3-minihelix model is adequate.

A tRNA- and Anticodon-Centric View of the Evolution of Aminoacyl-tRNA Synthetases, tRNAomes, and the Genetic Code.

Pathways of standard genetic code evolution remain conserved and apparent, particularly upon analysis of aminoacyl-tRNA synthetase (aaRS) lineages. Despite having incompatible active site folds, class I and class II aaRS are homologs by sequence. Specifically, structural class IA aaRS enzymes derive from class IIA aaRS enzymes by in-frame extension of the protein N-terminus and by an alternate fold nucleated by the N-terminal extension. The divergence of aaRS enzymes in the class I and class II clades was analyzed using the Phyre2 protein fold recognition server. The class I aaRS radiated from the class IA enzymes, and the class II aaRS radiated from the class IIA enzymes. The radiations of aaRS enzymes bolster the coevolution theory for evolution of the amino acids, tRNAomes, the genetic code, and aaRS enzymes and support a tRNA anticodon-centric perspective. We posit that second- and third-position tRNA anticodon sequence preference (C>(U~G)>A) powerfully selected the sectoring pathway for the code. GlyRS-IIA appears to have been the primordial aaRS from which all aaRS enzymes evolved, and glycine appears to have been the primordial amino acid around which the genetic code evolved.

Ribosome Structure, Function, and Early Evolution.

Ribosomes are among the largest and most dynamic molecular motors. The structure and dynamics of translation initiation and elongation are reviewed. Three ribosome motions have been identified for initiation and translocation. A swivel motion between the head/beak and the body of the 30S subunit was observed. A tilting dynamic of the head/beak versus the body of the 30S subunit was detected using simulations. A reversible ratcheting motion was seen between the 30S and the 50S subunits that slide relative to one another. The 30S⁻50S intersubunit contacts regulate translocation. IF2, EF-Tu, and EF-G are homologous G-protein GTPases that cycle on and off the same site on the ribosome. The ribosome, aminoacyl-tRNA synthetase (aaRS) enzymes, transfer ribonucleic acid (tRNA), and messenger ribonucleic acid (mRNA) form the core of information processing in cells and are coevolved. Surprisingly, class I and class II aaRS enzymes, with distinct and incompatible folds, are homologs. Divergence of class I and class II aaRS enzymes and coevolution of the genetic code are described by analysis of ancient archaeal species.

Type-II tRNAs and Evolution of Translation Systems and the Genetic Code.

Because tRNA is the core biological intellectual property that was necessary to evolve translation systems, tRNAomes, ribosomes, aminoacyl-tRNA synthetases, and the genetic code, the evolution of tRNA is the core story in evolution of life on earth. We have previously described the evolution of type-I tRNAs. Here, we use the same model to describe the evolution of type-II tRNAs, with expanded V loops. The models are strongly supported by inspection of typical tRNA diagrams, measuring lengths of V loop expansions, and analyzing the homology of V loop sequences to tRNA acceptor stems. Models for tRNA evolution provide a pathway for the inanimate-to-animate transition and for the evolution of translation systems, the genetic code, and cellular life.

Aminoacyl-tRNA synthetase evolution and sectoring of the genetic code.

The genetic code sectored via tRNA charging errors, and the code progressed toward closure and universality because of evolution of aminoacyl-tRNA synthetase (aaRS) fidelity and translational fidelity mechanisms. Class I and class II aaRS folds are identified as homologs. From sequence alignments, a structurally conserved Zn-binding domain common to class I and class II aaRS was identified. A model for the class I and class II aaRS alternate folding pathways is posited. Five mechanisms toward code closure are highlighted: 1) aaRS proofreading to remove mischarged amino acids from tRNA; 2) accurate aaRS active site specification of amino acid substrates; 3) aaRS-tRNA anticodon recognition; 4) conformational coupling proofreading of the anticodon-codon interaction; and 5) deamination of tRNA wobble adenine to inosine. In tRNA anticodons there is strong wobble sequence preference that results in a broader spectrum of contacts to synonymous mRNA codon wobble bases. Adenine is excluded from the anticodon wobble position of tRNA unless it is modified to inosine. Uracil is generally preferred to cytosine in the tRNA anticodon wobble position. Because of wobble ambiguity when tRNA reads mRNA, the maximal coding capacity of the three nucleotide code read by tRNA is 31 amino acids + stops.

Homology threading to generate RNA polymerase structures.

Homology threading is a powerful technology for generating structural models based on homologous structures. Here we use threading to generate four complex RNA polymerase models. The models appear to be as useful as x-ray crystal structures or cryo-electron microscopy structures to support research projects.

Rooted tRNAomes and evolution of the genetic code.

We advocate for a tRNA- rather than an mRNA-centric model for evolution of the genetic code. The mechanism for evolution of cloverleaf tRNA provides a root sequence for radiation of tRNAs and suggests a simplified understanding of code evolution. To analyze code sectoring, rooted tRNAomes were compared for several archaeal and one bacterial species. Rooting of tRNAome trees reveals conserved structures, indicating how the code was shaped during evolution and suggesting a model for evolution of a LUCA tRNAome tree. We propose the polyglycine hypothesis that the initial product of the genetic code may have been short chain polyglycine to stabilize protocells. In order to describe how anticodons were allotted in evolution, the sectoring-degeneracy hypothesis is proposed. Based on sectoring, a simple stepwise model is developed, in which the code sectors from a 1→4→8→∼16 letter code. At initial stages of code evolution, we posit strong positive selection for wobble base ambiguity, supporting convergence to 4-codon sectors and ∼16 letters. In a later stage, ∼5-6 letters, including stops, were added through innovating at the anticodon wobble position. In archaea and bacteria, tRNA wobble adenine is negatively selected, shrinking the maximum size of the primordial genetic code to 48 anticodons. Because 64 codons are recognized in mRNA, tRNA-mRNA coevolution requires tRNA wobble position ambiguity leading to degeneracy of the code.

Hinge action versus grip in translocation by RNA polymerase.

Based on molecular dynamics simulations and functional studies, a conformational mechanism is posited for forward translocation by RNA polymerase (RNAP). In a simulation of a ternary elongation complex, the clamp and downstream cleft were observed to close. Hinges within the bridge helix and trigger loop supported generation of translocation force against the RNA-DNA hybrid resulting in opening of the furthest upstream i-8 RNA-DNA bp, establishing conditions for RNAP sliding. The ß flap tip helix and the most N-terminal ß' Zn finger engage the RNA, indicating a path of RNA threading out of the exit channel. Because the ß flap tip connects to the RNAP active site through the ß subunit double-Ψ-ß-barrel and the associated sandwich barrel hybrid motif (also called the flap domain), the RNAP active site is coupled to the RNA exit channel and to the translocation of RNA-DNA. Using an exonuclease III assay to monitor translocation of RNAP elongation complexes, we show that K+ and Mg2+ and also an RNA 3'-OH or a 3'-H2 affect RNAP sliding. Because RNAP grip to template suggests a sticky translocation mechanism, and because grip is enhanced by increasing K+ and Mg2+concentration, biochemical assays are consistent with a conformational change that drives forward translocation as observed in simulations. Mutational analysis of the bridge helix indicates that 778-GARKGL-783 (Escherichia coli numbering) is a homeostatic hinge that undergoes multiple bends to compensate for complex conformational dynamics during phosphodiester bond formation and translocation.

tRNA structure and evolution and standardization to the three nucleotide genetic code.

Cloverleaf tRNA with a 75 nucleotide (nt) core is posited to have evolved from ligation of three 31 nt minihelices followed by symmetric internal deletions of 9 nt within ligated acceptor stems. Statistical tests strongly support the model. Although the tRNA anticodon loop and T loop are homologs, their U-turns have been treated as distinct motifs. An appropriate comparison, however, shows that intercalation of D loop G19 between T loop bases 4 and 5 causes elevation of T loop base 5 and flipping of T loop bases 6 and 7 out of the 7 nt loop. In the anticodon loop, by contrast, loop bases 3-7 stack tightly to form a stiff connection to mRNA. Furthermore, we identify ancient repeat sequences of 3 (GCG), 5 (UAGCC) and 17 nt (∼CCGGGUUCAAAACCCGG) that comprise 75 out of 75 nts of the tRNA cloverleaf core. To present a sufficiently stiff 3-nt anticodon, a 7-nt anticodon loop was necessary with a U-turn between loop positions 2 and 3. Cloverleaf tRNA, therefore, was a radical evolutionary innovation essential for the 3-nt code. Conservation of GCG and UAGCC repeat sequences indicates that cloverleaf tRNA is at the interface between a strange RNA repeat world and the first evolution of molecules that fold to assume biologic functions. We posit that cloverleaf tRNA was the molecular archetype around which translation systems evolved.

tRNA evolution from the proto-tRNA minihelix world.

Multiple models have been advanced for the evolution of cloverleaf tRNA. Here, the conserved archaeal tRNA core (75-nt) is posited to have evolved from ligation of three proto-tRNA minihelices (31-nt) and two-symmetrical 9-nt deletions within joined acceptor stems (93 - 18 = 75-nt). The primary evidence for this conclusion is that the 5-nt stem 7-nt anticodon loop and the 5-nt stem 7-nt T loop are structurally homologous and related by coding sequence. We posit that the D loop was generated from a third minihelix (31-nt) in which the stem and loop became rearranged after 9-nt acceptor stem deletions and cloverleaf folding. The most 3´-5-nt segment of the D loop and the 5-nt V loop are apparent remnants of the joined acceptor stems (14 - 9 = 5-nt). Before refolding in the tRNA cloverleaf, we posit that the 3'-5-nt segment of the D loop and the 5-nt V loop were paired, and, in the tRNA cloverleaf, frequent pairing of positions 29 (D loop) and 47 (V loop) remains (numbered on a 75-nt tRNA cloverleaf core). Amazingly, after >3.5 billion years of evolutionary pressure on the tRNA cloverleaf structure, a model can be constructed that convincingly describes the genesis of 75/75-nt conserved archaeal tRNA core positions. Judging from the tRNA structure, cloverleaf tRNA appears to represent at least a second-generation scheme (and possibly a third-generation scheme) that replaced a robust 31-nt minihelix protein-coding system, evidence for which is preserved in the cloverleaf structure. Understanding tRNA evolution provides insights into ribosome and rRNA evolution.

A model for genesis of transcription systems.

Repeating sequences generated from RNA gene fusions/ligations dominate ancient life, indicating central importance of building structural complexity in evolving biological systems. A simple and coherent story of life on earth is told from tracking repeating motifs that generate α/ß proteins, 2-double-Ψ-ß-barrel (DPBB) type RNA polymerases (RNAPs), general transcription factors (GTFs), and promoters. A general rule that emerges is that biological complexity that arises through generation of repeats is often bounded by solubility and closure (i.e., to form a pseudo-dimer or a barrel). Because the first DNA genomes were replicated by DNA template-dependent RNA synthesis followed by RNA template-dependent DNA synthesis via reverse transcriptase, the first DNA replication origins were initially 2-DPBB type RNAP promoters. A simplifying model for evolution of promoters/replication origins via repetition of core promoter elements is proposed. The model can explain why Pribnow boxes in bacterial transcription (i.e., (-12)TATAATG(-6)) so closely resemble TATA boxes (i.e., (-31)TATAAAAG(-24)) in archaeal/eukaryotic transcription. The evolution of anchor DNA sequences in bacterial (i.e., (-35)TTGACA(-30)) and archaeal (BRE(up); BRE for TFB recognition element) promoters is potentially explained. The evolution of BRE(down) elements of archaeal promoters is potentially explained.

Five checkpoints maintaining the fidelity of transcription by RNA polymerases in structural and energetic details.

Transcriptional fidelity, which prevents the misincorporation of incorrect nucleoside monophosphates in RNA, is essential for life. Results from molecular dynamics (MD) simulations of eukaryotic RNA polymerase (RNAP) II and bacterial RNAP with experimental data suggest that fidelity may involve as many as five checkpoints. Using MD simulations, the effects of different active site NTPs in both open and closed trigger loop (TL) structures of RNAPs are compared. Unfavorable initial binding of mismatched substrates in the active site with an open TL is proposed to be the first fidelity checkpoint. The leaving of an incorrect substrate is much easier than a correct one energetically from the umbrella sampling simulations. Then, the closing motion of the TL, required for catalysis, is hindered by the presence of mismatched NTPs. Mismatched NTPs also lead to conformational changes in the active site, which perturb the coordination of magnesium ions and likely affect the ability to proceed with catalysis. This step appears to be the most important checkpoint for deoxy-NTP discrimination. Finally, structural perturbations in the template DNA and the nascent RNA in the presence of mismatches likely hinder nucleotide addition and provide the structural foundation for backtracking followed by removing erroneously incorporated nucleotides during proofreading.

The σ enigma: bacterial σ factors, archaeal TFB and eukaryotic TFIIB are homologs.

Structural comparisons of initiating RNA polymerase complexes and structure-based amino acid sequence alignments of general transcription initiation factors (eukaryotic TFIIB, archaeal TFB and bacterial σ factors) show that these proteins are homologs. TFIIB and TFB each have two-five-helix cyclin-like repeats (CLRs) that include a C-terminal helix-turn-helix (HTH) motif (CLR/HTH domains). Four homologous HTH motifs are present in bacterial σ factors that are relics of CLR/HTH domains. Sequence similarities clarify models for σ factor and TFB/TFIIB evolution and function and suggest models for promoter evolution. Commitment to alternate modes for transcription initiation appears to be a major driver of the divergence of bacteria and archaea.

The Old and New Testaments of gene regulation. Evolution of multi-subunit RNA polymerases and co-evolution of eukaryote complexity with the RNAP II CTD.

I relate a story of genesis told from the point of view of multi-subunit RNA polymerases (RNAPs) including an Old Testament (core RNAP motifs in all cellular life) and a New Testament (the RNAP II heptad repeat carboxy terminal domain (CTD) and CTD interactome in eukarya). The Old Testament: at their active site, one class of eukaryotic interfering RNAP and ubiquitous multi-subunit RNAPs each have two-double psi ß barrel (DPBB) motifs (a distinct pattern for compact 6-ß sheet barrels). Between ß sheets 2 and 3 of the ß subunit type DPBB of all multi-subunit RNAPs is a sandwich barrel hybrid motif (SBHM) that interacts with conserved initiation and elongation factors required to utilize a DNA template. Analysis of RNAP core protein motifs, therefore, indicates that RNAP evolution can be traced from the RNA-protein world to LUCA (the last universal common ancestor) branching to LECA (the last eukaryotic common ancestor) and to the present day, spanning about 4 billion years. The New Testament: in the eukaryotic lineage, I posit that splitting RNAP functions into RNAPs I, II and III and innovations developed around the CTD heptad repeat of RNAP II and the extensive CTD interactome helps to describe how greater structural, cell cycle, epigenetic and signaling complexity co-evolved in eukaryotes relative to eubacteria and archaea.

Energetic and structural details of the trigger-loop closing transition in RNA polymerase II.

An evolutionarily conserved element in RNA polymerase II, the trigger loop (TL), has been suggested to play an important role in the elongation rate, fidelity of selection of the matched nucleoside triphosphate (NTP), catalysis of transcription elongation, and translocation in both eukaryotes and prokaryotes. In response to NTP binding, the TL undergoes large conformational changes to switch between distinct open and closed states to tighten the active site and avail catalysis. A computational strategy for characterizing the conformational transition pathway is presented to bridge the open and closed states of the TL. Information from a large number of independent all-atom molecular dynamics trajectories from Hamiltonian replica exchange and targeted molecular dynamics simulations is gathered together to assemble a connectivity map of the conformational transition. The results show that with a cognate NTP, TL closing should be a spontaneous process. One major intermediate state is identified along the conformational transition pathway, and the key structural features are characterized. The complete pathway from the open TL to the closed TL provides a clear picture of the TL closing.