RNA stem-loop enhanced expression of previously non-expressible genes

Probing the orthogonality and robustness of the mammalian RNA-binding protein Musashi-1 in Escherichia coli

RNA recognition motifs (RRMs) are widespread RNA-binding protein domains in eukaryotes, which represent promising synthetic biology tools due to their compact structure and efficient activity. Yet, their use in prokaryotes is limited and their functionality poorly characterized. Recently, we repurposed a mammalian Musashi protein containing two RRMs as a translation regulator in Escherichia coli. Here, employing high-throughput RNA sequencing, we explored the impact of Musashi expression on the transcriptomic and translatomic profiles of E. coli, revealing certain metabolic interference, induction of post-transcriptional regulatory processes, and spurious protein-RNA interactions. Engineered Musashi protein mutants displayed compromised regulatory activity, emphasizing the importance of both RRMs for specific and sensitive RNA binding. We found that a mutation known to impede allosteric regulation led to similar translation control activity. Evolutionary experiments disclosed a loss of function of the synthetic circuit in about 40 generations, with the gene coding for the Musashi protein showing a stability comparable to other heterologous genes. Overall, this work expands our understanding of RRMs for post-transcriptional regulation in prokaryotes and highlight their potential for biotechnological and biomedical applications.

Repurposing the mammalian RNA-binding protein Musashi-1 as an allosteric translation repressor in bacteria

The RNA recognition motif (RRM) is the most common RNA-binding protein domain identified in nature. However, RRM-containing proteins are only prevalent in eukaryotic phyla, in which they play central regulatory roles. Here, we engineered an orthogonal post-transcriptional control system of gene expression in the bacterium Escherichia coli with the mammalian RNA-binding protein Musashi-1, which is a stem cell marker with neurodevelopmental role that contains two canonical RRMs. In the circuit, Musashi-1 is regulated transcriptionally and works as an allosteric translation repressor thanks to a specific interaction with the N-terminal coding region of a messenger RNA and its structural plasticity to respond to fatty acids. We fully characterized the genetic system at the population and single-cell levels showing a significant fold change in reporter expression, and the underlying molecular mechanism by assessing the in vitro binding kinetics and in vivo functionality of a series of RNA mutants. The dynamic response of the system was well recapitulated by a bottom-up mathematical model. Moreover, we applied the post-transcriptional mechanism engineered with Musashi-1 to specifically regulate a gene within an operon, implement combinatorial regulation, and reduce protein expression noise. This work illustrates how RRM-based regulation can be adapted to simple organisms, thereby adding a new regulatory layer in prokaryotes for translation control.

A (GCC) repeat in SBF1 reveals a novel biological phenomenon in human and links to late onset neurocognitive disorder

The human SBF1 (SET binding factor 1) gene, alternatively known as MTMR5, is predominantly expressed in the brain, and its epigenetic dysregulation is linked to late-onset neurocognitive disorders (NCDs), such as Alzheimer’s disease. This gene contains a (GCC)-repeat at the interval between + 1 and + 60 of the transcription start site (SBF1-202 ENST00000380817.8). We sequenced the SBF1 (GCC)-repeat in a sample of 542 Iranian individuals, consisting of late-onset NCDs (N = 260) and controls (N = 282). While multiple alleles were detected at this locus, the 8 and 9 repeats were predominantly abundant, forming > 95% of the allele pool across the two groups. Among a number of anomalies, the allele distribution was significantly different in the NCD group versus controls (Fisher’s exact p = 0.006), primarily as a result of enrichment of the 8-repeat in the former. The genotype distribution departed from the Hardy–Weinberg principle in both groups (p < 0.001), and was significantly different between the two groups (Fisher’s exact p = 0.001). We detected significantly low frequency of the 8/9 genotype in both groups, higher frequency of this genotype in the NCD group, and reverse order of 8/8 versus 9/9 genotypes in the NCD group versus controls. Biased heterozygous/heterozygous ratios were also detected for the 6/8 versus 6/9 genotypes (in favor of 6/8) across the human samples studied (Fisher’s exact p = 0.0001). Bioinformatics studies revealed that the number of (GCC)-repeats may change the RNA secondary structure and interaction sites at least across human exon 1. This STR was specifically expanded beyond 2-repeats in primates. In conclusion, we report indication of a novel biological phenomenon, in which there is selection against certain heterozygous genotypes at a STR locus in human. We also report different allele and genotype distribution at this STR locus in late-onset NCD versus controls. In view of the location of this STR in the 5′ untranslated region, RNA/RNA or RNA/DNA heterodimer formation of the involved genotypes and alternative RNA processing and/or translation should be considered.

