Translational initiation on structured messengers. Another role for the Shine-Dalgarno interaction

Machine Learning Guided Batched Design of a Bacterial Ribosome Binding Site

Optimization of gene expression levels is an essential part of the organism design process. Fine control of this process can be achieved by engineering transcription and translation control elements, including the ribosome binding site (RBS). Unfortunately, the design of specific genetic parts remains challenging because of the lack of reliable design methods. To address this problem, we have created a machine learning guided Design-Build-Test-Learn (DBTL) cycle for the experimental design of bacterial RBSs to demonstrate how small genetic parts can be reliably designed using relatively small, high-quality data sets. We used Gaussian Process Regression for the Learn phase of the cycle and the Upper Confidence Bound multiarmed bandit algorithm for the Design of genetic variants to be tested in vivo. We have integrated these machine learning algorithms with laboratory automation and high-throughput processes for reliable data generation. Notably, by Testing a total of 450 RBS variants in four DBTL cycles, we have experimentally validated RBSs with high translation initiation rates equaling or exceeding our benchmark RBS by up to 34%. Overall, our results show that machine learning is a powerful tool for designing RBSs, and they pave the way toward more complicated genetic devices.

Comparative Analysis of anti-Shine- Dalgarno Function in Flavobacterium johnsoniae and Escherichia coli

The anti-Shine-Dalgarno (ASD) sequence of 16S rRNA is highly conserved across Bacteria, and yet usage of Shine-Dalgarno (SD) sequences in mRNA varies dramatically, depending on the lineage. Here, we compared the effects of ASD mutagenesis in Escherichia coli, a Gammaproteobacteria which commonly employs SD sequences, and Flavobacterium johnsoniae, a Bacteroidia which rarely does. In E. coli, 30S subunits carrying any single substitution at positions 1,535–1,539 confer dominant negative phenotypes, whereas subunits with mutations at positions 1,540–1,542 are sufficient to support cell growth. These data suggest that CCUCC (1,535–1,539) represents the functional core of the element in E. coli. In F. johnsoniae, deletion of three ribosomal RNA (rrn) operons slowed growth substantially, a phenotype largely rescued by a plasmid-borne copy of the rrn operon. Using this complementation system, we found that subunits with single mutations at positions 1,535–1,537 are as active as control subunits, in sharp contrast to the E. coli results. Moreover, subunits with quadruple substitution or complete replacement of the ASD retain substantial, albeit reduced, activity. Sedimentation analysis revealed that these mutant subunits are overrepresented in the subunit fractions and underrepresented in polysome fractions, suggesting some defect in 30S biogenesis and/or translation initiation. Nonetheless, our collective data indicate that the ASD plays a much smaller role in F. johnsoniae than in E. coli, consistent with SD usage in the two organisms.

Importance of the 5' regulatory region to bacterial synthetic biology applications

The field of synthetic biology is evolving at a fast pace. It is advancing beyond single-gene alterations in single hosts to the logical design of complex circuits and the development of integrated synthetic genomes. Recent breakthroughs in deep learning, which is increasingly used in de novo assembly of DNA components with predictable effects, are also aiding the discipline. Despite advances in computing, the field is still reliant on the availability of pre-characterized DNA parts, whether natural or synthetic, to regulate gene expression in bacteria and make valuable compounds. In this review, we discuss the different bacterial synthetic biology methodologies employed in the creation of 5' regulatory regions - promoters, untranslated regions and 5'-end of coding sequences. We summarize methodologies and discuss their significance for each of the functional DNA components, and highlight the key advances made in bacterial engineering by concentrating on their flaws and strengths. We end the review by outlining the issues that the discipline may face in the near future.

Expressing Cloned Genes for Protein Production, Purification, and Analysis

Obtaining high quantities of a specific protein directly from native sources is often challenging, particularly when dealing with human proteins. To overcome this obstacle, many researchers take advantage of heterologous expression systems by cloning genes into artificial vectors designed to operate within easily cultured cells, such as Escherichia coli, Pichia pastoris (yeast), and several varieties of insect and mammalian cells. Heterologous expression systems also allow for easy modification of the protein to optimize expression, mutational analysis of specific sites within the protein and facilitate their purification with engineered affinity tags. Some degree of purification of the target protein is usually required for functional analysis. Purification to near homogeneity is essential for characterization of protein structure by X-ray crystallography or nuclear magnetic resonance (NMR) and characterization of the biochemical and biophysical properties of a protein, because contaminating proteins almost always adversely affect the results. Methods for producing and purifying proteins in several different expression platforms and using a variety of vectors are introduced here.

Structural basis of sequestration of the anti-Shine-Dalgarno sequence in the Bacteroidetes ribosome

Genomic studies have indicated that certain bacterial lineages such as the Bacteroidetes lack Shine-Dalgarno (SD) sequences, and yet with few exceptions ribosomes of these organisms carry the canonical anti-SD (ASD) sequence. Here, we show that ribosomes purified from Flavobacterium johnsoniae, a representative of the Bacteroidetes, fail to recognize the SD sequence of mRNA in vitro. A cryo-electron microscopy structure of the complete 70S ribosome from F. johnsoniae at 2.8 Å resolution reveals that the ASD is sequestered by ribosomal proteins bS21, bS18 and bS6, explaining the basis of ASD inhibition. The structure also uncovers a novel ribosomal protein—bL38. Remarkably, in F. johnsoniae and many other Flavobacteriia, the gene encoding bS21 contains a strong SD, unlike virtually all other genes. A subset of Flavobacteriia have an alternative ASD, and in these organisms the fully complementary sequence lies upstream of the bS21 gene, indicative of natural covariation. In other Bacteroidetes classes, strong SDs are frequently found upstream of the genes for bS21 and/or bS18. We propose that these SDs are used as regulatory elements, enabling bS21 and bS18 to translationally control their own production.

Improving heterologous protein expression in Synechocystis sp. PCC 6803 for alpha-bisabolene production

Cyanobacterial biofuels have the potential to reduce the cost and climate impacts of biofuel production because primary carbon fixation and conversion to fuel are completed together in the cultivation of the cyanobacteria. Cyanobacterial biofuels, therefore, do not rely on costly organic carbon feedstocks that heterotrophs require, which reduces competition for agricultural resources such as arable land and freshwater. However, the published product titer achieved for most molecules of interest using cyanobacteria lag behind what has been achieved using yeast and Escherichia coli (E. coli) cultures. In Synechocystis sp. PCC 6803 (S. 6803), we attempted to increase the product titer of the sesquiterpene, bisabolene, which may be converted to bisabolane, a possible diesel replacement. We tested 19 strains of genetically modified S. 6803 with five different codon usage sequences of the bisabolene synthase from the grand fir tree (Abies grandis). At least three ribosome binding sites (most designed using the RBS Calculator) were tested for each codon usage sequence. We also tested strains with and without the farnesyl pyrophosphate synthase gene from E. coli. Bisabolene titers after five days of growth in continuous light ranged from un-detected to 7.8 ​mg/L. Bisabolene synthase abundance was measured and found to be well correlated with titer. Select strains were also tested in 12:12 light:dark cycles, where similar titers were reached after the same amount of light exposure time. One engineered strain was also tested in photobioreactors exposed to a simulated outdoor light pattern with maximum light intensity of 1600 ​μmol photons m-2 s-1. Here, the bisabolene titer reached 22.2 ​mg/L after 36 days of growth. Dramatic improvements in our ability to control gene expression in cyanobacteria such as S. 6803, and the co-utilization of additional metabolic engineering methods, are needed in order for these titers to improve to the levels reported for engineered E. coli.

