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

The Possible Crystallization Process in the Origin of Bacteria, Archaea, Viruses, and Mobile Elements

We propose a hypothesis for the simultaneous emergence of bacteria, archaea, viruses, and mobile elements by sequential and concrete biochemical pathways. The emergence process can be considered analogous to crystallization, where genetic and biochemical systems stabilize as organisms evolve from their common ancestor, the LUCA, which was a non-free-living pool of single operon type genomes including double-stranded (ds) DNA at an ancient submarine alkaline vent. Each dsDNA operon was transcribed by different systems in σ, TFIIB, or TBP genomes. Double-stranded DNA operons can fuse and stabilize through the action of specific transcription systems, leading to differentiation between the Bacteria (σ genome) and Archaea (TBP genome) domains. Error catastrophe can be overcome by the parallel gain of DNA replication and DNA repair mechanisms in both genomes. Enlarged DNA enabled efficient local biochemical reactions. Both genomes independently recruited lipids to facilitate reactions by forming coacervates at the chamber of the vent. Bilayer lipid membrane formation, proto-cell formation with a permeable membrane, proto-cell division, and the evolution of membrane-associated biochemistry are presented in detail. Simultaneous crystallization of systems in non-free-living bacteria and non-free-living archaea triggered the co-crystallization of primitive viruses and mobile elements. An arms race between non-free-living cells and primitive viruses finally led to free-living cells with a cell wall and mature viruses.

Transcription elongation mechanisms of RNA polymerases I, II, and III and their therapeutic implications

Transcription is a tightly regulated, complex, and essential cellular process in all living organisms. Transcription is comprised of three steps, transcription initiation, elongation, and termination. The distinct transcription initiation and termination mechanisms of eukaryotic RNA polymerases I, II, and III (Pols I, II, and III) have long been appreciated. Recent methodological advances have empowered high-resolution investigations of the Pols' transcription elongation mechanisms. Here, we review the kinetic similarities and differences in the individual steps of Pol I-, II-, and III-catalyzed transcription elongation, including NTP binding, bond formation, pyrophosphate release, and translocation. This review serves as an important summation of Saccharomyces cerevisiae (yeast) Pol I, II, and III kinetic investigations which reveal that transcription elongation by the Pols is governed by distinct mechanisms. Further, these studies illustrate how basic, biochemical investigations of the Pols can empower the development of chemotherapeutic compounds.

Distinct states of nucleolar stress induced by anticancer drugs

Ribosome biogenesis is a vital and highly energy-consuming cellular function occurring primarily in the nucleolus. Cancer cells have an elevated demand for ribosomes to sustain continuous proliferation. This study evaluated the impact of existing anticancer drugs on the nucleolus by screening a library of anticancer compounds for drugs that induce nucleolar stress. For a readout, a novel parameter termed 'nucleolar normality score' was developed that measures the ratio of the fibrillar center and granular component proteins in the nucleolus and nucleoplasm. Multiple classes of drugs were found to induce nucleolar stress, including DNA intercalators, inhibitors of mTOR/PI3K, heat shock proteins, proteasome, and cyclin-dependent kinases (CDKs). Each class of drugs induced morphologically and molecularly distinct states of nucleolar stress accompanied by changes in nucleolar biophysical properties. In-depth characterization focused on the nucleolar stress induced by inhibition of transcriptional CDKs, particularly CDK9, the main CDK that regulates RNA Pol II. Multiple CDK substrates were identified in the nucleolus, including RNA Pol I- recruiting protein Treacle, which was phosphorylated by CDK9 in vitro. These results revealed a concerted regulation of RNA Pol I and Pol II by transcriptional CDKs. Our findings exposed many classes of chemotherapy compounds that are capable of inducing nucleolar stress, and we recommend considering this in anticancer drug development.

Uncovering the mechanisms of transcription elongation by eukaryotic RNA polymerases I, II, and III

