Transcription attenuation: once viewed as a novel regulatory strategy

WhiB-like proteins: Diversity of structure, function and mechanism

The WhiB-Like (Wbl) proteins are a large family of iron-sulfur (Fe-S) cluster-containing transcription factors exclusively found in the phylum Actinobacteria, including the notable genera like Mycobacteria, Streptomycetes and Corynebacteria. These proteins play pivotal roles in diverse biological processes, such as cell development, redox stress response and antibiotic resistance. Members of the Wbl family exhibit remarkable diversity in their sequences, structures and functions, attracting great attention since their first discovery. This review highlights the most recent breakthroughs in understanding the structural and mechanistic aspects of Wbl-dependent transcriptional regulation.

Multiple Forms and Functions of Premature Termination by RNA Polymerase II

Eukaryotic genomes are widely transcribed by RNA polymerase II (pol II) both within genes and in intergenic regions. POL II elongation complexes comprising the polymerase, the DNA template and nascent RNA transcript must be extremely processive in order to transcribe the longest genes which are over 1 megabase long and take many hours to traverse. Dedicated termination mechanisms are required to disrupt these highly stable complexes. Transcription termination occurs not only at the 3' ends of genes once a full length transcript has been made, but also within genes and in promiscuously transcribed intergenic regions. Termination at these latter positions is termed "premature" because it is not triggered in response to a specific signal that marks the 3' end of a gene, like a polyA site. One purpose of premature termination is to remove polymerases from intergenic regions where they are "not wanted" because they may interfere with transcription of overlapping genes or the progress of replication forks. Premature termination has recently been appreciated to occur at surprisingly high rates within genes where it is speculated to serve regulatory or quality control functions. In this review I summarize current understanding of the different mechanisms of premature termination and its potential functions.

Optimizing the production of dsRNA biocontrols in microbial systems using multiple transcriptional terminators

Crop pests and pathogens annually cause over $220 billion in global crop damage, with insects consuming 5%-20% of major grain crops. Current crop pest and disease control strategies rely on insecticidal and fungicidal sprays, plant genetic resistance, transgenes, and agricultural practices. Double-stranded RNA (dsRNA) is emerging as a novel sustainable method of plant protection as an alternative to traditional chemical pesticides. Successful commercialization of dsRNA-based biocontrols requires the economical production of large quantities of dsRNA combined with suitable delivery methods to ensure RNAi efficacy against the target pest. In this study, we have optimized the design of plasmid DNA constructs to produce dsRNA biocontrols in Escherichia coli, by employing a wide range of alternative synthetic transcriptional terminators before measurement of dsRNA yield. We demonstrate that a 7.8-fold increase of dsRNA was achieved using triple synthetic transcriptional terminators within a dual T7 dsRNA production system compared to the absence of transcriptional terminators. Moreover, our data demonstrate that batch fermentation production dsRNA using multiple transcriptional terminators is scalable and generates significantly higher yields of dsRNA generated in the absence of transcriptional terminators at both small-scale batch culture and large-scale fermentation. In addition, we show that application of these dsRNA biocontrols expressed in E. coli cells results in increased insect mortality. Finally, novel mass spectrometry analysis was performed to determine the precise sites of transcriptional termination at the different transcriptional terminators providing important further mechanistic insight.

Riboswitches

Riboswitches are structured noncoding RNA domains that are typically found embedded in messenger RNAs, where they sense specific target molecules or elemental ions and regulate gene expression. These RNAs thus serve as genetic switches that can activate or repress gene expression in response to changing levels of their target ligand. To many observers, riboswitches might seem like rare oddities that are not as sophisticated as, or competitive with, the various protein factors that perform these same roles. However, as the number of experimentally validated riboswitch classes increases, and their true biochemical sophistication is recognized, it is becoming clearer that many species from all three domains of life entrust RNAs to make important chemical sensing and gene control decisions without the necessary participation of protein factors.

A Career's Work, the l-Arabinose Operon: How It Functions and How We Learned It

Very few labs have had the good fortune to have been able to focus for more than 50 years on a relatively narrow research topic and to be in a field in which both basic knowledge and the research technology and methods have progressed as rapidly as they have in molecular biology. My research group, first at Brandeis University and then at Johns Hopkins University, has had this opportunity. In this review, therefore, I will describe largely the work from my laboratory that has spanned this period and which was carried out by 40 plus graduate students, several postdoctoral associates, my technician, and me. In addition to presenting the scientific findings or results, I will place many of the topics in scientific context and, because we needed to develop a good many of the experimental methods behind our findings, I will also describe some of these methods and their importance. Also included will be occasional comments on how the research community or my research group functioned. Because a wide variety of approaches were used throughout our work, no ideal organization of this review is apparent. Therefore, I have chosen to use a hybrid structure in which there are six sections. Within each of the sections, experiments and findings will be described roughly in chronological order. Frequent cross references between parts and sections will be made because some findings and experimental approaches could logically have been described in more than one place.

A Molecular Dynamic Model of Tryptophan Overproduction in Escherichia coli

Several deterministic models simulate the main molecular biology interactions among the numerous mechanisms controlling the dynamics of the tryptophan operon in native strains. However, no models exist to investigate bacterial tryptophan production from a biotechnological point of view. Here, we modified tryptophan models for native production to propose a biotechnological working model that incorporates the activity of tryptophan secretion systems and genetic modifications made in two reported E. coli strains. The resultant deterministic model could emulate the production of tryptophan in the same order of magnitude as those quantified experimentally by the genetically engineered E. coli strains GPT1001 and GPT1002 in shake flasks. We hope this work may contribute to the rational development of biological models that define and include the main parameters and molecular components for designing and engineering efficient biotechnological chassis to produce valuable chemicals.

Macromolecular assemblies supporting transcription-translation coupling

Coordination between the molecular machineries that synthesize and decode prokaryotic mRNAs is an important layer of gene expression control known as transcription-translation coupling. While it has long been known that translation can regulate transcription and vice-versa, recent structural and biochemical work has shed light on the underlying mechanistic basis. Complexes of RNA polymerase linked to a trailing ribosome (expressomes) have been structurally characterized in a variety of states at near-atomic resolution, and also directly visualized in cells. These data are complemented by recent biochemical and biophysical analyses of transcription-translation systems and the individual components within them. Here, we review our improved understanding of the molecular basis of transcription-translation coupling. These insights are discussed in relation to our evolving understanding of the role of coupling in cells.