Protein unties the pseudoknot: S1-mediated unfolding of RNA higher order structure

Ribosomal protein S1 plays important roles in the translation initiation step of many Escherichia coli mRNAs, particularly those with weak Shine-Dalgarno sequences or structured 5′ UTRs, in addition to a variety of cellular processes beyond the ribosome. In all cases, the RNA-binding activity of S1 is a central feature of its function. While sequence determinants of S1 affinity and many elements of the interactions of S1 with simple secondary structures are known, mechanistic details of the protein's interactions with RNAs of more complex secondary and tertiary structure are less understood. Here, we investigate the interaction of S1 with the well-characterized H-type pseudoknot of a class-I translational preQ₁ riboswitch as a highly structured RNA model whose conformation and structural dynamics can be tuned by the addition of ligands of varying binding affinity, particularly preQ₁, guanine, and 2,6-diaminopurine. Combining biochemical and single molecule fluorescence approaches, we show that S1 preferentially interacts with the less folded form of the pseudoknot and promotes a dynamic, partially unfolded conformation. The ability of S1 to unfold the RNA is inversely correlated with the structural stability of the pseudoknot. These mechanistic insights delineate the scope and limitations of S1-chaperoned unfolding of structured RNAs.

Unstructured 5'-tails act through ribosome standby to override inhibitory structure at ribosome binding sites

Initiation is the rate-limiting step in translation. It is well-known that stable structure at a ribosome binding site (RBS) impedes initiation. The ribosome standby model of de Smit and van Duin, based on studies of the MS2 phage coat cistron, proposed how high translation rates can be reconciled with stable, inhibitory structures at an RBS. Here, we revisited the coat protein system and assessed the translation efficiency from its sequestered RBS by introducing standby mutations. Further experiments with gfp reporter constructs assessed the effects of 5′-tails—as standby sites—with respect to length and sequence contributions. In particular, combining in vivo and in vitro assays, we can show that tails of CA-dinucleotide repeats—and to a lesser extent, AU-repeats—dramatically increase translation rates. Tails of increasing length reach maximal rate-enhancing effects at 16–18 nucleotides. These standby tails are single-stranded and do not exert their effect by structure changes in the neighboring RBS stem–loop. In vitro translation and toeprinting assays furthermore demonstrate that standby effects are exerted at the level of translation initiation. Finally, as expected, destabilizing mutations within the coat RBS indicate an interplay with the effects of standby tails.

Conformational switch in the ribosomal protein S1 guides unfolding of structured RNAs for translation initiation

Initiation of bacterial translation requires that the ribosome-binding site in mRNAs adopts single-stranded conformations. In Gram-negative bacteria the ribosomal protein S1 (rS1) is a key player in resolving of structured elements in mRNAs. However, the exact mechanism of how rS1 unfolds persistent secondary structures in the translation initiation region (TIR) is still unknown. Here, we show by NMR spectroscopy that Vibrio vulnificus rS1 displays a unique architecture of its mRNA-binding domains, where domains D3 and D4 provide the mRNA-binding platform and cover the nucleotide binding length of the full-length rS1. D5 significantly increases rS1’s chaperone activity, although it displays structural heterogeneity both in isolation and in presence of the other domains, albeit to varying degrees. The heterogeneity is induced by the switch between the two equilibrium conformations and is triggered by an order-to-order transition of two mutually exclusive secondary structures (β-strand-to-α-helix) of the ‘AERERI’ sequence. The conformational switching is exploited for melting of structured 5′-UTR’s, as the conformational heterogeneity of D5 can compensate the entropic penalty of complex formation. Our data thus provides a detailed understanding of the intricate coupling of protein and RNA folding dynamics enabling translation initiation of structured mRNAs.

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.

Prokaryotic coding regions have little if any specific depletion of Shine-Dalgarno motifs

The Shine-Dalgarno motif occurs in front of prokaryotic start codons, and is complementary to the 3’ end of the 16S ribosomal RNA. Hybridization between the Shine-Dalgarno sequence and the anti-Shine-Dalgarno region of the16S rRNA (CCUCCU) directs the ribosome to the start AUG of the mRNA for translation. Shine-Dalgarno-like motifs (AGGAGG in E. coli) are depleted from open reading frames of most prokaryotes. This may be because hybridization of the 16S rRNA at Shine-Dalgarnos inside genes would slow translation or induce internal initiation. However, we analyzed 128 species from diverse phyla where the 16S rRNA gene(s) lack the anti-Shine-Dalgarno sequence, and so the 16S rRNA is incapable of interacting with Shine-Dalgarno-like sequences. Despite this lack of an anti-Shine-Dalgarno, half of these species still displayed depletion of Shine-Dalgarno-like sequences when analyzed by previous methods. Depletion of the same G-rich sequences was seen by these methods even in eukaryotes, which do not use the Shine-Dalgarno mechanism. We suggest previous methods are partly detecting a non-specific depletion of G-rich sequences. Alternative informatics approaches show that most prokaryotes have only slight, if any, specific depletion of Shine-Dalgarno-like sequences from open reading frames. Together with recent evidence that ribosomes do not pause at ORF-internal Shine-Dalgarno motifs, these results suggest the presence of ORF-internal Shine-Dalgarno-like motifs may be inconsequential, perhaps because internal regions of prokaryotic mRNAs may be structurally “shielded” from translation initiation.

Re-annotation of 12,495 prokaryotic 16S rRNA 3' ends and analysis of Shine-Dalgarno and anti-Shine-Dalgarno sequences