Guidelines for nucleic acid template design for optimal cell-free protein synthesis using an Escherichia coli reconstituted system or a lysate-based system

Cell-free protein synthesis is an attractive method for generating enzyme/protein variants for simplified functional analysis as both in vitro protein expression and analysis may often be performed in a single vial or well. Today, researchers may choose from multiple commercial cell lysate products or reconstituted systems which are compatible with either mRNA, linear DNA or plasmid DNA templates. Here we provide guidance for optimal design of the genetic elements within linear and plasmid DNA templates which are required to reliably practice cell-free protein synthesis. Protocols are presented for generating linear DNA templates, and data are presented to show that linear DNA templates may in many cases provide robust protein yields even when employing an Escherichia coli lysate for protein synthesis. Finally, the use of linear DNA templates makes it possible to bypass all cell cultivation steps and proceed from PCR amplification of synthetic DNA to generation of target protein in a matter of hours.

On the role of mRNA secondary structure in bacterial translation

Messenger RNA (mRNA) is no longer considered as a mere informational molecule whose sole function is to convey the genetic information specified by DNA to the ribosome. Beyond this primary function, mRNA also contains additional instructions that influence the way and the extent to which this message is translated by the ribosome into protein(s). Indeed, owing to its intrinsic propensity to quickly and dynamically fold and form higher order structures, mRNA exhibits a second layer of structural information specified by the sequence itself. Besides influencing transcription and mRNA stability, this additional information also affects translation, and more precisely the frequency of translation initiation, the choice of open reading frame by recoding, the elongation speed, and the folding of the nascent protein. Many studies in bacteria have shown that mRNA secondary structure participates to the rapid adaptation of these versatile organisms to changing environmental conditions by efficiently tuning translation in response to diverse signals, such as the presence of ligands, regulatory proteins, or small RNAs. This article is categorized under: Regulatory RNAs/RNAi/Riboswitches > Regulatory RNAs RNA Structure and Dynamics > Influence of RNA Structure in Biological Systems Translation > Translation Regulation.

Aliivibrio wodanis as a production host: development of genetic tools for expression of cold-active enzymes

Background

Heterologous production of cold-adapted proteins currently represents one of the greatest bottlenecks in the ongoing bioprospecting efforts to find new enzymes from low-temperature environments, such as, the polar oceans that represent essentially untapped resources in this respect. In mesophilic expression hosts such as Escherichia coli, cold-adapted enzymes often form inactive aggregates. Therefore it is necessary to develop new low-temperature expression systems, including identification of new host organisms and complementary genetic tools. Psychrophilic bacteria, including Pseudoalteromonas haloplanktis, Shewanella and Rhodococcus erythropolis have all been explored as candidates for such applications. However to date none of these have found widespread use as efficient expression systems, or are commercially available. In the present work we explored the use of the sub-Arctic bacterium Aliivibrio wodanis as a potential host for heterologous expression of cold-active enzymes.

Results