Origin of translational control by eIF2α phosphorylation: insights from genome-wide translational profiling studies in fission yeast

During amino acid limitation, the protein kinase Gcn2 phosphorylates the α subunit of eIF2, thereby regulating mRNA translation. In yeast Saccharomyces cerevisiae and mammals, eIF2α phosphorylation regulates translation of related transcription factors Gcn4 and Atf4 through upstream open reading frames (uORFs) to activate transcription genome wide. However, mammals encode three more eIF2α kinases activated by distinct stimuli. Did the translational control system involving eIF2α phosphorylation evolve from so simple (as found in yeast S. cerevisiae) to complex (as found in humans)? Recent genome-wide translational profiling studies of amino acid starvation response in the fission yeast Schizosaccharomyces pombe provide an unexpected answer to this question.

Conserved DNA sequence features underlie pervasive RNA polymerase pausing

Pausing of transcribing RNA polymerase is regulated and creates opportunities to control gene expression. Research in metazoans has so far mainly focused on RNA polymerase II (Pol II) promoter-proximal pausing leaving the pervasive nature of pausing and its regulatory potential in mammalian cells unclear. Here, we developed a pause detecting algorithm (PDA) for nucleotide-resolution occupancy data and a new native elongating transcript sequencing approach, termed nested NET-seq, that strongly reduces artifactual peaks commonly misinterpreted as pausing sites. Leveraging PDA and nested NET-seq reveal widespread genome-wide Pol II pausing at single-nucleotide resolution in human cells. Notably, the majority of Pol II pauses occur outside of promoter-proximal gene regions primarily along the gene-body of transcribed genes. Sequence analysis combined with machine learning modeling reveals DNA sequence properties underlying widespread transcriptional pausing including a new pause motif. Interestingly, key sequence determinants of RNA polymerase pausing are conserved between human cells and bacteria. These studies indicate pervasive sequence-induced transcriptional pausing in human cells and the knowledge of exact pause locations implies potential functional roles in gene expression.

Simulation of Folding Kinetics for Aligned RNAs

Studying the folding kinetics of an RNA can provide insight into its function and is thus a valuable method for RNA analyses. Computational approaches to the simulation of folding kinetics suffer from the exponentially large folding space that needs to be evaluated. Here, we present a new approach that combines structure abstraction with evolutionary conservation to restrict the analysis to common parts of folding spaces of related RNAs. The resulting algorithm can recapitulate the folding kinetics known for single RNAs and is able to analyse even long RNAs in reasonable time. Our program RNAliHiKinetics is the first algorithm for the simulation of consensus folding kinetics and addresses a long-standing problem in a new and unique way.

Genetic circuit characterization by inferring RNA polymerase movement and ribosome usage

To perform their computational function, genetic circuits change states through a symphony of genetic parts that turn regulator expression on and off. Debugging is frustrated by an inability to characterize parts in the context of the circuit and identify the origins of failures. Here, we take snapshots of a large genetic circuit in different states: RNA-seq is used to visualize circuit function as a changing pattern of RNA polymerase (RNAP) flux along the DNA. Together with ribosome profiling, all 54 genetic parts (promoters, ribozymes, RBSs, terminators) are parameterized and used to inform a mathematical model that can predict circuit performance, dynamics, and robustness. The circuit behaves as designed; however, it is riddled with genetic errors, including cryptic sense/antisense promoters and translation, attenuation, incorrect start codons, and a failed gate. While not impacting the expected Boolean logic, they reduce the prediction accuracy and could lead to failures when the parts are used in other designs. Finally, the cellular power (RNAP and ribosome usage) required to maintain a circuit state is calculated. This work demonstrates the use of a small number of measurements to fully parameterize a regulatory circuit and quantify its impact on host.

Gcn2 eIF2α kinase mediates combinatorial translational regulation through nucleotide motifs and uORFs in target mRNAs

The protein kinase Gcn2 is a central transducer of nutritional stress signaling important for stress adaptation by normal cells and the survival of cancer cells. In response to nutrient deprivation, Gcn2 phosphorylates eIF2α, thereby repressing general translation while enhancing translation of specific mRNAs with upstream ORFs (uORFs) situated in their 5'-leader regions. Here we performed genome-wide measurements of mRNA translation during histidine starvation in fission yeast Schizosaccharomyces pombe. Polysome analyses were combined with microarray measurements to identify gene transcripts whose translation was up-regulated in response to the stress in a Gcn2-dependent manner. We determined that translation is reprogrammed to enhance RNA metabolism and chromatin regulation and repress ribosome synthesis. Interestingly, translation of intron-containing mRNAs was up-regulated. The products of the regulated genes include additional eIF2α kinase Hri2 amplifying the stress signaling and Gcn5 histone acetyl transferase and transcription factors, together altering genome-wide transcription. Unique dipeptide-coding uORFs and nucleotide motifs, such as '5'-UGA(C/G)GG-3', are found in 5' leader regions of regulated genes and shown to be responsible for translational control.

RNA Polymerase's Relationship with the Ribosome: Not So Physical, Most of the Time

In bacteria, the rates of transcription elongation and translation elongation are coordinated, changing together in response to growth conditions. It has been proposed that this is due to physical coupling of RNA polymerase and the lead ribosome on nascent mRNA, an interaction important for preventing premature transcription termination by Rho factor. Recent studies challenge this view and provide evidence that coordination is indirect, mediated in Escherichia coli by the alarmone (p)ppGpp. Here, we discuss these new findings and how they shape our understanding of the functional relationship between RNA polymerase and the ribosome as well as the basis of transcriptional polarity.

Automated Design of Diverse Stand-Alone Riboswitches

Riboswitches that couple binding of ligands to conformational changes offer sensors and control elements for RNA synthetic biology and medical biotechnology. However, design of these riboswitches has required expert intuition or software specialized to transcription or translation outputs; design has been particularly challenging for applications in which the riboswitch output cannot be amplified by other molecular machinery. We present a fully automated design method called RiboLogic for such “stand-alone” riboswitches and test it via high-throughput experiments on 2875 molecules using RNA-MaP (RNA on a massively parallel array) technology. These molecules consistently modulate their affinity to the MS2 bacteriophage coat protein upon binding of flavin mononucleotide, tryptophan, theophylline, and microRNA miR-208a, achieving activation ratios of up to 20 and significantly better performance than control designs. By encompassing a wide diversity of stand-alone switches and highly quantitative data, the resulting ribologic-solves experimental data set provides a rich resource for further improvement of riboswitch models and design methods.