Eukaryotes express three nuclear RNA polymerases (Pols I, II, and III) that are essential for cell survival. Despite extensive investigation of the three Pols, significant knowledge gaps regarding their biochemical properties remain because each Pol has been evaluated independently under disparate experimental conditions and methodologies. To advance our understanding of the Pols, we employed identical in vitro transcription assays for direct comparison of their elongation rates, elongation complex (EC) stabilities, and fidelities. Pol I is the fastest, most likely to misincorporate, forms the least stable EC, and is most sensitive to alterations in reaction buffers. Pol II is the slowest of the Pols, forms the most stable EC, and negligibly misincorporated an incorrect nucleotide. The enzymatic properties of Pol III were intermediate between Pols I and II in all assays examined. These results reveal unique enzymatic characteristics of the Pols that provide new insights into their evolutionary divergence.

Cell Cycle-Dependent Transcription: The Cyclin Dependent Kinase Cdk1 Is a Direct Regulator of Basal Transcription Machineries

The cyclin-dependent kinase Cdk1 is best known for its function as master regulator of the cell cycle. It phosphorylates several key proteins to control progression through the different phases of the cell cycle. However, studies conducted several decades ago with mammalian cells revealed that Cdk1 also directly regulates the basal transcription machinery, most notably RNA polymerase II. More recent studies in the budding yeast Saccharomyces cerevisiae have revisited this function of Cdk1 and also revealed that Cdk1 directly controls RNA polymerase III activity. These studies have also provided novel insight into the physiological relevance of this process. For instance, cell cycle-stage-dependent activity of these complexes may be important for meeting the increased demand for various proteins involved in housekeeping, metabolism, and protein synthesis. Recent work also indicates that direct regulation of the RNA polymerase II machinery promotes cell cycle entry. Here, we provide an overview of the regulation of basal transcription by Cdk1, and we hypothesize that the original function of the primordial cell-cycle CDK was to regulate RNAPII and that it later evolved into specialized kinases that govern various aspects of the transcription machinery and the cell cycle.

Seven Amino Acid Types Suffice to Create the Core Fold of RNA Polymerase

The extant complex proteins must have evolved from ancient short and simple ancestors. The double-ψ β-barrel (DPBB) is one of the oldest protein folds and conserved in various fundamental enzymes, such as the core domain of RNA polymerase. Here, by reverse engineering a modern DPBB domain, we reconstructed its plausible evolutionary pathway started by "interlacing homodimerization" of a half-size peptide, followed by gene duplication and fusion. Furthermore, by simplifying the amino acid repertoire of the peptide, we successfully created the DPBB fold with only seven amino acid types (Ala, Asp, Glu, Gly, Lys, Arg, and Val), which can be coded by only GNN and ARR (R = A or G) codons in the modern translation system. Thus, the DPBB fold could have been materialized by the early translation system and genetic code.

Removing quote marks from the RNA polymerase II CTD 'code'

In eukaryotes, RNA polymerase II (Pol II) is responsible for the synthesis of all mRNAs and myriads of short and long untranslated RNAs, whose fabrication involves close spatiotemporal coordination between transcription, RNA processing and chromatin modification. Crucial for such a coordination is an unusual C-terminal domain (CTD) of the Pol II largest subunit, made of tandem repetitions (26 in yeast, 52 in chordates) of the heptapeptide with the consensus sequence YSPTSPS. Although largely unstructured and with poor sequence content, the Pol II CTD derives its extraordinary functional versatility from the fact that each amino acid in the heptapeptide can be posttranslationally modified, and that different combinations of CTD covalent marks are specifically recognized by different protein binding partners. These features have led to propose the existence of a Pol II CTD code, but this expression is generally used by authors with some caution, revealed by the frequent use of quote marks for the word 'code'. Based on the theoretical framework of code biology, it is argued here that the Pol II CTD modification system meets the requirements of a true organic code, where different CTD modification states represent organic signs whose organic meanings are biological reactions contributing to the many facets of RNA biogenesis in coordination with RNA synthesis by Pol II. Importantly, the Pol II CTD code is instantiated by adaptor proteins possessing at least two distinct domains, one of which devoted to specific recognition of CTD modification profiles. Furthermore, code rules can be altered by experimental interchange of CTD recognition domains of different adaptor proteins, a fact arguing in favor of the arbitrariness, and thus bona fide character, of the Pol II CTD code. Since the growing family of CTD adaptors includes RNA binding proteins and histone modification complexes, the Pol II CTD code is by its nature integrated with other organic codes, in particular the splicing code and the histone code. These issues will be discussed taking into account fascinating developments in Pol II CTD research, like the discovery of novel modifications at non-consensus sites, the recently recognized CTD physicochemical properties favoring liquid-liquid phase separation, and the discovery that the Pol II CTD, originated before the divergence of most extant eukaryotic taxa, has expanded and diversified with developmental complexity in animals and plants.

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.