We tested 12 bacterial strains, as well as available vectors, promoters and reporter systems. We used RNA-sequencing to determine the most highly expressed genes and their intrinsic promoters in A. wodanis. In addition we examined a novel 5′-fusion to stimulate protein production and solubility. Finally we tested production of a set of “difficult-to-produce” enzymes originating from various bacteria and one Archaea. Our results show that cold-adapted enzymes can be produced in soluble and active form, even in cases when protein production failed in E. coli due to the formation of inclusion bodies. Moreover, we identified a 60-bp/20-aa fragment from the 5′-end of the AW0309160_00174 gene that stimulates expression of Green Fluorescent Protein and improves production of cold-active enzymes when used as a 5′-fusion. A 25-aa peptide from the same protein enhanced secretion of a 25-aa-sfGFP fusion.

Conclusions

Our results indicate the use of A. wodanis and associated genetic tools for low-temperature protein production and indicate that A. wodanis represents an interesting platform for further development of a protein production system that can promote further cold-enzyme discoveries.

Thermosensitive pilus production by FCT type 3 Streptococcus pyogenes controlled by Nra regulator translational efficiency

Streptococcus pyogenes produces a diverse variety of pili in a serotype‐dependent manner and thermosensitive expression of pilus biogenesis genes was previously observed in a serotype M49 strain. However, the precise mechanism and biological significance remain unclear. Herein, the pilus expression analysis revealed the thermosensitive pilus production only in strains possessing the transcriptional regulator Nra. Experimental data obtained for nra deletion and conditional nra‐expressing strains in the background of an M49 strain and the Lactococcus heterologous expression system, indicated that Nra is a positive regulator of pilus genes and also highlighted the importance of the level of intracellular Nra for the thermoregulation of pilus expression. While the nra mRNA level was not significantly influenced by a temperature shift, the Nra protein level was concomitantly increased when the culture temperature was decreased. Intriguingly, a putative stem‐loop structure within the coding region of nra mRNA was a factor related to the post‐transcriptional efficiency of nra mRNA translation. Either deletion of the stem‐loop structure or introduction of silent chromosomal mutations designed to melt the structure attenuated Nra levels, resulting in decreased pilus production. Consequently, the temperature‐dependent translational efficacy of nra mRNA influenced pilus thermoregulation, thereby potentially contributing to the fitness of nra‐positive S. pyogenes in human tissues.

Multiple upstream start codons (AUG) in 5' untranslated region enhance translation efficiency of cry2Ac11 without helper protein

Cry2Ac11, a 65 kDa insecticidal protein produced by Bacillus thuringiensis, shows toxicity against dipteran and lepidopteran larvae. It is encoded by cry2Ac11 gene ( orf3), which is part of an operon comprising orf1, orf2, and orf3. Orf2, a helper protein, helps in proper folding and prevents aberrant aggregation of newly produced molecules. In this study, we have elucidated the effect of different mutations in translation initiation region (TIR), particularly the ribosomal binding site and the start codon (RBS-ATG) on cry2Ac11 gene expression without helper protein. All recombinant constructs were expressed in acrystalliferous B. thuringiensis subsp israelensis 4Q7 under the control of strong chimeric promoter cyt1AP/STAB. Of all the mutants, mut/RBS2, with two consecutive AUGs after the spacer region in TIR, exhibited 89- and 2246-fold higher transcript levels compared with 4Q7-operSalI/RBS ( cry2Ac11 operon) and 4Q7-w-RBS ( cry2Ac11 gene), respectively. The analysis of mut/RBS2 messenger RNA (mRNA) structure in the RBS-AUG region showed the presence of RBS in the single-stranded part of the moderately stable hairpin loop. The high expression efficiency of Cry2Ac11 mutant without helper protein is a cumulative and cooperative result of chimeric promoter cyt1AP/STAB-SD with the optimal context of RBS-AUG region provided by multiple AUGs and stabilizer sequence at 3' ends.

An in Vivo Binding Assay for RNA-Binding Proteins Based on Repression of a Reporter Gene