The ancient alarmone ZTP and zinc homeostasis in Bacillus subtilis

In Bacillus subtilis a sophisticated regulatory circuit that involves Z nucleoside triphosphate (ZTP) is recruited to optimize cellular zinc distribution when cytoplasmic zinc is scarce. This process uses enzymatic reactions to measure the pool of available zinc ions and amplifies this signal to control the activity of zinc chaperones. The ZTP-dependent regulatory circuit that is exploited for zinc homeostasis controls purine and folate biosynthesis, which starts with GTP as initial substrate. Low concentrations of formyl-tetrahydrofolate (fTHF) lead to accumulation of the intermediate 5'-phosphoribosyl-4-carboxyamide-5-aminoimidazole (AICAR or ZMP), which is pyrophosphorylated by another intermediate to ZTP. This alarmone activates expression of genes using a ZTP-dependent riboswitch in many bacterial strains. In this way, the cellular folate concentration controls folate biosynthesis via the enzymatic activity of the fTHF-dependent AICAR-transforming reaction. Zinc distribution control is layered onto this circuit. The 'sensor' is the activity of the initial reaction of folate synthesis from GTP, which is catalyzed by a zinc-dependent enzyme FolE_IA or its metal-cambialistic paralog FolE_IB . Consequently, low zinc lowers folate levels, causing AICAR accumulation and ZTP formation. In addition to the riboswitch, ZTP activates the zinc chaperone ZagA of the COG0523 protein family, which efficiently allocate zinc to zinc-dependent enzymes such as FolE_IA .

Two Old Dogs, One New Trick: A Review of RNA Polymerase and Ribosome Interactions during Transcription-Translation Coupling

The coupling of transcription and translation is more than mere translation of an mRNA that is still being transcribed. The discovery of physical interactions between RNA polymerase and ribosomes has spurred renewed interest into this long-standing paradigm of bacterial molecular biology. Here, we provide a concise presentation of recent insights gained from super-resolution microscopy, biochemical, and structural work, including cryo-EM studies. Based on the presented data, we put forward a dynamic model for the interaction between RNA polymerase and ribosomes, in which the interactions are repeatedly formed and broken. Furthermore, we propose that long intervening nascent RNA will loop out and away during the forming the interactions between the RNA polymerase and ribosomes. By comparing the effect of the direct interactions between RNA polymerase and ribosomes with those that transcription factors NusG and RfaH mediate, we submit that two distinct modes of coupling exist: Factor-free and factor-mediated coupling. Finally, we provide a possible framework for transcription-translation coupling and elude to some open questions in the field.

Ribosome-induced RNA conformational changes in a viral 3'-UTR sense and regulate translation levels

Structured RNA elements, programmed RNA conformational changes, and interactions between different RNA domains underlie many modes of regulating gene expression, mandating studies to understand the foundational principles that govern these phenomena. Exploring the structured 3' untranslated region (UTR) of a viral RNA, we discovered that different contexts of the 3'-UTR confer different abilities to enhance translation of an associated open reading frame. In one context, ribosome-induced conformational changes in a 'sensor' RNA domain affect a separate RNA 'functional' domain, altering translation efficiency. The structure of the entire 3'-UTR reveals that structurally distinct domains use a spine of continuously stacked bases and a strut-like linker to create a conduit for communication within the higher-order architecture. Thus, this 3'-UTR RNA illustrates how RNA can use programmed conformational changes to sense the translation status of an upstream open reading frame, then create a tuned functional response by communicating that information to other RNA elements.

Riboswitches and Translation Control

A growing collection of bacterial riboswitch classes is being discovered that sense central metabolites, coenzymes, and signaling molecules. Included among the various mechanisms of gene regulation exploited by these RNA regulatory elements are several that modulate messenger RNA (mRNA) translation. In this review, the mechanisms of riboswitch-mediated translation control are summarized to highlight both their diversity and potential ancient origins. These mechanisms include ligand-gated presentation or occlusion of ribosome-binding sites, control of alternative splicing of mRNAs, and the regulation of mRNA stability. Moreover, speculation on the potential for novel riboswitch discoveries is presented, including a discussion on the potential for the discovery of a greater diversity of mechanisms for translation control.

Measures of single- versus multiple-round translation argue against a mechanism to ensure coupling of transcription and translation

In prokaryotes, the synthesis of RNA and protein occurs simultaneously in the cytoplasm. A number of studies indicate that translation can strongly impact transcription, a phenomenon often attributed to physical coupling between RNA polymerase (RNAP) and the lead ribosome on the nascent mRNA. Whether there generally exists a mechanism to ensure or promote RNAP-ribosome coupling remains unclear. Here, we used an efficient hammerhead ribozyme and developed a reporter system to measure single- versus multiple-round translation in Escherichia coli Six pairs of cotranscribed and differentially translated genes were analyzed. For five of them, the stoichiometry of the two protein products came no closer to unity (1:1) when the rounds of translation were severely reduced in wild-type cells. Introduction of mutation rpoB(I572N), which slows RNAP elongation, could promote coupling, as indicated by stoichiometric SspA and SspB products in the single-round assay. These data are consistent with models of stochastic coupling in which the probability of coupling depends on the relative rates of transcription and translation and suggest that RNAP often transcribes without a linked ribosome.

Transcription elongation rate affects nascent histone pre-mRNA folding and 3' end processing

In this study, Saldi et al. investigated how transcription elongation rate influences cotranscriptional pre-mRNA maturation. Their findings show that regulation of transcription speed can modulate pre-mRNA processing by changing nascent RNA structure and suggest a mechanism by which alternative processing could be controlled.