Eukaryotic RNA Polymerases: The Many Ways to Transcribe a Gene

In eukaryotic cells, three nuclear RNA polymerases (RNA pols) carry out the transcription from DNA to RNA, and they all seem to have evolved from a single enzyme present in the common ancestor with archaea. The multiplicity of eukaryotic RNA pols allows each one to remain specialized in the synthesis of a subset of transcripts, which are different in the function, length, cell abundance, diversity, and promoter organization of the corresponding genes. We hypothesize that this specialization of RNA pols has conditioned the evolution of the regulatory mechanisms used to transcribe each gene subset to cope with environmental changes. We herein present the example of the homeostatic regulation of transcript levels versus changes in cell volume. We propose that the diversity and instability of messenger RNAs, transcribed by RNA polymerase II, have conditioned the appearance of regulatory mechanisms based on different gene promoter strength and mRNA stability. However, for the regulation of ribosomal RNA levels, which are very stable and transcribed mainly by RNA polymerase I from only one promoter, different mechanisms act based on gene copy variation, and a much simpler regulation of the synthesis rate.

Evolution of the genetic code

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

The negative elongation factor NELF promotes induced transcriptional response of Drosophila ecdysone-dependent genes

For many years it was believed that promoter-proximal RNA-polymerase II (Pol II) pausing manages the transcription of genes in Drosophila development by controlling spatiotemporal properties of their activation and repression. But the exact proteins that cooperate to stall Pol II in promoter-proximal regions of developmental genes are still largely unknown. The current work describes the molecular mechanism employed by the Negative ELongation Factor (NELF) to control the Pol II pause at genes whose transcription is induced by 20-hydroxyecdysone (20E). According to our data, the NELF complex is recruited to the promoters and enhancers of 20E-dependent genes. Its presence at the regulatory sites of 20E-dependent genes correlates with observed interaction between the NELF-A subunit and the ecdysone receptor (EcR). The complete NELF complex is formed at the 20E-dependent promoters and participates in both their induced transcriptional response and maintenance of the uninduced state to keep them ready for the forthcoming transcription. NELF depletion causes a significant decrease in transcription induced by 20E, which is associated with the disruption of Pol II elongation complexes. A considerable reduction in the promoter-bound level of the Spt5 subunit of transcription elongation factor DSIF was observed at the 20E-dependent genes upon NELF depletion. We presume that an important function of NELF is to participate in stabilizing the Pol II-DSIF complex, resulting in a significant impact on transcription of its target genes. In order to directly link NELF to regulation of 20E-dependent genes in development, we show the presence of NELF at the promoters of 20E-dependent genes during their active transcription in both embryogenesis and metamorphosis. We also demonstrate that 20E-dependent promoters, while temporarily inactive at the larval stage, preserve a Pol II paused state and bind NELF complex.

Defining the divergent enzymatic properties of RNA polymerases I and II

Eukaryotes express at least three nuclear DNA-dependent RNA polymerases (Pols) responsible for synthesizing all RNA required by the cell. Despite sharing structural homology, they have functionally diverged to suit their distinct cellular roles. Although the Pols have been studied extensively, direct comparison of their enzymatic properties is difficult because studies are often conducted under disparate experimental conditions and techniques. Here, we directly compare and reveal functional differences between Saccharomyces cerevisiae Pols I and II using a series of quantitative in vitro transcription assays. We find that Pol I single-nucleotide and multinucleotide addition rate constants are faster than those of Pol II. Pol I elongation complexes are less stable than Pol II elongation complexes, and Pol I is more error prone than Pol II. Collectively, these data show that the enzymatic properties of the Pols have diverged over the course of evolution, optimizing these enzymes for their unique cellular responsibilities.

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.