We study translation repression in bacteria by engineering a regulatory circuit that functions as a binding assay for RNA binding proteins (RBP) in vivo. We do so by inducing expression of a fluorescent protein-RBP chimera, together with encoding its binding site at various positions within the ribosomal initiation region (+11-13 nt from the AUG) of a reporter module. We show that when bound by their cognate RBPs, the phage coat proteins for PP7 (PCP) and Qβ (QCP), strong repression is observed for all hairpin positions within the initiation region. Yet, a sharp transition to no-effect is observed when positioned in the elongation region, at a single-nucleotide resolution. Employing in vivo Selective 2'-hydroxyl acylation analyzed by primer extension followed by sequencing (SHAPE-seq) for a representative construct, established that in the translationally active state the mRNA molecule is nonstructured, while in the repressed state a structured signature was detected. We then utilize this regulatory phenomena to quantify the binding affinity of the coat proteins of phages MS2, PP7, GA, and Qβ to 14 cognate and noncognate binding sites in vivo. Using our circuit, we demonstrate qualitative differences between in vitro to in vivo binding characteristics for various variants when comparing to past studies. Furthermore, by introducing a simple mutation to the loop region for the Qβ-wt site, MCP binding is abolished, creating the first high-affinity QCP site that is completely orthogonal to MCP. Consequently, we demonstrate that our hybrid transcriptional-post-transcriptional circuit can be utilized as a binding assay to quantify RNA-RBP interactions in vivo.

Design and application of genetically-encoded malonyl-CoA biosensors for metabolic engineering of microbial cell factories

Malonyl-CoA is the basic building block for synthesizing a range of important compounds including fatty acids, phenylpropanoids, flavonoids and non-ribosomal polyketides. Centering around malonyl-CoA, we summarized here the various metabolic engineering strategies employed recently to regulate and control malonyl-CoA metabolism and improve cellular productivity. Effective metabolic engineering of microorganisms requires the introduction of heterologous pathways and dynamically rerouting metabolic flux towards products of interest. Transcriptional factor-based biosensors translate an internal cellular signal to a transcriptional output and drive the expression of the designed genetic/biomolecular circuits to compensate the activity loss of the engineered biosystem. Recent development of genetically-encoded malonyl-CoA sensor has stood out as a classical example to dynamically reprogram cell metabolism for various biotechnological applications. Here, we reviewed the design principles of constructing a transcriptional factor-based malonyl-CoA sensor with superior detection limit, high sensitivity and broad dynamic range. We discussed various synthetic biology strategies to remove pathway bottleneck and how genetically-encoded metabolite sensor could be deployed to improve pathway efficiency. Particularly, we emphasized that integration of malonyl-CoA sensing capability with biocatalytic function would be critical to engineer efficient microbial cell factory. Biosensors have also advanced beyond its classical function of a sensor actuator for in situ monitoring of intracellular metabolite concentration. Applications of malonyl-CoA biosensors as a sensor-invertor for negative feedback regulation of metabolic flux, a metabolic switch for oscillatory balancing of malonyl-CoA sink pathway and source pathway and a screening tool for engineering more efficient biocatalyst are also presented in this review. We envision the genetically-encoded malonyl-CoA sensor will be an indispensable tool to optimize cell metabolism and cost-competitively manufacture malonyl-CoA-derived compounds.

The interplay of mRNA stimulatory signals required for AUU-mediated initiation and programmed -1 ribosomal frameshifting in decoding of transposable element IS911

The IS911 bacterial transposable element uses -1 programmed translational frameshifting to generate the protein required for its mobility: translation initiated in one gene (orfA) shifts to the -1 frame and continues in a second overlapping gene (orfB), thus generating the OrfAB transposase. The A-AAA-AAG frameshift site of IS911 is flanked by two stimulatory elements, an upstream Shine-Dalgarno sequence and a downstream stem-loop. We show here that, while they can act independently, these stimulators have a synergistic effect when combined. Mutagenic analyses revealed features of the complex stem-loop that make it a low-efficiency stimulator. They also revealed the dual role of the upstream Shine-Dalgarno sequence as (i) a stimulator of frameshifting, by itself more potent than the stem-loop, and (ii) a mandatory determinant of initiation of OrfB protein synthesis on an AUU codon directly preceding the A6G motif. Both roles rely on transient base pairing of the Shine-Dalgarno sequence with the 3' end of 16S rRNA. Because of its effect on frameshifting, the Shine-Dalgarno sequence is an important determinant of the level of transposase in IS911-containing cells, and hence of the frequency of transposition.