We examined 20,648 prokaryotic unique taxids with respect to the annotation of the 3' end of the 16S rRNA, which contains the anti-Shine-Dalgarno sequence. We used the sequence of highly conserved helix 45 of the 16S rRNA as a guide. By this criterion, 8,153 annotated 3' ends correctly included the anti-Shine-Dalgarno sequence, but 12,495 were foreshortened or otherwise mis-annotated, missing part or all of the anti-Shine-Dalgarno sequence, which immediately follows helix 45. We re-annotated, giving a total of 20,648 16S rRNA 3' ends. The vast majority indeed contained a consensus anti-Shine-Dalgarno sequence, embedded in a highly conserved 13 base "tail". However, 128 exceptional organisms had either a variant anti-Shine-Dalgarno, or no recognizable anti-Shine-Dalgarno, in their 16S rRNA(s). For organisms both with and without an anti-Shine-Dalgarno, we identified the Shine-Dalgarno motifs actually enriched in front of each organism's open reading frames. This showed to what extent the Shine-Dalgarno motifs correlated with anti-Shine Dalgarno motifs. In general, organisms whose rRNAs lacked a perfect anti-Shine-Dalgarno motif also lacked a recognizable Shine-Dalgarno. For organisms whose 16S rRNAs contained a perfect anti-Shine-Dalgarno motif, a variety of results were obtained. We found one genus, Alteromonas, where several taxids apparently maintain two different types of 16S rRNA genes, with different, but conserved, antiSDs. The fact that some organisms do not seem to have or use Shine-Dalgarno motifs supports the idea that prokaryotes have other robust mechanisms for recognizing start codons for translation.

Genome scale analysis of Escherichia coli with a comprehensive prokaryotic sequence-based biophysical model of translation initiation and elongation

Translation initiation in prokaryotes is affected by the mRNA folding and interaction of the ribosome binding site with the ribosomal RNA. The elongation rate is affected, among other factors, by the local biophysical properties of the coding regions, the decoding rates of different codons, and the interactions among ribosomes. Currently, there is no comprehensive biophysical model of translation that enables the prediction of mRNA translation dynamics based only on the transcript sequence and while considering all of these fundamental aspects of translation. In this study, we provide, for the first time, a computational simulative biophysical model of both translation initiation and elongation with all aspects mentioned above. We demonstrate our model performance and advantages focusing on Escherichia coli genes. We further show that the model enables prediction of translation rate, protein levels, and ribosome densities. In addition, our model enables quantifying the effect of silent mutations on translation rate in different parts of the transcript, the relative effect of mutations on translation initiation and elongation, and the effect of mutations on ribosome traffic jams. Thus, unlike previous models, the proposed one provides comprehensive information, facilitating future research in disciplines such as molecular evolution, synthetic biology, and functional genomics. A toolkit to estimate translation dynamics of transcripts is available at: https://www.cs.tau.ac.il/∼tamirtul/transim.

Design of a Temperature-Responsive Transcription Terminator

RNA structures regulate various steps in gene expression. Transcription in bacteria is typically terminated by stable hairpin structures. Translation initiation can be modulated by metabolite- or temperature-sensitive RNA structures, called riboswitches or RNA thermometers (RNATs), respectively. RNATs control translation initiation by occlusion of the ribosome binding site at low temperatures. Increasing temperatures destabilize the RNA structure and facilitate ribosome access. In this study, we exploited temperature-responsive RNAT structures to design regulatory elements that control transcription termination instead of translation initiation in Escherichia coli. In order to mimic the structure of factor-independent intrinsic terminators, naturally occurring RNAT hairpins were genetically engineered to be followed by a U-stretch. Functional temperature-responsive terminators (thermoterms) prevented mRNA synthesis at low temperatures but resumed transcription after a temperature upshift. The successful design of temperature-controlled terminators highlights the potential of RNA structures as versatile gene expression control elements.

Characterizing the Structure-Function Relationship of a Naturally Occurring RNA Thermometer

A large number of bacteria have been found to govern virulence and heat shock responses using temperature-sensing RNAs known as RNA thermometers. A prime example is the agsA thermometer known to regulate the production of the AgsA heat shock protein in Salmonella enterica using a "fourU" structural motif. Using the SHAPE-Seq RNA structure-probing method in vivo and in vitro, we found that the regulator functions by a subtle shift in equilibrium RNA structure populations that leads to a partial melting of the helix containing the ribosome binding site. We also demonstrate that binding of the ribosome to the agsA mRNA causes changes to the thermometer structure that appear to facilitate thermometer helix unwinding. These results demonstrate how subtle RNA structural changes can govern gene expression and illuminate the function of an important bacterial regulatory motif.

Elucidating the 16S rRNA 3' boundaries and defining optimal SD/aSD pairing in Escherichia coli and Bacillus subtilis using RNA-Seq data

Bacterial translation initiation is influenced by base pairing between the Shine-Dalgarno (SD) sequence in the 5' UTR of mRNA and the anti-SD (aSD) sequence at the free 3' end of the 16S rRNA (3' TAIL) due to: 1) the SD/aSD sequence binding location and 2) SD/aSD binding affinity. In order to understand what makes an SD/aSD interaction optimal, we must define: 1) terminus of the 3' TAIL and 2) extent of the core aSD sequence within the 3' TAIL. Our approach to characterize these components in Escherichia coli and Bacillus subtilis involves 1) mapping the 3' boundary of the mature 16S rRNA using high-throughput RNA sequencing (RNA-Seq), and 2) identifying the segment within the 3' TAIL that is strongly preferred in SD/aSD pairing. Using RNA-Seq data, we resolve previous discrepancies in the reported 3' TAIL in B. subtilis and recovered the established 3' TAIL in E. coli. Furthermore, we extend previous studies to suggest that both highly and lowly expressed genes favor SD sequences with intermediate binding affinity, but this trend is exclusive to SD sequences that complement the core aSD sequences defined herein.

Diversity of Translation Initiation Mechanisms across Bacterial Species Is Driven by Environmental Conditions and Growth Demands

The Shine-Dalgarno (SD) sequence motif is frequently found upstream of protein coding genes and is thought to be the dominant mechanism of translation initiation used by bacteria. Experimental studies have shown that the SD sequence facilitates start codon recognition and enhances translation initiation by directly interacting with the highly conserved anti-SD sequence on the 30S ribosomal subunit. However, the proportion of SD-led genes within a genome varies across species and the factors governing this variation in translation initiation mechanisms remain largely unknown. Here, we conduct a phylogenetically informed analysis and find that species capable of rapid growth contain a higher proportion of SD-led genes throughout their genomes. We show that SD sequence utilization covaries with a suite of genomic features that are important for efficient translation initiation and elongation. In addition to these endogenous genomic factors, we further show that exogenous environmental factors may influence the evolution of translation initiation mechanisms by finding that thermophilic species contain significantly more SD-led genes than mesophiles. Our results demonstrate that variation in translation initiation mechanisms across bacterial species is predictable and is a consequence of differential life-history strategies related to maximum growth rate and environmental-specific constraints.

Leveraging genome-wide datasets to quantify the functional role of the anti-Shine-Dalgarno sequence in regulating translation efficiency

Studies dating back to the 1970s established that sequence complementarity between the anti-Shine–Dalgarno (aSD) sequence on prokaryotic ribosomes and the 5′ untranslated region of mRNAs helps to facilitate translation initiation. The optimal location of aSD sequence binding relative to the start codon, the full extents of the aSD sequence and the functional form of the relationship between aSD sequence complementarity and translation efficiency have not been fully resolved. Here, we investigate these relationships by leveraging the sequence diversity of endogenous genes and recently available genome-wide estimates of translation efficiency. We show that—after accounting for predicted mRNA structure—aSD sequence complementarity increases the translation of endogenous mRNAs by roughly 50%. Further, we observe that this relationship is nonlinear, with translation efficiency maximized for mRNAs with intermediate levels of aSD sequence complementarity. The mechanistic insights that we observe are highly robust: we find nearly identical results in multiple datasets spanning three distantly related bacteria. Further, we verify our main conclusions by re-analysing a controlled experimental dataset.