On the Origin of Biomolecular Networks

Biomolecular networks have already found great utility in characterizing complex biological systems arising from pairwise interactions amongst biomolecules. Here, we explore the important and hitherto neglected role of information asymmetry in the genesis and evolution of such pairwise biomolecular interactions. Information asymmetry between sender and receiver genes is identified as a key feature distinguishing early biochemical reactions from abiotic chemistry, and a driver of network topology as biomolecular systems become more complex. In this context, we review how graph theoretical approaches can be applied not only for a better understanding of various proximate (mechanistic) relations, but also, ultimate (evolutionary) structures encoded in such networks from among all types of variations they induce. Among many possible variations, we emphasize particularly the essential role of gene duplication in terms of signaling game theory, whereby sender and receiver gene players accrue benefit from gene duplication, leading to a preferential attachment mode of network growth. The study of the resulting dynamics suggests many mathematical/computational problems, the majority of which are intractable yet yield to efficient approximation algorithms, when studied through an algebraic graph theoretic lens. We relegate for future work the role of other possible generalizations, additionally involving horizontal gene transfer, sexual recombination, endo-symbiosis, etc., which enrich the underlying graph theory even further.

Archaea, from obscurity to superhero microbes: 40 years of surprises and critical biological insights

This issue of Emerging Topics in the Life Sciences highlights current areas of research in the field of archaeal biology and the following introductory editorial sets the stage by considering some of the key developments over the last four decades since the initial identification of the archaea as a unique form of life. Emerging topics from this vibrant and rapidly expanding field of research are considered and detailed further in the articles within this issue.

G4PromFinder: an algorithm for predicting transcription promoters in GC-rich bacterial genomes based on AT-rich elements and G-quadruplex motifs

Background

Over the last few decades, computational genomics has tremendously contributed to decipher biology from genome sequences and related data. Considerable effort has been devoted to the prediction of transcription promoter and terminator sites that represent the essential “punctuation marks” for DNA transcription. Computational prediction of promoters in prokaryotes is a problem whose solution is far from being determined in computational genomics. The majority of published bacterial promoter prediction tools are based on a consensus-sequences search and they were designed specifically for vegetative σ⁷⁰ promoters and, therefore, not suitable for promoter prediction in bacteria encoding a lot of σ factors, like actinomycetes.

Results

In this study we investigated the possibility to identify putative promoters in prokaryotes based on evolutionarily conserved motifs, and focused our attention on GC-rich bacteria in which promoter prediction with conventional, consensus-based algorithms is often not-exhaustive. Here, we introduce G4PromFinder, a novel algorithm that predicts putative promoters based on AT-rich elements and G-quadruplex DNA motifs. We tested its performances by using available genomic and transcriptomic data of the model microorganisms Streptomyces coelicolor A3(2) and Pseudomonas aeruginosa PA14. We compared our results with those obtained by three currently available promoter predicting algorithms: the σ⁷⁰consensus-based PePPER, the σ factors consensus-based bTSSfinder, and PromPredict which is based on double-helix DNA stability. Our results demonstrated that G4PromFinder is more suitable than the three reference tools for both the genomes. In fact our algorithm achieved the higher accuracy (F₁-scores 0.61 and 0.53 in the two genomes) as compared to the next best tool that is PromPredict (F₁-scores 0.46 and 0.48). Consensus-based algorithms produced lower performances with the analyzed GC-rich genomes.

Conclusions

Our analysis shows that G4PromFinder is a powerful tool for promoter search in GC-rich bacteria, especially for bacteria coding for a lot of σ factors, such as the model microorganism S. coelicolor A3(2). Moreover consensus-based tools and, in general, tools that are based on specific features of bacterial σ factors seem to be less performing for promoter prediction in these types of bacterial genomes.

Electronic supplementary material

The online version of this article (10.1186/s12859-018-2049-x) contains supplementary material, which is available to authorized users.

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.

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 Mg²⁺ and also an RNA 3'-OH or a 3'-H₂ affect RNAP sliding. Because RNAP grip to template suggests a sticky translocation mechanism, and because grip is enhanced by increasing K⁺ and Mg²⁺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.