Synonymous but not the same: the causes and consequences of codon bias

Despite their name, synonymous mutations have significant consequences for cellular processes in all taxa. As a result, an understanding of codon bias is central to fields as diverse as molecular evolution and biotechnology. Although recent advances in sequencing and synthetic biology have helped to resolve longstanding questions about codon bias, they have also uncovered striking patterns that suggest new hypotheses about protein synthesis. Ongoing work to quantify the dynamics of initiation and elongation is as important for understanding natural synonymous variation as it is for designing transgenes in applied contexts.

Enhanced expression, detection and purification of recombinant proteins using RNA stem loop and tandem fusion tags

The creation of a double His-tag fusion that forms a RNA stem loop in the mRNA encoding the N-terminus of the target protein is a novel approach for the enhancement of expression, purification, and detection of a recombinant protein. Compared to a single His-tag fusion, a tandem His-tag fusion RNA stem loop, located downstream of the constitutive groE and Ch promoters, enhanced heterologous gene expression in Brucella, Salmonella, and Escherichia. We demonstrated one-step detection and purification of recombinant green fluorescence protein (GFP) directly from Brucella spp. without using Escherichia coli as an expression host. The amount of purified GFP using the tandem His-tag RNA stem loop increased more than threefold; moreover, the sensitivity of detection increased more than fourfold in comparison to the single His-tag fusion form. This method has the potential to significantly improve heterologous gene expression and high-throughput protein synthesis and purification.

RNAs everywhere: genome-wide annotation of structured RNAs

Starting with the discovery of microRNAs and the advent of genome-wide transcriptomics, non-protein-coding transcripts have moved from a fringe topic to a central field research in molecular biology. In this contribution we review the state of the art of "computational RNomics", i.e., the bioinformatics approaches to genome-wide RNA annotation. Instead of rehashing results from recently published surveys in detail, we focus here on the open problem in the field, namely (functional) annotation of the plethora of putative RNAs. A series of exploratory studies are used to provide non-trivial examples for the discussion of some of the difficulties.

mRNA secondary structure at start AUG codon is a key limiting factor for human protein expression in Escherichia coli

Codon usage and thermodynamic optimization of the 5'-end of mRNA have been applied to improve the efficiency of human protein production in Escherichia coli. However, high level expression of human protein in E. coli is still a challenge that virtually depends upon each individual target genes. Using human interleukin 10 (huIL-10) and interferon alpha (huIFN-alpha) coding sequences, we systematically analyzed the influence of several major factors on expression of human protein in E. coli. The results from huIL-10 and reinforced by huIFN-alpha showed that exposing AUG initiator codon from base-paired structure within mRNA itself significantly improved the translation of target protein, which resulted in a 10-fold higher protein expression than the wild-type genes. It was also noted that translation process was not affected by the retained short-range stem-loop structure at Shine-Dalgarno (SD) sequences. On the other hand, codon-optimized constructs of huIL-10 showed unimproved levels of protein expression, on the contrary, led to a remarkable RNA degradation. Our study demonstrates that exposure of AUG initiator codon from long-range intra-strand secondary structure at 5'-end of mRNA may be used as a general strategy for human protein production in E. coli.

Silent mutations at the 5'-end of the cDNA of actinoporins from the sea anemone Stichodactyla helianthus allow their heterologous overproduction in Escherichia coli