Transcription elongation rate influences cotranscriptional pre-mRNA maturation, but how such kinetic coupling works is poorly understood. The formation of nonadenylated histone mRNA 3′ ends requires recognition of an RNA structure by stem–loop-binding protein (SLBP). We report that slow transcription by mutant RNA polymerase II (Pol II) caused accumulation of polyadenylated histone mRNAs that extend past the stem–loop processing site. UV irradiation, which decelerates Pol II elongation, also induced long poly(A)⁺ histone transcripts. Inhibition of 3′ processing by slow Pol II correlates with failure to recruit SLBP to histone genes. Chemical probing of nascent RNA structure showed that the stem–loop fails to fold in transcripts made by slow Pol II, thereby explaining the absence of SLBP and failure to process 3′ ends. These results show that regulation of transcription speed can modulate pre-mRNA processing by changing nascent RNA structure and suggest a mechanism by which alternative processing could be controlled.

Allosteric mechanism of the V. vulnificus adenine riboswitch resolved by four-dimensional chemical mapping

The structural interconversions that mediate the gene regulatory functions of RNA molecules may be different from classic models of allostery, but the relevant structural correlations have remained elusive in even intensively studied systems. Here, we present a four-dimensional expansion of chemical mapping called lock-mutate-map-rescue (LM²R), which integrates multiple layers of mutation with nucleotide-resolution chemical mapping. This technique resolves the core mechanism of the adenine-responsive V. vulnificus add riboswitch, a paradigmatic system for which both Monod-Wyman-Changeux (MWC) conformational selection models and non-MWC alternatives have been proposed. To discriminate amongst these models, we locked each functionally important helix through designed mutations and assessed formation or depletion of other helices via compensatory rescue evaluated by chemical mapping. These LM²R measurements give strong support to the pre-existing correlations predicted by MWC models, disfavor alternative models, and suggest additional structural heterogeneities that may be general across ligand-free riboswitches.

Transcription-translation coupling: direct interactions of RNA polymerase with ribosomes and ribosomal subunits

In prokaryotes, RNA polymerase and ribosomes can bind concurrently to the same RNA transcript, leading to the functional coupling of transcription and translation. The interactions between RNA polymerase and ribosomes are crucial for the coordination of transcription with translation. Here, we report that RNA polymerase directly binds ribosomes and isolated large and small ribosomal subunits. RNA polymerase and ribosomes form a one-to-one complex with a micromolar dissociation constant. The formation of the complex is modulated by the conformational and functional states of RNA polymerase and the ribosome. The binding interface on the large ribosomal subunit is buried by the small subunit during protein synthesis, whereas that on the small subunit remains solvent-accessible. The RNA polymerase binding site on the ribosome includes that of the isolated small ribosomal subunit. This direct interaction between RNA polymerase and ribosomes may contribute to the coupling of transcription to translation.

Long-Range Interactions in Riboswitch Control of Gene Expression

Riboswitches are widespread RNA motifs that regulate gene expression in response to fluctuating metabolite concentrations. Known primarily from bacteria, riboswitches couple specific ligand binding and changes in RNA structure to mRNA expression in cis. Crystal structures of the ligand binding domains of most of the phylogenetically widespread classes of riboswitches, each specific to a particular metabolite or ion, are now available. Thus, the bound states-one end point-have been thoroughly characterized, but the unbound states have been more elusive. Consequently, it is less clear how the unbound, sensing riboswitch refolds into the ligand binding-induced output state. The ligand recognition mechanisms of riboswitches are diverse, but we find that they share a common structural strategy in positioning their binding sites at the point of the RNA three-dimensional fold where the residues farthest from one another in sequence meet. We review how riboswitch folds adhere to this fundamental strategy and propose future research directions for understanding and harnessing their ability to specifically control gene expression.

Global analysis of translation termination in E. coli

Terminating protein translation accurately and efficiently is critical for both protein fidelity and ribosome recycling for continued translation. The three bacterial release factors (RFs) play key roles: RF1 and 2 recognize stop codons and terminate translation; and RF3 promotes disassociation of bound release factors. Probing release factors mutations with reporter constructs containing programmed frameshifting sequences or premature stop codons had revealed a propensity for readthrough or frameshifting at these specific sites, but their effects on translation genome-wide have not been examined. We performed ribosome profiling on a set of isogenic strains with well-characterized release factor mutations to determine how they alter translation globally. Consistent with their known defects, strains with increasingly severe release factor defects exhibit increasingly severe accumulation of ribosomes over stop codons, indicative of an increased duration of the termination/release phase of translation. Release factor mutant strains also exhibit increased occupancy in the region following the stop codon at a significant number of genes. Our global analysis revealed that, as expected, translation termination is generally efficient and accurate, but that at a significant number of genes (≥ 50) the ribosome signature after the stop codon is suggestive of translation past the stop codon. Even native E. coli K-12 exhibits the ribosome signature suggestive of protein extension, especially at UGA codons, which rely exclusively on the reduced function RF2 variant of the K-12 strain for termination. Deletion of RF3 increases the severity of the defect. We unambiguously demonstrate readthrough and frameshifting protein extensions and their further accumulation in mutant strains for a few select cases. In addition to enhancing recoding, ribosome accumulation over stop codons disrupts attenuation control of biosynthetic operons, and may alter expression of some overlapping genes. Together, these functional alterations may either augment the protein repertoire or produce deleterious proteins.

Proteins are the cellular workhorses, performing essentially all of the functions required for cell and organismal survival. But, it takes a great deal of energy to make proteins, making it critical that proteins are made accurately and in the proper time frame. After a ribosome synthesizes a protein, release factors catalyze the accurate and timely release of the finished protein from the ribosome, a process called termination. Ribosomes are then recycled and start the next protein. We utilized ribosome profiling, a method that allows us to follow the position of every ribosome that is making a protein, to globally investigate and strengthen insights on termination fidelity for cells with and without mutant release factors. We find that as we decrease release factor function, the time to terminate/release a protein increases across the genome. We observe that the accuracy of terminating a protein at the correct place decreases on a global scale. Using this metric we identify genes with inherently low termination efficiency and confirm two novel events resulting in extended protein products. In addition we find that beyond disrupting accurate protein synthesis, release factor mutations can alter expression of genes involved in the production of key amino acids.

Insights into the Mechanisms of Basal Coordination of Transcription Using a Genome-Reduced Bacterium