Precise quantification of translation inhibition by mRNA structures that overlap with the ribosomal footprint in N-terminal coding sequences

A mRNA's translation rate is controlled by several sequence determinants, including the presence of RNA structures within the N-terminal regions of its coding sequences. However, the physical rules that govern when such mRNA structures will inhibit translation remain unclear. Here, we introduced systematically designed RNA hairpins into the N-terminal coding region of a reporter protein with steadily increasing distances from the start codon, followed by characterization of their mRNA and expression levels in Escherichia coli. We found that the mRNAs' translation rates were repressed, by up to 530-fold, when mRNA structures overlapped with the ribosome's footprint. In contrast, when the mRNA structure was located outside the ribosome's footprint, translation was repressed by <2-fold. By combining our measurements with biophysical modeling, we determined that the ribosomal footprint extends 13 nucleotides into the N-terminal coding region and, when a mRNA structure overlaps or partially overlaps with the ribosomal footprint, the free energy to unfold only the overlapping structure controlled the extent of translation repression. Overall, our results provide precise quantification of the rules governing translation initiation at N-terminal coding regions, improving the predictive design of post-transcriptional regulatory elements that regulate translation rate.

DAMBE6: New Tools for Microbial Genomics, Phylogenetics, and Molecular Evolution

DAMBE is a comprehensive software workbench for data analysis in molecular biology, phylogenetics, and evolution. Several important new functions have been added since version 5 of DAMBE: 1) comprehensive genomic profiling of translation initiation efficiency of different genes in different prokaryotic species, 2) a new index of translation elongation (I_TE) that takes into account both tRNA-mediated selection and background mutation on codon–anticodon adaptation, 3) a new and accurate phylogenetic approach based on pairwise alignment only, which is useful for highly divergent sequences from which a reliable multiple sequence alignment is difficult to obtain. Many other functions have been updated and improved including PWM for motif characterization, Gibbs sampler for de novo motif discovery, hidden Markov models for protein secondary structure prediction, self-organizing map for nonlinear clustering of transcriptomic data, comprehensive sequence alignment, and phylogenetic functions. DAMBE features a graphic, user-friendly and intuitive interface, and is freely available from http://dambe.bio.uottawa.ca.

Novel Translation Initiation Regulation Mechanism in Escherichia coli ptrB Mediated by a 5'-Terminal AUG