Wild-type actinoporins StnI and StnII from the sea anemone Stichodactyla helianthus, as well as their NH(2)-terminal six-His tagged versions, have been overproduced in Escherichia coli. Overproduction of both wild-type proteins was only possible after introducing silent mutations within the 5'-end of their original cDNA sequences. These mutations would prevent the formation of RNA secondary structures blocking the ribosome-binding site and the initiation codon. The four recombinant proteins were purified to homogeneity in milligrams amount and characterized from spectroscopic and functional points of view. All the isolated proteins behaved as the corresponding natural ones although the six-His tagged variants exhibited a decreased lytic activity. The strategy described will be useful to allow the production of mutant variants of these proteins and probably of other actinoporins.

Partition function and base pairing probabilities of RNA heterodimers

BACKGROUND

RNA has been recognized as a key player in cellular regulation in recent years. In many cases, non-coding RNAs exert their function by binding to other nucleic acids, as in the case of microRNAs and snoRNAs. The specificity of these interactions derives from the stability of inter-molecular base pairing. The accurate computational treatment of RNA-RNA binding therefore lies at the heart of target prediction algorithms.

METHODS

The standard dynamic programming algorithms for computing secondary structures of linear single-stranded RNA molecules are extended to the co-folding of two interacting RNAs.

RESULTS

We present a program, RNAcofold, that computes the hybridization energy and base pairing pattern of a pair of interacting RNA molecules. In contrast to earlier approaches, complex internal structures in both RNAs are fully taken into account. RNAcofold supports the calculation of the minimum energy structure and of a complete set of suboptimal structures in an energy band above the ground state. Furthermore, it provides an extension of McCaskill's partition function algorithm to compute base pairing probabilities, realistic interaction energies, and equilibrium concentrations of duplex structures.

Regulation of translation via mRNA structure in prokaryotes and eukaryotes

The mechanism of initiation of translation differs between prokaryotes and eukaryotes, and the strategies used for regulation differ accordingly. Translation in prokaryotes is usually regulated by blocking access to the initiation site. This is accomplished via base-paired structures (within the mRNA itself, or between the mRNA and a small trans-acting RNA) or via mRNA-binding proteins. Classic examples of each mechanism are described. The polycistronic structure of mRNAs is an important aspect of translational control in prokaryotes, but polycistronic mRNAs are not usable (and usually not produced) in eukaryotes. Four structural elements in eukaryotic mRNAs are important for regulating translation: (i) the m7G cap; (ii) sequences flanking the AUG start codon; (iii) the position of the AUG codon relative to the 5' end of the mRNA; and (iv) secondary structure within the mRNA leader sequence. The scanning model provides a framework for understanding these effects. The scanning mechanism also explains how small open reading frames near the 5' end of the mRNA can down-regulate translation. This constraint is sometimes abrogated by changing the structure of the mRNA, sometimes with clinical consequences. Examples are described. Some mistaken ideas about regulation of translation that have found their way into textbooks are pointed out and corrected.

The effect of RNA secondary structures on RNA-ligand binding and the modifier RNA mechanism: a quantitative model

RNA-ligand binding often depends crucially on the local RNA secondary structure at the binding site. We develop here a model that quantitatively predicts the effect of RNA secondary structure on effective RNA-ligand binding activities based on equilibrium thermodynamics and the explicit computations of partition functions for the RNA structures. A statistical test for the impact of a particular structural feature on the binding affinities follows directly from this approach. The formalism is extended to describing the effects of hybridizing small "modifier RNAs" to a target RNA molecule outside its ligand binding site. We illustrate the applicability of our approach by quantitatively describing the interaction of the mRNA stabilizing protein HuR with AU-rich elements. We discuss our model and recent experimental findings demonstrating the effectivity of modifier RNAs in vitro in the context of the current research activities in the field of non-coding RNAs. We speculate that modifier RNAs might also exist in nature; if so, they present an additional regulatory layer for fine-tuning gene expression that could evolve rapidly, leaving no obvious traces in the genomic DNA sequences.