Coordination of transcription in bacteria occurs at supra-operonic scales, but the extent, specificity, and mechanisms of such regulation are poorly understood. Here, we tackle this problem by profiling the transcriptome of the model organism Mycoplasma pneumoniae across 115 growth conditions. We identify three qualitatively different levels of co-expression corresponding to distinct relative orientations and intergenic properties of adjacent genes. We reveal that the degree of co-expression between co-directional adjacent operons, and more generally between genes, is tightly related to their capacity to be transcribed en bloc into the same mRNA. We further show that this genome-wide pervasive transcription of adjacent genes and operons is specifically repressed by DNA regions preferentially bound by RNA polymerases, by intrinsic terminators, and by large intergenic distances. Taken together, our findings suggest that the basal coordination of transcription is mediated by the physical entities and mechanical properties of the transcription process itself, and that operon-like behaviors may strongly vary from condition to condition.

Terminator Operon Reporter: combining a transcription termination switch with reporter technology for improved gene synthesis and synthetic biology applications

Synthetic biology is characterized by the development of novel and powerful DNA fabrication methods and by the application of engineering principles to biology. The current study describes Terminator Operon Reporter (TOR), a new gene assembly technology based on the conditional activation of a reporter gene in response to sequence errors occurring at the assembly stage of the synthetic element. These errors are monitored by a transcription terminator that is placed between the synthetic gene and reporter gene. Switching of this terminator between active and inactive states dictates the transcription status of the downstream reporter gene to provide a rapid and facile readout of the accuracy of synthetic assembly. Designed specifically and uniquely for the synthesis of protein coding genes in bacteria, TOR allows the rapid and cost-effective fabrication of synthetic constructs by employing oligonucleotides at the most basic purification level (desalted) and without the need for costly and time-consuming post-synthesis correction methods. Thus, TOR streamlines gene assembly approaches, which are central to the future development of synthetic biology.

Coupling of RNA Polymerase II Transcription Elongation with Pre-mRNA Splicing

Pre-mRNA maturation frequently occurs at the same time and place as transcription by RNA polymerase II. The co-transcriptionality of mRNA processing has permitted the evolution of mechanisms that functionally couple transcription elongation with diverse events that occur on the nascent RNA. This review summarizes the current understanding of the relationship between transcriptional elongation through a chromatin template and co-transcriptional splicing including alternative splicing decisions that affect the expression of most human genes.

Biosynthesis and Regulation of the Branched-Chain Amino Acids†

This review focuses on more recent studies concerning the systems biology of branched-chain amino acid biosynthesis, that is, the pathway-specific and global metabolic and genetic regulatory networks that enable the cell to adjust branched-chain amino acid synthesis rates to changing nutritional and environmental conditions. It begins with an overview of the enzymatic steps and metabolic regulatory mechanisms of the pathways and descriptions of the genetic regulatory mechanisms of the individual operons of the isoleucine-leucine-valine (ilv) regulon. This is followed by more-detailed discussions of recent evidence that global control mechanisms that coordinate the expression of the operons of this regulon with one another and the growth conditions of the cell are mediated by changes in DNA supercoiling that occur in response to changes in cellular energy charge levels that, in turn, are modulated by nutrient and environmental signals. Since the parallel pathways for isoleucine and valine biosynthesis are catalyzed by a single set of enzymes, and because the AHAS-catalyzed reaction is the first step specific for valine biosynthesis but the second step of isoleucine biosynthesis, valine inhibition of a single enzyme for this enzymatic step might compromise the cell for isoleucine or result in the accumulation of toxic intermediates. The operon-specific regulatory mechanisms of the operons of the ilv regulon are discussed in the review followed by a consideration and brief review of global regulatory proteins such as integration host factor (IHF), Lrp, and CAP (CRP) that affect the expression of these operons.

Mathematical modeling of the apo and holo transcriptional regulation in Escherichia coli

Transcription factors can bind to DNA either with their effector bound (holo conformation), or as free proteins (apo conformation).

Synergetic regulation of translational reading-frame switch by ligand-responsive RNAs in mammalian cells

Distinct translational initiation mechanisms between prokaryotes and eukaryotes limit the exploitation of prokaryotic riboswitch repertoire for regulatory RNA circuit construction in mammalian application. Here, we explored programmed ribosomal frameshifting (PRF) as the regulatory gene expression platform for engineered ligand-responsive RNA devices in higher eukaryotes. Regulation was enabled by designed ligand-dependent conformational rearrangements of the two cis-acting RNA motifs of opposite activity in -1 PRF. Particularly, RNA elements responsive to trans-acting ligands can be tailored to modify co-translational RNA refolding dynamics of a hairpin upstream of frameshifting site to achieve reversible and adjustable -1 PRF attenuating activity. Combined with a ligand-responsive stimulator, synthetic RNA devices for synergetic translational-elongation control of gene expression can be constructed. Due to the similarity between co-transcriptional RNA hairpin folding and co-translational RNA hairpin refolding, the RNA-responsive ligand repertoire provided in prokaryotic systems thus becomes accessible to gene-regulatory circuit construction for synthetic biology application in mammalian cells.

Novel HIV-1 miRNAs stimulate TNFα release in human macrophages via TLR8 signaling pathway

Purpose

To determine whether HIV-1 produces microRNAs and elucidate whether these miRNAs can induce inflammatory response in macrophages (independent of the conventional miRNA function in RNA interference) leading to chronic immune activation.

Methods

Using sensitive quantitative Real Time RT-PCR and sequencing, we detected novel HIV-derived miRNAs in the sera of HIV+ persons, and associated with exosomes. Release of TNFα by macrophages challenged with HIV miRNAs was measured by ELISA.

Results

HIV infection of primary alveolar macrophages produced elevated levels of viral microRNAs vmiR88, vmiR99 and vmiR-TAR in cell extracts and in exosome preparations from conditioned medium. Furthermore, these miRNAs were also detected in exosome fraction of sera from HIV-infected persons. Importantly, vmiR88 and vmiR99 (but not vmiR-TAR) stimulated human macrophage TNFα release, which is dependent on macrophage TLR8 expression. These data support a potential role for HIV-derived vmiRNAs released from infected macrophages as contributing to chronic immune activation in HIV-infected persons, and may represent a novel therapeutic target to limit AIDS pathogenesis.

Conclusion