Alternative translation initiation mechanisms, distinct from the Shine-Dalgarno (SD) sequence-dependent mechanism, are more prevalent in bacteria than once anticipated. Translation of Escherichia coliptrB instead requires an AUG triplet at the 5' terminus of its mRNA. The 5'-terminal AUG (5'-uAUG) acts as a ribosomal recognition signal to attract ribosomes to the ptrB mRNA rather than functioning as an initiation codon to support translation of an upstream open reading frame. ptrB expression exhibits a stronger dependence on the 5'-uAUG than the predicted SD sequence; however, strengthening the predicted ptrB SD sequence relieves the necessity for the 5'-uAUG. Additional sequences within the ptrB 5' untranslated region (5'-UTR) work cumulatively with the 5'-uAUG to control expression of the downstream ptrB coding sequence (CDS), thereby compensating for the weak SD sequence. Replacement of 5'-UTRs from other mRNAs with the ptrB 5'-UTR sequence showed a similar dependence on the 5'-uAUG for CDS expression, suggesting that the regulatory features contained within the ptrB 5'-UTR are sufficient to control the expression of other E. coli CDSs. Demonstration that the 5'-uAUG present on the ptrB leader mRNA is involved in ribosome binding and expression of the downstream ptrB CDS revealed a novel form of translational regulation. Due to the abundance of AUG triplets at the 5' termini of E. coli mRNAs and the ability of ptrB 5'-UTR regulation to function independently of gene context, the regulatory effects of 5'-uAUGs on downstream CDSs may be widespread throughout the E. coli genome.IMPORTANCE As the field of synthetic biology continues to grow, a complete understanding of basic biological principles will be necessary. The increasing complexity of the synthetic systems highlights the gaps in our current knowledge of RNA regulation. This study demonstrates that there are novel ways to regulate canonical Shine-Dalgarno-led mRNAs in Escherichia coli, illustrating that our understanding of the fundamental processes of translation and RNA regulation is still incomplete. Even for E. coli, one of the most-studied model organisms, genes with translation initiation mechanisms that do not fit the canonical Shine-Dalgarno sequence paradigm are being revealed. Uncovering diverse mechanisms that control translational expression will allow synthetic biologists to finely tune protein production of desired gene products.

Translation Initiation is Controlled by RNA Folding Kinetics via a Ribosome Drafting Mechanism

RNA folding plays an important role in controlling protein synthesis as well as other cellular processes. Existing models have focused on how RNA folding energetics control translation initiation rate under equilibrium conditions but have largely ignored the effects of nonequilibrium RNA folding. We introduce a new mechanism, called "ribosome drafting", that explains how a mRNA's folding kinetics and the ribosome's binding rate collectively control its translation initiation rate. During cycles of translation, ribosome drafting emerges whenever successive ribosomes bind to a mRNA faster than the mRNA can refold, maintaining it in a nonequilibrium state with an acceleration of protein synthesis. Using computational design, time-correlated single photon counting, and expression measurements, we demonstrate that slow-folding and fast-folding RNA structures with equivalent folding energetics can vary protein synthesis rates by 1000-fold. We determine the necessary conditions for ribosome drafting by characterizing mRNAs with rationally designed ribosome binding rates, folding kinetics, and folding energetics, confirming the predictions of a nonequilibrium Markov model of translation. Our results have widespread implications, illustrating how competitive folding and assembly kinetics can shape the gene expression machinery's sequence-structure-function relationship inside cells.

The MazF-regulon: a toolbox for the post-transcriptional stress response in Escherichia coli

Flexible adaptation to environmental stress is vital for bacteria. An energy-efficient post-transcriptional stress response mechanism in Escherichia coli is governed by the toxin MazF. After stress-induced activation the endoribonuclease MazF processes a distinct subset of transcripts as well as the 16S ribosomal RNA in the context of mature ribosomes. As these 'stress-ribosomes' are specific for the MazF-processed mRNAs, the translational program is changed. To identify this 'MazF-regulon' we employed Poly-seq (polysome fractionation coupled with RNA-seq analysis) and analyzed alterations introduced into the transcriptome and translatome after mazF overexpression. Unexpectedly, our results reveal that the corresponding protein products are involved in all cellular processes and do not particularly contribute to the general stress response. Moreover, our findings suggest that translational reprogramming serves as a fast-track reaction to harsh stress and highlight the so far underestimated significance of selective translation as a global regulatory mechanism in gene expression. Considering the reported implication of toxin-antitoxin (TA) systems in persistence, our results indicate that MazF acts as a prime effector during harsh stress that potentially introduces translational heterogeneity within a bacterial population thereby stimulating persister cell formation.

A universal strategy for regulating mRNA translation in prokaryotic and eukaryotic cells

We describe a simple strategy to control mRNA translation in both prokaryotic and eukaryotic cells which relies on a unique protein–RNA interaction. Specifically, we used the Pumilio/FBF (PUF) protein to repress translation by binding in between the ribosome binding site (RBS) and the start codon (in Escherichia coli), or by binding to the 5′ untranslated region of target mRNAs (in mammalian cells). The design principle is straightforward, the extent of translational repression can be tuned and the regulator is genetically encoded, enabling the construction of artificial signal cascades. We demonstrate that this approach can also be used to regulate polycistronic mRNAs; such regulation has rarely been achieved in previous reports. Since the regulator used in this study is a modular RNA-binding protein, which can be engineered to target different 8-nucleotide RNA sequences, our strategy could be used in the future to target endogenous mRNAs for regulating metabolic flows and signaling pathways in both prokaryotic and eukaryotic cells.

Escherichia coli and Staphylococcus phages: effect of translation initiation efficiency on differential codon adaptation mediated by virulent and temperate lifestyles

Rapid biosynthesis is key to the success of bacteria and viruses. Highly expressed genes in bacteria exhibit a strong codon bias corresponding to the differential availability of tRNAs. However, a large clade of lambdoid coliphages exhibits relatively poor codon adaptation to the host translation machinery, in contrast to other coliphages that exhibit strong codon adaptation to the host. Three possible explanations were previously proposed but dismissed: (1) the phage-borne tRNA genes that reduce the dependence of phage translation on host tRNAs, (2) lack of time needed for evolving codon adaptation due to recent host switching, and (3) strong strand asymmetry with biased mutation disrupting codon adaptation. Here, we examined the possibility that phages with relatively poor codon adaptation have poor translation initiation which would weaken the selection on codon adaptation. We measured translation initiation by: (1) the strength and position of the Shine-Dalgarno (SD) sequence, and (2) the stability of the secondary structure of sequences flanking the SD and start codon known to affect accessibility of the SD sequence and start codon. Phage genes with strong codon adaptation had significantly stronger SD sequences than those with poor codon adaptation. The former also had significantly weaker secondary structure in sequences flanking the SD sequence and start codon than the latter. Thus, lambdoid phages do not exhibit strong codon adaptation because they have relatively inefficient translation initiation and would benefit little from increased elongation efficiency. We also provided evidence suggesting that phage lifestyle (virulent versus temperate) affected selection intensity on the efficiency of translation initiation and elongation.

Translation rate is controlled by coupled trade-offs between site accessibility, selective RNA unfolding and sliding at upstream standby sites

The ribosome's interactions with mRNA govern its translation rate and the effects of post-transcriptional regulation. Long, structured 5' untranslated regions (5' UTRs) are commonly found in bacterial mRNAs, though the physical mechanisms that determine how the ribosome binds these upstream regions remain poorly defined. Here, we systematically investigate the ribosome's interactions with structured standby sites, upstream of Shine-Dalgarno sequences, and show that these interactions can modulate translation initiation rates by over 100-fold. We find that an mRNA's translation initiation rate is controlled by the amount of single-stranded surface area, the partial unfolding of RNA structures to minimize the ribosome's binding free energy penalty, the absence of cooperative binding and the potential for ribosomal sliding. We develop a biophysical model employing thermodynamic first principles and a four-parameter free energy model to accurately predict the ribosome's translation initiation rates for 136 synthetic 5' UTRs with large structures, diverse shapes and multiple standby site modules. The model predicts and experiments confirm that the ribosome can readily bind distant standby site modules that support high translation rates, providing a physical mechanism for observed context effects and long-range post-transcriptional regulation.

Translation enhancer improves the ribosome liberation from translation initiation

For translation initiation in bacteria, the Shine-Dalgarno (SD) and anti-SD sequence of the 30S subunit play key roles for specific interactions between ribosomes and mRNAs to determine the exact position of the translation initiation region. However, ribosomes also must dissociate from the translation initiation region to slide toward the downstream sequence during mRNA translation. Translation enhancers upstream of the SD sequences of mRNAs, which likely contribute to a direct interaction with ribosome protein S1, enhance the yields of protein biosynthesis. Nevertheless, the mechanism of the effect of translation enhancers to initiate the translation is still unknown. In this paper, we investigated the effects of the SD and enhancer sequences on the binding kinetics of the 30S ribosomal subunits to mRNAs and their translation efficiencies. mRNAs with both the SD and translation enhancers promoted the amount of protein synthesis but destabilized the interaction between the 30S subunit and mRNA by increasing the dissociation rate constant (koff) of the 30S subunit. Based on a model for kinetic parameters, a 16-fold translation efficiency could be achieved by introducing a tandem repeat of adenine sequences (A20) between the SD and translation enhancer sequences. Considering the results of this study, translation enhancers with an SD sequence regulate ribosomal liberation from translation initiation to determine the translation efficiency of the downstream coding region.

Short ROSE-like RNA thermometers control IbpA synthesis in Pseudomonas species

The bacterial small heat shock protein IbpA protects client proteins from aggregation. Due to redundancy in the cellular chaperone network, deletion of the ibpA gene often leads to only a mild or no phenotypic defect. In this study, we show that a Pseudomonas putida ibpA deletion mutant has a severe growth defect under heat stress conditions and reduced survival during recovery revealing a critical role of IbpA in heat tolerance. Transcription of the ibpA gene depends on the alternative heat shock sigma factor σ³². Production of IbpA protein only at heat shock temperatures suggested additional translational control. We conducted a comprehensive structural and functional analysis of the 5′ untranslated regions of the ibpA genes from P. putida and Pseudomonas aeruginosa. Both contain a ROSE-type RNA thermometer that is substantially shorter and simpler than previously reported ROSE elements. Comprised of two hairpin structures only, they inhibit translation at low temperature and permit translation initiation after a temperature upshift. Both elements regulate reporter gene expression in Escherichia coli and ribosome binding in vitro in a temperature-dependent manner. Structure probing revealed local melting of the second hairpin whereas the first hairpin remained unaffected. High sequence and structure conservation of pseudomonal ibpA untranslated regions and their ability to confer thermoregulation in vivo suggest that short ROSE-like thermometers are commonly used to control IbpA synthesis in Pseudomonas species.

A comparative analysis of the properties of regulated promoter systems commonly used for recombinant gene expression in Escherichia coli

Background

Production of recombinant proteins in bacteria for academic and commercial purposes is a well established field; however the outcomes of process developments for specific proteins are still often unpredictable. One reason is the limited understanding of the performance of expression cassettes relative to each other due to different genetic contexts. Here we report the results of a systematic study aiming at exclusively comparing commonly used regulator/promoter systems by standardizing the designs of the replicon backbones.

Results

The vectors used in this study are based on either the RK2- or the pMB1- origin of replication and contain the regulator/promoter regions of XylS/Pm (wild-type), XylS/Pm ML1-17 (a Pm variant), LacI/P_T7lac, LacI/P_trc and AraC/P_BAD to control expression of different proteins with various origins. Generally and not unexpected high expression levels correlate with high replicon copy number and the LacI/P_T7lac system generates more transcript than all the four other cassettes. However, this transcriptional feature does not always lead to a correspondingly more efficient protein production, particularly if protein functionality is considered. In most cases the XylS/Pm ML1-17 and LacI/P_T7lac systems gave rise to the highest amounts of functional protein production, and the XylS/Pm ML1-17 is the most flexible in the sense that it does not require any specific features of the host. The AraC/P_BAD system is very good with respect to tightness, and a commonly used bioinformatics prediction tool (RBS calculator) suggested that it has the most translation-efficient UTR. Expression was also studied by flow cytometry in individual cells, and the results indicate that cell to cell heterogeneity is very relevant for understanding protein production at the population level.

Conclusions

The choice of expression system needs to be evaluated for each specific case, but we believe that the standardized vectors developed for this study can be used to more easily identify the nature of case-specific bottlenecks. By then taking into account the relevant characteristics of each expression cassette it will be easier to make the best choice with respect to the goal of achieving high levels of protein expression in functional or non-functional form.

Translational efficiency of rpoS mRNA from Borrelia burgdorferi: effects of the length and sequence of the mRNA leader region

Regulation of the enzootic cycle in Borrelia burgdorferi requires a shift to the RNA polymerase alternative sigma factor, RpoS. We used in vitro and in vivo assays to assess the relative importance of the putative Shine-Dalgarno sequence and its sequestration for the translational efficiency of rpoS. We created mutant leader regions in which we either removed the Shine-Dalgarno sequence, disrupted the secondary structure or both. Binding assays and toeprint assays demonstrated that both the presence and the availability of the Shine-Dalgarno sequence are important to the efficiency and specificity of ribosome binding. Adding a DsrABb mimic in the form of a single-stranded DNA oligonucleotide increased the level and specificity of binding ribosomes to the transcript with an extended leader, presumably by making the Shine-Dalgarno sequence available for binding. In in vivo assays we confirmed that the Shine-Dalgarno sequence must be both present and un-sequestered in order for translation to proceed efficiently. The longer transcript was significantly better translated in B. burgdorferi at 37 °C than at 26 °C, lending support to the hypothesis that DsrABb acts as a temperature-dependent stimulator of translation. These studies demonstrate that translational regulation of gene expression in B. burgdorferi may be an important mechanism for responding to environmental signals important in the enzootic cycle.

Ribosomal protein S1 unwinds double-stranded RNA in multiple steps

The sequence and secondary structure of the 5'-end of mRNAs regulate translation by controlling ribosome initiation on the mRNA. Ribosomal protein S1 is crucial for ribosome initiation on many natural mRNAs, particularly for those with structured 5'-ends, or with no or weak Shine-Dalgarno sequences. Besides a critical role in translation, S1 has been implicated in several other cellular processes, such as transcription recycling, and the rescuing of stalled ribosomes by tmRNA. The mechanisms of S1 functions are still elusive but have been widely considered to be linked to the affinity of S1 for single-stranded RNA and its corresponding destabilization of mRNA secondary structures. Here, using optical tweezers techniques, we demonstrate that S1 promotes RNA unwinding by binding to the single-stranded RNA formed transiently during the thermal breathing of the RNA base pairs and that S1 dissociation results in RNA rezipping. We measured the dependence of the RNA unwinding and rezipping rates on S1 concentration, and the force applied to the ends of the RNA. We found that each S1 binds 10 nucleotides of RNA in a multistep fashion implying that S1 can facilitate ribosome initiation on structured mRNA by first binding to the single strand next to an RNA duplex structure ("stand-by site") before subsequent binding leads to RNA unwinding. Unwinding by multiple small substeps is much less rate limited by thermal breathing than unwinding in a single step. Thus, a multistep scheme greatly expedites S1 unwinding of an RNA structure compared to a single-step mode.

Kinetic control of translation initiation in bacteria

Translation initiation is a crucial step of protein synthesis which largely defines how the composition of the cellular transcriptome is converted to the proteome and controls the response and adaptation to environmental stimuli. The efficiency of translation of individual mRNAs, and hence the basal shape of the proteome, is defined by the structures of the mRNA translation initiation regions. Initiation efficiency can be regulated by small molecules, proteins, or antisense RNAs, underscoring its importance in translational control. Although initiation has been studied in bacteria for decades, many aspects remain poorly understood. Recent evidence has suggested an unexpected diversity of pathways by which mRNAs can be recruited to the bacterial ribosome, the importance of structural dynamics of initiation intermediates, and the complexity of checkpoints for mRNA selection. In this review, we discuss how the ribosome shapes the landscape of translation initiation by non-linear kinetic processing of the transcriptome information. We summarize the major pathways by which mRNAs enter the ribosome depending on the structure of their 5' untranslated regions, the assembly and the structure of initiation intermediates, the individual and synergistic roles of initiation factors, and the mechanisms of mRNA and initiator tRNA selection.

Robust translation of the nucleoid protein Fis requires a remote upstream AU element and is enhanced by RNA secondary structure

Synthesis of the Fis nucleoid protein rapidly increases in response to nutrient upshifts, and Fis is one of the most abundant DNA binding proteins in Escherichia coli under nutrient-rich growth conditions. Previous work has shown that control of Fis synthesis occurs at transcription initiation of the dusB-fis operon. We show here that while translation of the dihydrouridine synthase gene dusB is low, unusual mechanisms operate to enable robust translation of fis. At least two RNA sequence elements located within the dusB coding region are responsible for high fis translation. The most important is an AU element centered 35 nucleotides (nt) upstream of the fis AUG, which may function as a binding site for ribosomal protein S1. In addition, a 44-nt segment located upstream of the AU element and predicted to form a stem-loop secondary structure plays a prominent role in enhancing fis translation. On the other hand, mutations close to the AUG, including over a potential Shine-Dalgarno sequence, have little effect on Fis protein levels. The AU element and stem-loop regions are phylogenetically conserved within dusB-fis operons of representative enteric bacteria.

Quality control of mRNA decoding on the bacterial ribosome

The ribosome is a major player in providing accurate gene expression in the cell. The fidelity of substrate selection is tightly controlled throughout the translation process, including the initiation, elongation, and termination phases. Although each phase of translation involves different players, that is, translation factors and tRNAs, the general principles of selection appear surprisingly similar for very different substrates. At essentially every step of translation, differences in complex stabilities as well as induced fit are sources of selectivity. A view starts to emerge of how the ribosome uses local and global conformational switches to govern induced-fit mechanisms that ensure fidelity. This review describes the mechanisms of tRNA and mRNA selection at all phases of protein synthesis in bacteria.

Inability of Prevotella bryantii to form a functional Shine-Dalgarno interaction reflects unique evolution of ribosome binding sites in Bacteroidetes

The Shine-Dalgarno (SD) sequence is a key element directing the translation to initiate at the authentic start codons and also enabling translation initiation to proceed in 5′ untranslated mRNA regions (5′-UTRs) containing moderately strong secondary structures. Bioinformatic analysis of almost forty genomes from the major bacterial phylum Bacteroidetes revealed, however, a general absence of SD sequence, drop in GC content and consequently reduced tendency to form secondary structures in 5′-UTRs. The experiments using the Prevotella bryantii TC1-1 expression system were in agreement with these findings: neither addition nor omission of SD sequence in the unstructured 5′-UTR affected the level of the reporter protein, non-specific nuclease NucB. Further, NucB level in P. bryantii TC1-1, contrary to hMGFP level in Escherichia coli, was five times lower when SD sequence formed part of the secondary structure with a folding energy -5,2 kcal/mol. Also, the extended SD sequences did not affect protein levels as in E. coli. It seems therefore that a functional SD interaction does not take place during the translation initiation in P. bryanttii TC1-1 and possibly other members of phylum Bacteroidetes although the anti SD sequence is present in 16S rRNA genes of their genomes. We thus propose that in the absence of the SD sequence interaction, the selection of genuine start codons in Bacteroidetes is accomplished by binding of ribosomal protein S1 to unstructured 5′-UTR as opposed to coding region which is inaccessible due to mRNA secondary structure. Additionally, we found that sequence logos of region preceding the start codons may be used as taxonomical markers. Depending on whether complete sequence logo or only part of it, such as information content and base proportion at specific positions, is used, bacterial genera or families and in some cases even bacterial phyla can be distinguished.

Determinants of translation efficiency and accuracy

A given protein sequence can be encoded by an astronomical number of alternative nucleotide sequences. Recent research has revealed that this flexibility provides evolution with multiple ways to tune the efficiency and fidelity of protein translation and folding.

Proper functioning of biological cells requires that the process of protein expression be carried out with high efficiency and fidelity. Given an amino-acid sequence of a protein, multiple degrees of freedom still remain that may allow evolution to tune efficiency and fidelity for each gene under various conditions and cell types. Particularly, the redundancy of the genetic code allows the choice between alternative codons for the same amino acid, which, although ‘synonymous,' may exert dramatic effects on the process of translation. Here we review modern developments in genomics and systems biology that have revolutionized our understanding of the multiple means by which translation is regulated. We suggest new means to model the process of translation in a richer framework that will incorporate information about gene sequences, the tRNA pool of the organism and the thermodynamic stability of the mRNA transcripts. A practical demonstration of a better understanding of the process would be a more accurate prediction of the proteome, given the transcriptome at a diversity of biological conditions.

The expression of recombinant genes in Escherichia coli can be strongly stimulated at the transcript production level by mutating the DNA-region corresponding to the 5'-untranslated part of mRNA

Secondary structures and the short Shine-Dalgarno sequence in the 5'-untranslated region of bacterial mRNAs (UTR) are known to affect gene expression at the level of translation. Here we report the use of random combinatorial DNA sequence libraries to study UTR function, applying the strong, σ(32)/σ(38)-dependent, and positively regulated Pm promoter as a model. All mutations in the libraries are located at least 8 bp downstream of the transcriptional start site. The libraries were screened using the ampicillin-resistance gene (bla) as reporter, allowing easy identification of UTR mutants that display high levels of expression (up to 20-fold increase relative to the wild-type at the protein level). Studies of the two UTR mutants identified by a modified screening procedure showed that their expression is stimulated to a similar extent at both the transcript and protein product levels. For one such mutant a model analysis of the transcription kinetics showed significant evidence of a difference in the transcription rate (about 18-fold higher than the wild type), while there was no evidence of a difference in transcript stability. The two UTR sequences also stimulated expression from a constitutive σ(70)-dependent promoter (P1/P(anti-tet)), demonstrating that the UTR at the DNA or RNA level has a hitherto unrecognized role in transcription.

Shine-Dalgarno sequence of bacteriophage T4: GAGG prevails in early genes

Translation initiation is governed by a limited number of mRNA sequence motifs within the translation initiation region (TIR). In bacteria and bacteriophages, one of the most important determinants is a Shine-Dalgarno (SD) sequence that base pairs with the anti-SD sequence GAUCACCUCCUUA localized in the 3' end of 16S rRNA. This work assesses a diversity of TIR features in phage T4, focusing on the SD sequence, its spacing to the start codon and relationship to gene expression and essentiality patterns. Analysis shows that GAGG is predominant of all core SD motifs in T4 and its related phages, particularly in early genes. Possible implication of the RegB activity is discussed.

Optimization of synthetic operons using libraries of post-transcriptional regulatory elements

Constructing polycistronic operons is an advantageous strategy for coordinating the expression of -multiple genes in a prokaryotic host. Unfortunately, a basic construct consisting of an inducible promoter and genes cloned in series does not generally lead to optimal results. Here, a combinatorial approach for tuning relative gene expression in operons is presented. The method constructs libraries of post--transcriptional regulatory elements that can be cloned into the noncoding sequence between genes. Libraries can be screened to identify sequences that optimize expression of metabolic pathways, multisubunit proteins, or other situations where precise stoichiometric ratios of proteins are desired.

Translation on demand by a simple RNA-based thermosensor

Structured RNA regions are important gene control elements in prokaryotes and eukaryotes. Here, we show that the mRNA of a cyanobacterial heat shock gene contains a built-in thermosensor critical for photosynthetic activity under stress conditions. The exceptionally short 5′-untranslated region is comprised of a single hairpin with an internal asymmetric loop. It inhibits translation of the Synechocystis hsp17 transcript at normal growth conditions, permits translation initiation under stress conditions and shuts down Hsp17 production in the recovery phase. Point mutations that stabilized or destabilized the RNA structure deregulated reporter gene expression in vivo and ribosome binding in vitro. Introduction of such point mutations into the Synechocystis genome produced severe phenotypic defects. Reversible formation of the open and closed structure was beneficial for viability, integrity of the photosystem and oxygen evolution. Continuous production of Hsp17 was detrimental when the stress declined indicating that shutting-off heat shock protein production is an important, previously unrecognized function of RNA thermometers. We discovered a simple biosensor that strictly adjusts the cellular level of a molecular chaperone to the physiological need.

Mathematical modeling of translation initiation for the estimation of its efficiency to computationally design mRNA sequences with desired expression levels in prokaryotes

Background

Within the emerging field of synthetic biology, engineering paradigms have recently been used to design biological systems with novel functionalities. One of the essential challenges hampering the construction of such systems is the need to precisely optimize protein expression levels for robust operation. However, it is difficult to design mRNA sequences for expression at targeted protein levels, since even a few nucleotide modifications around the start codon may alter translational efficiency and dramatically (up to 250-fold) change protein expression. Previous studies have used ad hoc approaches (e.g., random mutagenesis) to obtain the desired translational efficiencies for mRNA sequences. Hence, the development of a mathematical methodology capable of estimating translational efficiency would greatly facilitate the future design of mRNA sequences aimed at yielding desired protein expression levels.

Results

We herein propose a mathematical model that focuses on translation initiation, which is the rate-limiting step in translation. The model uses mRNA-folding dynamics and ribosome-binding dynamics to estimate translational efficiencies solely from mRNA sequence information. We confirmed the feasibility of our model using previously reported expression data on the MS2 coat protein. For further confirmation, we used our model to design 22 luxR mRNA sequences predicted to have diverse translation efficiencies ranging from 10^-5 to 1. The expression levels of these sequences were measured in Escherichia coli and found to be highly correlated (R² = 0.87) with their estimated translational efficiencies. Moreover, we used our computational method to successfully transform a low-expressing DsRed2 mRNA sequence into a high-expressing mRNA sequence by maximizing its translational efficiency through the modification of only eight nucleotides upstream of the start codon.

Conclusions

We herein describe a mathematical model that uses mRNA sequence information to estimate translational efficiency. This model could be used to design best-fit mRNA sequences having a desired protein expression level, thereby facilitating protein over-production in biotechnology or the protein expression-level optimization necessary for the construction of robust networks in synthetic biology.

Mechanism of regulating the expression of lambdaN gene by ribosomal protein at translational level

In ribosomal protein S12 mutant or L24 mutant the expression of lambdaN gene was depressed at translational level. To study its mechanism the lambdaN gene region of lambdaN -lacZ gene fusion was trimmed from its 5' end to 3' end with DNA exonuclease III (DNA exoIII) in order to alter the TIR (translational initiation region) and the coding region of lambdaN gene. After DNA sequencing 23 species of different lambdaN-lacZ fused genes were obtained. The beta-galactosidase activities of these deletants in ribosomal protein mutant were compared with that in wild type strain. The result indicated that (i) S12 mutant could affect 305 subunit's binding to the TIR of lambdaN gene messenger and cause the difficulty in forming 30s initiation complex and then decrease the efficiency of translational initiation; (ii) in S12 mutant the coding region of lambdaN gene also affected the expression lambdaN gene; (iii) in L24 mutant the inhibition of lambdaN gene expression was not related to translational initiation and the 5' end of the coding region of lambdaN gene, but related to the 3' end of lambdaN gene.

Translational activation of rpoS mRNA by the non-coding RNA DsrA and Hfq does not require ribosome binding

At low temperature, translational activation of rpoS mRNA, encoding the stationary phase sigma-factor, σ^S, involves the small regulatory RNA (sRNA) DsrA and the RNA chaperone Hfq. The Hfq-mediated DsrA-rpoS interaction relieves an intramolecular secondary structure that impedes ribosome access to the rpoS ribosome binding site. In addition, DsrA/rpoS duplex formation creates an RNase III cleavage site within the duplex. Previous biochemical studies suggested that DsrA and Hfq associate with the 30S ribosomal subunit protein S1, which implied a role for the ribosome in sRNA-mediated post-transcriptional regulation. Here, we show by ribosome profiling that Hfq partitions with the cytoplasmic fraction rather than with 30S subunits. Besides, by employing immunological techniques, no evidence for a physical interaction between Hfq and S1 was obtained. Similarly, in vitro studies did not reveal a direct interaction between DsrA and S1. By employing a ribosome binding deficient rpoS mRNA, and by using the RNase III clevage in the DsrA/rpoS duplex as a diagnostic marker, we provide in vivo evidence that the Hfq-mediated DsrA/rpoS interaction, and consequently the structural changes in rpoS mRNA precede ribosome binding. These data suggest a simple mechanistic model in which translational activation by DsrA provides a translationally competent rpoS mRNA to which 30S subunits can readily bind.

Adaptation of the short intergenic spacers between co-directional genes to the Shine-Dalgarno motif among prokaryote genomes

Background

In prokaryote genomes most of the co-directional genes are in close proximity. Even the coding sequence or the stop codon of a gene can overlap with the Shine-Dalgarno (SD) sequence of the downstream co-directional gene. In this paper we analyze how the presence of SD may influence the stop codon usage or the spacing lengths between co-directional genes.

Results

The SD sequences for 530 prokaryote genomes have been predicted using computer calculations of the base-pairing free energy between translation initiation regions and the 16S rRNA 3' tail. Genomes with a large number of genes with the SD sequence concentrate this regulatory motif from 4 to 11 bps before the start codon. However, not all genes seem to have the SD sequence. Genes separated from 1 to 4 bps from a co-directional upstream gene show a high SD presence, though this regulatory signal is located towards the 3' end of the coding sequence of the upstream gene. Genes separated from 9 to 15 bps show the highest SD presence as they accommodate the SD sequence within an intergenic region. However, genes separated from around 5 to 8 bps have a lower percentage of SD presence and when the SD is present, the stop codon usage of the upstream gene changes to accommodate the overlap between the SD sequence and the stop codon.

Conclusion

The SD presence makes the intergenic lengths from 5 to 8 bps less frequent and causes an adaptation of the stop codon usage. Our results introduce new elements to the discussion of which factors affect the intergenic lengths, which cannot be totally explained by the pressure to compact the prokaryote genomes.

Ribosome-binding site interference caused by Shine-Dalgarno-like nucleotide sequences in Escherichia coli cells

Two-cistronic expression plasmids are useful for high-level expression of heterologous genes in Escherichia coli cells by preventing the inhibition of translational initiation. In the process of constructing a two-cistronic expression plasmid pCbSTCR-4 containing the fragments of the porcine cytochrome b(5) (Psb5) and NADPH-cytochrome P450 reductase (PsCPR) genes as the first and second cistrons, respectively, the presence of a specific region in the first cistron that lowered the accumulation level of the PsCPR was suggested [Kimura, S., et al. (2005) J. Biochem. 137, 523-533]. In this study, a disturbing nucleotide sequence similar to a Shine-Dalgarno (SD) sequence (SD-like sequence), AGGAG, was identified at the 5'-upstream region near the SD sequence for the second cistron. Silent mutations in the SD-like sequence that lowered the similarity to a typical SD sequence increased the accumulation level of PsCPR. SD-like sequences introduced into mono-cistronic expression plasmids for the Psb5 and PsCPR genes also decreased the accumulation level of these proteins. The SD-like sequence also decreased the accumulation level of the insoluble PsCPR protein. This type of ribosome-binding site interference is useful not only for precise control of protein accumulation but also for increasing the soluble form of recombinant proteins in E. coli cells.