Prepartal Energy Intake Alters Blood Polymorphonuclear Leukocyte Transcriptome During the Peripartal Period in Holstein Cows

In the dairy industry, cow health and farmer profits depend on the balance between diet (ie, nutrient composition, daily intake) and metabolism. This is especially true during the transition period, where dramatic physiological changes foster vulnerability to immunosuppression, negative energy balance, and clinical and subclinical disorders. Using an Agilent microarray platform, this study examined changes in the transcriptome of bovine polymorphonuclear leukocytes (PMNLs) due to prepartal dietary intake. Holstein cows were fed a high-straw, control-energy diet (CON; NE_L = 1.34 Mcal/kg) or overfed a moderate-energy diet (OVE; NE_L = 1.62 Mcal/kg) during the dry period. Blood for PMNL isolation and metabolite analysis was collected at −14 and +7 days relative to parturition. At an analysis of variance false discovery rate <0.05, energy intake (OVE vs CON) influenced 1806 genes. Dynamic Impact Approach bioinformatics analysis classified treatment effects on Kyoto Encyclopedia of Genes and Genomes pathways, including activated oxidative phosphorylation and biosynthesis of unsaturated fatty acids and inhibited RNA polymerase, proteasome, and toll-like receptor signaling pathway. This analysis indicates that processes critical for energy metabolism and cellular and immune function were affected with mixed results. However, overall interpretation of the transcriptome data agreed in part with literature documenting a potentially detrimental, chronic activation of PMNL in response to overfeeding. The widespread, transcriptome-level changes captured here confirm the importance of dietary energy adjustments around calving on the immune system.

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.

Flexibility-rigidity index for protein-nucleic acid flexibility and fluctuation analysis

Protein-nucleic acid complexes are important for many cellular processes including the most essential functions such as transcription and translation. For many protein-nucleic acid complexes, flexibility of both macromolecules has been shown to be critical for specificity and/or function. The flexibility-rigidity index (FRI) has been proposed as an accurate and efficient approach for protein flexibility analysis. In this article, we introduce FRI for the flexibility analysis of protein-nucleic acid complexes. We demonstrate that a multiscale strategy, which incorporates multiple kernels to capture various length scales in biomolecular collective motions, is able to significantly improve the state of art in the flexibility analysis of protein-nucleic acid complexes. We take the advantage of the high accuracy and O(N) computational complexity of our multiscale FRI method to investigate the flexibility of ribosomal subunits, which are difficult to analyze by alternative approaches. An anisotropic FRI approach, which involves localized Hessian matrices, is utilized to study the translocation dynamics in an RNA polymerase.

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 identification of putative RNA polymerase II C-terminal domain associated proteins in red and green algae

A tandemly repeated C-terminal domain (CTD) of the largest subunit of RNA polymerase II is functionally essential and strongly conserved in many organisms, including animal, yeast and plant models. Although present in simple, ancestral red algae, CTD tandem repeats have undergone extensive modifications and degeneration during the evolutionary transition to developmentally complex rhodophytes. In contrast, CTD repeats are conserved in both green algae and their more complex land plant relatives. Understanding the mechanistic differences that underlie these variant patterns of CTD evolution requires knowledge of CTD-associated proteins in these 2 lineages. To provide an initial baseline comparison, we bound potential phospho-CTD associated proteins (PCAPs) to artificially synthesized and phosphorylated CTD repeats from the unicellular red alga Cyanidioschyzon merolae and green alga Chlamydomonas reinhardtii. Our results indicate that red and green algae share a number of PCAPs, including kinases and proteins involved in mRNA export. There also are important taxon-specific differences, including mRNA splicing-related PCAPs recovered from Chlamydomonas but not Cyanidioschyzon, consistent with the relative intron densities in green and red algae. Our results also offer the first experimental indication that different proteins bind 2 distinct types of repeats in Cyanidioschyzon, suggesting a division of function between the proximal and distal CTD, similar to patterns identified in more developmentally complex model organisms.