Novel HIV vmiR88 and vmiR99 are present in the systemic circulation of HIV+ persons and could exhibit biological function (independent of gene silencing) as ligands for TLR8 signaling that promote macrophage TNFα release, and may contribute to chronic immune activation. Targeting novel HIV-derived miRNAs may represent a therapeutic strategy to limit chronic immune activation and AIDS progression.

Optimal performance of the tryptophan operon of E. coli: a stochastic, dynamical, mathematical-modeling approach

In this work, we develop a detailed, stochastic, dynamical model for the tryptophan operon of E. coli, and estimate all of the model parameters from reported experimental data. We further employ the model to study the system performance, considering the amount of biochemical noise in the trp level, the system rise time after a nutritional shift, and the amount of repressor molecules necessary to maintain an adequate level of repression, as indicators of the system performance regime. We demonstrate that the level of cooperativity between repressor molecules bound to the first two operators in the trp promoter affects all of the above enlisted performance characteristics. Moreover, the cooperativity level found in the wild-type bacterial strain optimizes a cost-benefit function involving low biochemical noise in the tryptophan level, short rise time after a nutritional shift, and low number of regulatory molecules.

Regulation of gene expression by macrolide-induced ribosomal frameshifting

The expression of many genes is controlled by upstream ORFs (uORFs). Typically, the progression of the ribosome through a regulatory uORF, which depends on the physiological state of the cell, influences the expression of the downstream gene. In the classic mechanism of induction of macrolide resistance genes, antibiotics promote translation arrest within the uORF, and the static ribosome induces a conformational change in mRNA, resulting in the activation of translation of the resistance cistron. We show that ketolide antibiotics, which do not induce ribosome stalling at the uORF of the ermC resistance gene, trigger its expression via a unique mechanism. Ketolides promote frameshifting at the uORF, allowing the translating ribosome to invade the intergenic spacer. The dynamic unfolding of the mRNA structure leads to the activation of resistance. Conceptually similar mechanisms may control other cellular genes. The identified property of ketolides to reduce the fidelity of reading frame maintenance may have medical implications.

Measurement and modeling of intrinsic transcription terminators

The reliable forward engineering of genetic systems remains limited by the ad hoc reuse of many types of basic genetic elements. Although a few intrinsic prokaryotic transcription terminators are used routinely, termination efficiencies have not been studied systematically. Here, we developed and validated a genetic architecture that enables reliable measurement of termination efficiencies. We then assembled a collection of 61 natural and synthetic terminators that collectively encode termination efficiencies across an ∼800-fold dynamic range within Escherichia coli. We simulated co-transcriptional RNA folding dynamics to identify competing secondary structures that might interfere with terminator folding kinetics or impact termination activity. We found that structures extending beyond the core terminator stem are likely to increase terminator activity. By excluding terminators encoding such context-confounding elements, we were able to develop a linear sequence-function model that can be used to estimate termination efficiencies (r = 0.9, n = 31) better than models trained on all terminators (r = 0.67, n = 54). The resulting systematically measured collection of terminators should improve the engineering of synthetic genetic systems and also advance quantitative modeling of transcription termination.

Bacterial translational regulations: high diversity between all mRNAs and major role in gene expression

BACKGROUND

In bacteria, the weak correlations at the genome scale between mRNA and protein levels suggest that not all mRNAs are translated with the same efficiency. To experimentally explore mRNA translational level regulation at the systemic level, the detailed translational status (translatome) of all mRNAs was measured in the model bacterium Lactococcus lactis in exponential phase growth.

RESULTS

Results demonstrated that only part of the entire population of each mRNA species was engaged in translation. For transcripts involved in translation, the polysome size reached a maximum of 18 ribosomes. The fraction of mRNA engaged in translation (ribosome occupancy) and ribosome density were not constant for all genes. This high degree of variability was analyzed by bioinformatics and statistical modeling in order to identify general rules of translational regulation. For most of the genes, the ribosome density was lower than the maximum value revealing major control of translation by initiation. Gene function was a major translational regulatory determinant. Both ribosome occupancy and ribosome density were particularly high for transcriptional regulators, demonstrating the positive role of translational regulation in the coordination of transcriptional networks. mRNA stability was a negative regulatory factor of ribosome occupancy and ribosome density, suggesting antagonistic regulation of translation and mRNA stability. Furthermore, ribosome occupancy was identified as a key component of intracellular protein levels underlining the importance of translational regulation.

CONCLUSIONS

We have determined, for the first time in a bacterium, the detailed translational status for all mRNAs present in the cell. We have demonstrated experimentally the high diversity of translational states allowing individual gene differentiation and the importance of translation-level regulation in the complex process linking gene expression to protein synthesis.

mRNA decapping factors and the exonuclease Xrn2 function in widespread premature termination of RNA polymerase II transcription

We report a function of human mRNA decapping factors in control of transcription by RNA polymerase II. Decapping proteins Edc3, Dcp1a, and Dcp2 and the termination factor TTF2 coimmunoprecipitate with Xrn2, the nuclear 5'-3' exonuclease "torpedo" that facilitates transcription termination at the 3' ends of genes. Dcp1a, Xrn2, and TTF2 localize near transcription start sites (TSSs) by ChIP-seq. At genes with 5' peaks of paused pol II, knockdown of decapping or termination factors Xrn2 and TTF2 shifted polymerase away from the TSS toward upstream and downstream distal positions. This redistribution of pol II is similar in magnitude to that caused by depletion of the elongation factor Spt5. We propose that coupled decapping of nascent transcripts and premature termination by the "torpedo" mechanism is a widespread mechanism that limits bidirectional pol II elongation. Regulated cotranscriptional decapping near promoter-proximal pause sites followed by premature termination could control productive pol II elongation.

Control of Transcriptional Elongation by RNA Polymerase II: A Retrospective

The origins of our current understanding of control of transcription elongation lie in pioneering experiments that mapped RNA polymerase II on viral and cellular genes. These studies first uncovered the surprising excess of polymerase molecules that we now know to be situated at the at the 5' ends of most genes in multicellular organisms. The pileup of pol II near transcription start sites reflects a ubiquitous bottle-neck that limits elongation right at the start of the transcription elongation. Subsequent seminal work identified conserved protein factors that positively and negatively control the flux of polymerase through this bottle-neck, and make a major contribution to control of gene expression.

Identification of prokaryotic small proteins using a comparative genomic approach

MOTIVATION

Accurate prediction of genes encoding small proteins (on the order of 50 amino acids or less) remains an elusive open problem in bioinformatics. Some of the best methods for gene prediction use either sequence composition analysis or sequence similarity to a known protein coding sequence. These methods often fail for small proteins, however, either due to a lack of experimentally verified small protein coding genes or due to the limited statistical significance of statistics on small sequences. Our approach is based upon the hypothesis that true small proteins will be under selective pressure for encoding the particular amino acid sequence, for ease of translation by the ribosome and for structural stability. This stability can be achieved either independently or as part of a larger protein complex. Given this assumption, it follows that small proteins should display conserved local protein structure properties much like larger proteins. Our method incorporates neural-net predictions for three local structure alphabets within a comparative genomic approach using a genomic alignment of 22 closely related bacteria genomes to generate predictions for whether or not a given open reading frame (ORF) encodes for a small protein.

RESULTS

We have applied this method to the complete genome for Escherichia coli strain K12 and looked at how well our method performed on a set of 60 experimentally verified small proteins from this organism. Out of a total of 11 407 possible ORFs, we found that 6 of the top 10 and 27 of the top 100 predictions belonged to the set of 60 experimentally verified small proteins. We found 35 of all the true small proteins within the top 200 predictions. We compared our method to Glimmer, using a default Glimmer protocol and a modified small ORF Glimmer protocol with a lower minimum size cutoff. The default Glimmer protocol identified 16 of the true small proteins (all in the top 200 predictions), but failed to predict on 34 due to size cutoffs. The small ORF Glimmer protocol made predictions for all the experimentally verified small proteins but only contained 9 of the 60 true small proteins within the top 200 predictions.

CONTACT

jsamayoa@jhu.edu

TRAP binding to the Bacillus subtilis trp leader region RNA causes efficient transcription termination at a weak intrinsic terminator

The Bacillus subtilis trpEDCFBA operon is regulated by a transcription attenuation mechanism controlled by the trp RNA-binding attenuation protein (TRAP). TRAP binds to 11 (G/U)AG repeats in the trp leader transcript and prevents formation of an antiterminator, which allows formation of an intrinsic terminator (attenuator). Previously, formation of the attenuator RNA structure was believed to be solely responsible for signaling RNA polymerase (RNAP) to halt transcription. However, base substitutions that prevent formation of the antiterminator, and thus allow the attenuator structure to form constitutively, do not result in efficient transcription termination. The observation that the attenuator requires the presence of TRAP bound to the nascent RNA to cause efficient transcription termination suggests TRAP has an additional role in causing termination at the attenuator. We show that the trp attenuator is a weak intrinsic terminator due to low GC content of the hairpin stem and interruptions in the U-stretch following the hairpin. We also provide evidence that termination at the trp attenuator requires forward translocation of RNA polymerase and that TRAP binding to the nascent transcript can induce this activity.

BarMap: RNA folding on dynamic energy landscapes

Dynamical changes of RNA secondary structures play an important role in the function of many regulatory RNAs. Such kinetic effects, especially in time-variable and externally triggered systems, are usually investigated by means of extensive and expensive simulations of large sets of individual folding trajectories. Here we describe the theoretical foundations of a generic approach that not only allows the direct computation of approximate population densities but also reduces the efforts required to analyze the folding energy landscapes to a one-time preprocessing step. The basic idea is to consider the kinetics on individual landscapes and to model external triggers and environmental changes as small but discrete changes in the landscapes. A "barmap" links macrostates of temporally adjacent landscapes and defines the transfer of population densities from one "snapshot" to the next. Implemented in the BarMap software, this approach makes it feasible to study folding processes at the level of basins, saddle points, and barriers for many nonstationary scenarios, including temperature changes, cotranscriptional folding, refolding in consequence to degradation, and mechanically constrained kinetics, as in the case of the translocation of a polymer through a pore.

TRANSAT-- method for detecting the conserved helices of functional RNA structures, including transient, pseudo-knotted and alternative structures

The prediction of functional RNA structures has attracted increased interest, as it allows us to study the potential functional roles of many genes. RNA structure prediction methods, however, assume that there is a unique functional RNA structure and also do not predict functional features required for in vivo folding. In order to understand how functional RNA structures form in vivo, we require sophisticated experiments or reliable prediction methods. So far, there exist only a few, experimentally validated transient RNA structures. On the computational side, there exist several computer programs which aim to predict the co-transcriptional folding pathway in vivo, but these make a range of simplifying assumptions and do not capture all features known to influence RNA folding in vivo. We want to investigate if evolutionarily related RNA genes fold in a similar way in vivo. To this end, we have developed a new computational method, Transat, which detects conserved helices of high statistical significance. We introduce the method, present a comprehensive performance evaluation and show that Transat is able to predict the structural features of known reference structures including pseudo-knotted ones as well as those of known alternative structural configurations. Transat can also identify unstructured sub-sequences bound by other molecules and provides evidence for new helices which may define folding pathways, supporting the notion that homologous RNA sequence not only assume a similar reference RNA structure, but also fold similarly. Finally, we show that the structural features predicted by Transat differ from those assuming thermodynamic equilibrium. Unlike the existing methods for predicting folding pathways, our method works in a comparative way. This has the disadvantage of not being able to predict features as function of time, but has the considerable advantage of highlighting conserved features and of not requiring a detailed knowledge of the cellular environment.

Many non-coding genes exert their function via an RNA structure which starts emerging while the RNA sequence is being transcribed from the genome. The resulting folding pathway is known to depend on a variety of features such as the transcription speed, the concentration of various ions and the binding of proteins and other molecules. Not all of these influences can be adequately captured by the existing computational methods which try to replicate what happens in vivo. So far, it has been challenging to experimentally investigate co-transcriptional folding pathways in vivo and only little data from in vitro experiments exists. In order to investigate if functionally similar RNA sequences from different organisms fold in a similar way, we have developed a new computational method, called Transat, which does not require the detailed computational modeling of the cellular environment. We show in a comprehensive analysis that our method is capable of detecting known structural features and provide evidence that structural features of the in vivo folding pathways have been conserved for several biologically interesting classes of RNA sequences.

Functional genomics via metabolic footprinting: monitoring metabolite secretion by Escherichia coli tryptophan metabolism mutants using FT-IR and direct injection electrospray mass spectrometry

We sought to test the hypothesis that mutant bacterial strains could be discriminated from each other on the basis of the metabolites they secrete into the medium (their ‘metabolic footprint’), using two methods of ‘global’ metabolite analysis (FT–IR and direct injection electrospray mass spectrometry). The biological system used was based on a published study of Escherichia coli tryptophan mutants that had been analysed and discriminated by Yanofsky and colleagues using transcriptome analysis. Wild-type strains supplemented with tryptophan or analogues could be discriminated from controls using FT–IR of 24 h broths, as could each of the mutant strains in both minimal and supplemented media. Direct injection electrospray mass spectrometry with unit mass resolution could also be used to discriminate the strains from each other, and had the advantage that the discrimination required the use of just two or three masses in each case. These were determined via a genetic algorithm. Both methods are rapid, reagentless, reproducible and cheap, and might beneficially be extended to the analysis of gene knockout libraries.

Cycling expression and cooperative operator interaction in the trp operon of Escherichia coli

Oscillatory behaviour in the tryptophan operon of an Escherichia coli mutant strain lacking the enzyme-inhibition regulatory mechanism has been observed by Bliss et al. but not confirmed by others. This behaviour could be important from the standpoint of synthetic biology, whose goals include the engineering of intracellular genetic oscillators. This work is devoted to investigating, from a mathematical modelling point of view, the possibility that the trp operon of the E. coli inhibition-free strain expresses cyclically. For that we extend a previously introduced model for the regulatory pathway of the tryptophan operon in Escherichia coli to account for the observed multiplicity and cooperativity of repressor binding sites. Thereafter we investigate the model dynamics using deterministic numeric solutions, stochastic simulations, and analytic studies. Our results suggest that a quasi-periodic behaviour could be observed in the trp operon expression level of single bacteria.

Translational control of the antibiotic inducibility of the PA5471 gene required for mexXY multidrug efflux gene expression in Pseudomonas aeruginosa

The PA5471 gene required for induction of the MexXY multidrug efflux system in response to ribosome-targeting antimicrobials was itself shown to be inducible by ribosome-targeting antimicrobials (Y. Morita, M. L. Sobel, and K. Poole, J. Bacteriol. 188:1847-1855, 2006). Using a lacZ transcriptional reporter, drug inducibility of PA5471 was shown to require the entirety of the 367-bp PA5472-PA5471 intergenic region. A constitutive promoter activity was, however, localized to the first 75 bp of this region, within which a single PA5471 transcription initiation site was mapped. That 3' sequences of the intergenic region blocked PA5471 expression and made it antibiotic dependent was suggestive of an attenuation mechanism of control. A 13-amino-acid leader peptide (LP)-encoding open reading frame preceded by a Shine-Dalgarno sequence was identified ca. 250 bp upstream of the PA5471 coding sequence, and its expression and translation were confirmed using a lacZ translational reporter. Alteration of the initiation codon (M1T) or introduction of translational stop signals at codons 3 (Q3Am) and 8 (C8Op) of this LP sequence (PA5471.1) yielded high-level constitutive expression of PA5471, suggesting that interference with LP translation was linked to PA5471 gene expression. Consistent with this, a Q3K mutation in the LP sequence maintained the drug inducibility of PA5471 expression. Introduction of the LP Q3Am mutation into the chromosome of Pseudomonas aeruginosa yielded stronger expression of PA5471 than did antibiotic (chloramphenicol) exposure of wild-type P. aeruginosa, in agreement with lacZ transcriptional fusion data. Still, the Q3Am mutation yielded modest expression of mexXY, less than that seen for antibiotic-treated wild-type P. aeruginosa. These data suggest that PA5471 is not sufficient for MexXY recruitment in response to antibiotic exposure and that additional antibiotic-dependent effects are needed.

Biochemical features and functional implications of the RNA-based T-box regulatory mechanism

The T-box mechanism is a common regulatory strategy used for modulating the expression of genes of amino acid metabolism-related operons in gram-positive bacteria, especially members of the Firmicutes. T-box regulation is usually based on a transcription attenuation mechanism in which an interaction between a specific uncharged tRNA and the 5' region of the transcript stabilizes an antiterminator structure in preference to a terminator structure, thereby preventing transcription termination. Although single T-box regulatory elements are common, double or triple T-box arrangements are also observed, expanding the regulatory range of these elements. In the present study, we predict the functional implications of T-box regulation in genes encoding aminoacyl-tRNA synthetases, proteins of amino acid biosynthetic pathways, transporters, and regulatory proteins. We also consider the global impact of the use of this regulatory mechanism on cell physiology. Novel biochemical relationships between regulated genes and their corresponding metabolic pathways were revealed. Some of the genes identified, such as the quorum-sensing gene luxS, in members of the Lactobacillaceae were not previously predicted to be regulated by the T-box mechanism. Our analyses also predict an imbalance in tRNA sensing during the regulation of operons containing multiple aminoacyl-tRNA synthetase genes or biosynthetic genes involved in pathways common to more than one amino acid. Based on the distribution of T-box regulatory elements, we propose that this regulatory mechanism originated in a common ancestor of members of the Firmicutes, Chloroflexi, Deinococcus-Thermus group, and Actinobacteria and was transferred into the Deltaproteobacteria by horizontal gene transfer.

Small membrane proteins found by comparative genomics and ribosome binding site models

The correct annotation of genes encoding the smallest proteins is one of the biggest challenges of genome annotation, and perhaps more importantly, few annotated short open reading frames have been confirmed to correspond to synthesized proteins. We used sequence conservation and ribosome binding site models to predict genes encoding small proteins, defined as having 16-50 amino acids, in the intergenic regions of the Escherichia coli genome. We tested expression of these predicted as well as previously annotated genes by integrating the sequential peptide affinity tag directly upstream of the stop codon on the chromosome and assaying for synthesis using immunoblot assays. This approach confirmed that 20 previously annotated and 18 newly discovered proteins of 16-50 amino acids are synthesized. We summarize the properties of these small proteins; remarkably more than half of the proteins are predicted to be single-transmembrane proteins, nine of which we show co-fractionate with cell membranes.