Quantitative analysis of ribosome binding sites in E.coli

Controlled production of lipopolysaccharides increases immune activation in Salmonella treatments of cancer

Immunotherapies have revolutionized cancer treatment. These treatments rely on immune cell activation in tumours, which limits the number of patients that respond. Inflammatory molecules, like lipopolysaccharides (LPS), can activate innate immune cells, which convert tumour microenvironments from cold to hot, and increase therapeutic efficacy. However, systemic delivery of lipopolysaccharides (LPS) can induce cytokine storm. In this work, we developed immune-controlling Salmonella (ICS) that only produce LPS in tumours after colonization and systemic clearance. We tuned the expression of msbB, which controls production of immunogenic LPS, by optimizing its ribosomal binding sites and protein degradation tags. This genetic system induced a controllable inflammatory response and increased dendritic cell cross-presentation in vitro. The strong off state did not induce TNFα production and prevented adverse events when injected into mice. The accumulation of ICS in tumours after intravenous injection focused immune responses specifically to tumours. Tumour-specific expression of msbB increased infiltration of immune cells, activated monocytes and neutrophils, increased tumour levels of IL-6, and activated CD8 T cells in draining lymph nodes. These immune responses reduced tumour growth and increased mouse survival. By increasing the efficacy of bacterial anti-cancer therapy, localized production of LPS could provide increased options to patients with immune-resistant cancers.

The Influence of the Nucleotide Composition of Genes and Gene Regulatory Elements on the Efficiency of Protein Expression in Escherichia coli

Recombinant proteins expressed in Escherichia coli are widely used in biochemical research and industrial processes. At the same time, achieving higher protein expression levels and correct protein folding still remains the key problem, since optimization of nutrient media, growth conditions, and methods for induction of protein synthesis do not always lead to the desired result. Often, low protein expression is determined by the sequences of the expressed genes and their regulatory regions. The genetic code is degenerated; 18 out of 20 amino acids are encoded by more than one codon. Choosing between synonymous codons in the coding sequence can significantly affect the level of protein expression and protein folding due to the influence of the gene nucleotide composition on the probability of formation of secondary mRNA structures that affect the ribosome binding at the translation initiation phase, as well as the ribosome movement along the mRNA during elongation, which, in turn, influences the mRNA degradation and the folding of the nascent protein. The nucleotide composition of the mRNA untranslated regions, in particular the promoter and Shine-Dalgarno sequences, also affects the efficiency of mRNA transcription, translation, and degradation. In this review, we describe the genetic principles that determine the efficiency of protein production in Escherichia coli.

A Kinetic Framework for Modeling Oleochemical Biosynthesis in E. coli

Microorganisms build fatty acids with biocatalytic assembly lines, or fatty acid synthases (FASs), that can be repurposed to produce a broad set of fuels and chemicals. Despite their versatility, the product profiles of FAS-based pathways are challenging to adjust without experimental iteration, and off-target products are common. This study uses a detailed kinetic model of the E. coli FAS as a foundation to model nine oleochemical pathways. These models provide good fits to experimental data and help explain unexpected results from in vivo studies. An analysis of pathways for alkanes and fatty acid ethyl esters, for example, suggests that reductions in titer caused by enzyme overexpression-an experimentally consistent phenomenon-can result from shifts in metabolite pools that are incompatible with the substrate specificities of downstream enzymes, and a focused examination of multiple alcohol pathways indicates that coordinated shifts in enzyme concentrations provide a general means of tuning the product profiles of pathways with promiscuous components. The study concludes by integrating all models into a graphical user interface. The models supplied by this work provide a versatile kinetic framework for studying oleochemical pathways in different biochemical contexts. This article is protected by copyright. All rights reserved.

Deciphering the Rules of Ribosome Binding Site Differentiation in Context Dependence

The ribosome binding site (RBS) is a crucial element regulating translation. However, the activity of RBS is poorly predictable, because it is strongly affected by the local possible secondary structure, that is, context dependence. By the Flowseq technique, over 20 000 RBS variants were sorted and sequenced, and the translation of multiple genes under the same RBS was quantitatively characterized to evaluate the context dependence of each RBS variant in E. coli. Two regions, (-7 to -2) and (-17 to -12), of RBS were predicted with a higher possibility to pair with each other to slow down the translation initiation. Associations between phenotypes and the intrinsic factors suspected to affect translation efficiency and context dependence of the RBS, including nucleotide bias at each position, free energy, and conservation, were disentangled. The results showed that translation efficiency was influenced more significantly by conservation of the SD region (-16 to -8), while an AC-rich spacer region (-7 to -1) was associated with low context dependence. We confirmed these characteristics using a series of synthesized RBSs. The average correlation between multiple reporters was significantly higher for RBSs with an AC-rich spacer (0.714) compared with a GU-rich spacer (0.286). Overall, we proposed general design criteria to improve programmability and minimize context dependence of RBS. The characteristics unraveled here can be adapted to other bacteria for fine-tuning target-gene expression.

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.

Reducing metabolic burden in the PACEmid evolver system by remastering high-copy phagemid vectors

Orthogonal or non-cross-reacting transcription factors are used in synthetic biology as components of genetic circuits. Brödel et al. (2016) engineered 12 such cIλ transcription factor variants using a directed evolution 'PACEmid' system. The variants operate as dual activator/repressors and expand gene circuit construction possibilities. However, the high-copy phagemid vectors carrying the cIλ variants imposed high metabolic burden upon cells. Here, the authors 'remaster' the phagemid backbones to relieve their burden substantially, exhibited by a recovery in Escherichia coli growth. The remastered phagemids' ability to function within the PACEmid evolver system is maintained, as is the cIλ transcription factors' activity within these vectors. The low-burden phagemid versions are more suitable for use in PACEmid experiments and synthetic gene circuits; the authors have, therefore, replaced the original high-burden phagemids on the Addgene repository. The authors' work emphasises the importance of understanding metabolic burden and incorporating it into design steps in future synthetic biology ventures.

StartLink and StartLink+: Prediction of Gene Starts in Prokaryotic Genomes

State-of-the-art algorithms of ab initio gene prediction for prokaryotic genomes were shown to be sufficiently accurate. A pair of algorithms would agree on predictions of gene 3'ends. Nonetheless, predictions of gene starts would not match for 15-25% of genes in a genome. This discrepancy is a serious issue that is difficult to be resolved due to the absence of sufficiently large sets of genes with experimentally verified starts. We have introduced StartLink that infers gene starts from conservation patterns revealed by multiple alignments of homologous nucleotide sequences. We also have introduced StartLink+ combining both ab initio and alignment-based methods. The ability of StartLink to predict the start of a given gene is restricted by the availability of homologs in a database. We observed that StartLink made predictions for 85% of genes per genome on average. The StartLink+ accuracy was shown to be 98-99% on the sets of genes with experimentally verified starts. In comparison with database annotations, we observed that the annotated gene starts deviated from the StartLink+ predictions for ∼5% of genes in AT-rich genomes and for 10-15% of genes in GC-rich genomes on average. The use of StartLink+ has a potential to significantly improve gene start annotation in genomic databases.

Dynamic Control of Gene Expression with Riboregulated Switchable Feedback Promoters

One major challenge in synthetic biology is the deleterious impacts of cellular stress caused by expression of heterologous pathways, sensors, and circuits. Feedback control and dynamic regulation are broadly proposed strategies to mitigate this cellular stress by optimizing gene expression levels temporally and in response to biological cues. While a variety of approaches for feedback implementation exist, they are often complex and cannot be easily manipulated. Here, we report a strategy that uses RNA transcriptional regulators to integrate additional layers of control over the output of natural and engineered feedback responsive circuits. Called riboregulated switchable feedback promoters (rSFPs), these gene expression cassettes can be modularly activated using multiple mechanisms, from manual induction to autonomous quorum sensing, allowing control over the timing, magnitude, and autonomy of expression. We develop rSFPs in Escherichia coli to regulate multiple feedback networks and apply them to control the output of two metabolic pathways. We envision that rSFPs will become a valuable tool for flexible and dynamic control of gene expression in metabolic engineering, biological therapeutic production, and many other applications.

New Set of Yeast Vectors for Shuttle Expression in Escherichia coli

Promoters that play an essential role in the gene regulation are of particular interest to the synthetic biology communities. Recent advances in high-throughput DNA sequencing have greatly increased the breadth of new genetic parts. The development of promoters with broad host properties could enable rapid phenotyping of genetic constructs in different hosts. In this study, we discovered that the GAL1/10 bidirectional promoter from Saccharomyces cerevisiae could be used for shuttle expression in Escherichia coli. Further investigation revealed that the GAL1/10 bidirectional promoter is subjected to catabolite repression in E. coli. We next constructed a set of Golden-Gate assembly vectors for shuttle expression between S. cerevisiae and E. coli. The utility of shuttle vectors was demonstrated for rapid phenotyping of a multigene pathway for cinnamyl alcohol production. Taken together, our work opens a new avenue for the future development of broad host expression systems between prokaryotic and eukaryotic hosts.

Fine-Tuning Multi-Gene Clusters via Well-Characterized Gene Expression Regulatory Elements: Case Study of the Arginine Synthesis Pathway in C. glutamicum

Promoters and ribosome binding sites (RBSs) are routinely applied in gene expression regulation, but their orthogonality and combinatorial effects have not yet been systematically studied in Corynebacterium glutamicum. Here, 17 core promoters and 29 RBSs in C. glutamicum were characterized, which exhibited 470-fold and 430-fold in transcriptional and translational activity, respectively. By comparing the expression of two reporter genes regulated by multiple RBSs, the RBS efficacy showed significant dependence on the gene context, besides the RBSs' strength, reflecting the poor orthogonality of RBSs. Bicistron-modified RBS (referred as bc-RBS) was adapted to C. glutamicum, which improved RBS reliability. By coupling a series of promoters with RBSs/bc-RBSs, a much broader regulation range that spanned 4 orders of magnitude was observed compared with that of a sole element, and the contribution to gene expression of RBS was more than that of promoter. Finally, promoters and RBSs were applied as built-in elements to fine-tune the gene cluster in the arginine synthesis pathway in C. glutamicum. Compared with the original strain, more arginine (1.61-fold) or citrulline (2.35-fold) was accumulated in a 7 L bioreactor by strains with the gene expression regulation system rationally engineered. We demonstrated that, via combination of well-characterized gene elements, and overall consideration for both transcription and translation, the biosynthesis pathway can be effectively balanced, and the yield of a target metabolite can be further improved.

Tunable Gene Expression System Independent of Downstream Coding Sequence

Fine control of the expression levels of proteins constitutes a major challenge in synthetic biology and metabolic engineering. However, the dependence of translation initiation on the downstream coding sequence (CDS) obscures accurate prediction of the protein expression levels from mRNA sequences. Here, we present a tunable gene-expression system comprising 24 expression cassettes that produce predefined relative expression levels of proteins ranging from 0.001 to 1 without being influenced by the downstream CDS. To validate the practical utility of the tunable expression system, it was applied to a synthetic circuit displaying three states of fluorescence depending on the difference in protein expression levels. To demonstrate the suitability of application to metabolic engineering, this system was used to diversify the levels of key metabolic enzymes. As a result, expression-optimized strains were capable of producing 2.25 g/L of cadaverine, 2.59 g/L of L-proline, and 95.7 mg/L of 1-propanol. Collectively, the tunable expression system could be utilized to optimize genetic circuits for desired operation and to optimize metabolic fluxes through biosynthetic pathways for enhancing production yields of bioproducts. This tunable system will be useful for studying basic and applied biological sciences in addition to applications in synthetic biology and metabolic engineering.

pheS AG Based Rapid and Efficient Markerless Mutagenesis in Methylotuvimicrobium

Due to their fast growth rate and robustness, some haloalkalitolerant methanotrophs from the genus Methylotuvimicrobium have recently become not only promising biocatalysts for methane conversion but also favorable materials for obtaining fundamental knowledge on methanotrophs. Here, to realize unmarked genome modification in Methylotuvimicrobium bacteria, a counterselectable marker (CSM) was developed based on pheS, which encodes the α-subunit of phenylalanyl-tRNA synthetase. Two-point mutations (T252A and A306G) were introduced into PheS in Methylotuvimicrobium buryatense 5GB1C, generating PheS AG , which can recognize p-chloro-phenylalanine (p-Cl-Phe) as a substrate. Theoretically, the expression of PheS AG in a cell will result in the incorporation of p-Cl-Phe into proteins, leading to cell death. The P tac promoter and the ribosome-binding site region of mmoX were employed to control pheS AG , producing the pheS AG -3 CSM. M. buryatense 5GB1C harboring pheS AG -3 was extremely sensitive to 0.5 mM p-Cl-Phe. Then, a positive and counterselection cassette, PZ (only 1.5 kb in length), was constructed by combining pheS AG -3 and the zeocin resistance gene. A PZ- and PCR-based strategy was used to create the unmarked deletion of glgA1 or the whole smmo operon in M. buryatense 5GB1C and Methylotuvimicrobium alcaliphilum 20Z. The positive rates were over 92%, and the process could be accomplished in as few as eight days.

Restriction enzymes use a 24 dimensional coding space to recognize 6 base long DNA sequences

Restriction enzymes recognize and bind to specific sequences on invading bacteriophage DNA. Like a key in a lock, these proteins require many contacts to specify the correct DNA sequence. Using information theory we develop an equation that defines the number of independent contacts, which is the dimensionality of the binding. We show that EcoRI, which binds to the sequence GAATTC, functions in 24 dimensions. Information theory represents messages as spheres in high dimensional spaces. Better sphere packing leads to better communications systems. The densest known packing of hyperspheres occurs on the Leech lattice in 24 dimensions. We suggest that the single protein EcoRI molecule employs a Leech lattice in its operation. Optimizing density of sphere packing explains why 6 base restriction enzymes are so common.

Within-Gene Shine-Dalgarno Sequences Are Not Selected for Function

The Shine-Dalgarno (SD) sequence motif facilitates translation initiation and is frequently found upstream of bacterial start codons. However, thousands of instances of this motif occur throughout the middle of protein coding genes in a typical bacterial genome. Here, we use comparative evolutionary analysis to test whether SD sequences located within genes are functionally constrained. We measure the conservation of SD sequences across Enterobacteriales, and find that they are significantly less conserved than expected. Further, the strongest SD sequences are the least conserved whereas we find evidence of conservation for the weakest possible SD sequences given amino acid constraints. Our findings indicate that most SD sequences within genes are likely to be deleterious and removed via selection. To illustrate the origin of these deleterious costs, we show that ATG start codons are significantly depleted downstream of SD sequences within genes, highlighting the constraint that these sequences impose on the surrounding nucleotides to minimize the potential for erroneous translation initiation.

Identification and Mutational Analysis of Escherichia coli Sorbitol-Enhanced Glucose-Repressed srlA Promoter Expressed in LB Medium by Using Homologous Recombination and One-Round PCR Products

Escherichia coli has been used for recombinant protein production for many years. However, no native E. coli promoters have been found for constitutive expression in LB medium. To obtain high-expression E. coli promoters active in LB medium, we inserted various promoter regions upstream of eEmRFP that encodes a red fluorescent protein. Among the selected promoters, only colonies of srlA promoter transformants turned red on LB plate. srlA is a gene that regulates sorbitol utilization. The addition of sorbitol enhanced eEmRFP expression but glucose and other sugars repressed, indicating that srlAp is a sorbitol-enhanced glucose-repressed promoter. To analyze the srlAp sequence, a novel site-directed mutagenesis method was developed. Since we demonstrated that homologous recombination in E. coli could occur between 12-bp sequences, 12-bp overlapping sequences were attached to the set of primers that were designed to produce a full-length plasmid, denoted "one-round PCR product." Using this method, we identified that the srlA promoter region was 100 bp. Further, the sequence adjacent to the start codon was found to be essential for high expression, suggesting that the traditionally used restriction enzyme sites for cloning in the promoter region have hindered expression. The srlA-driven expression system and DNA manipulation with one-round PCR products are useful tools in E. coli genetic engineering.

Genetic analysis of the basis of translation in the -1 frame of an unusual non-ORF sequence isolated from phage display

An unusual peptide-encoding sequence, called H10, and several derivatives of this sequence were previously isolated from a random peptide library screened by phage display during drug discovery protocols. The H10 family of sequences had the unusual property of being expressed despite the absence of an open reading frame. When these sequences were fused to a reporter lacZ gene in all three frames, beta-galactosidase was expressed not only from the parental non-open reading frame, consistent with the original isolations, but also from the frame -1 to the parental. This unexpected translation in a second reading frame could result from either a recoding event or from an internal translation initiation event. In order to elucidate which type of event, a genetic approach was selected to eliminate a potential downstream initiator site within the H10 sequence. This report provides strong evidence that translation in the -1 frame in this family of sequences is indeed originating from a downstream translation initiation event. Unexpectedly, the mutation eliminating the downstream initiation event in the -1 frame simultaneously elevated expression in the original non-open reading frame.

Modeling leaderless transcription and atypical genes results in more accurate gene prediction in prokaryotes

In a conventional view of the prokaryotic genome organization, promoters precede operons and ribosome binding sites (RBSs) with Shine-Dalgarno consensus precede genes. However, recent experimental research suggesting a more diverse view motivated us to develop an algorithm with improved gene-finding accuracy. We describe GeneMarkS-2, an ab initio algorithm that uses a model derived by self-training for finding species-specific (native) genes, along with an array of precomputed “heuristic” models designed to identify harder-to-detect genes (likely horizontally transferred). Importantly, we designed GeneMarkS-2 to identify several types of distinct sequence patterns (signals) involved in gene expression control, among them the patterns characteristic for leaderless transcription as well as noncanonical RBS patterns. To assess the accuracy of GeneMarkS-2, we used genes validated by COG (Clusters of Orthologous Groups) annotation, proteomics experiments, and N-terminal protein sequencing. We observed that GeneMarkS-2 performed better on average in all accuracy measures when compared with the current state-of-the-art gene prediction tools. Furthermore, the screening of ∼5000 representative prokaryotic genomes made by GeneMarkS-2 predicted frequent leaderless transcription in both archaea and bacteria. We also observed that the RBS sites in some species with leadered transcription did not necessarily exhibit the Shine-Dalgarno consensus. The modeling of different types of sequence motifs regulating gene expression prompted a division of prokaryotic genomes into five categories with distinct sequence patterns around the gene starts.

70ProPred: a predictor for discovering sigma70 promoters based on combining multiple features

BACKGROUND

Promoter is an important sequence regulation element, which is in charge of gene transcription initiation. In prokaryotes, σ70 promoters regulate the transcription of most genes. The promoter recognition has been a crucial part of gene structure recognition. It's also the core issue of constructing gene transcriptional regulation network. With the successfully completion of genome sequencing from an increasing number of microbe species, the accurate identification of σ70 promoter regions in DNA sequence is not easy.

RESULTS

In order to improve the prediction accuracy of sigma70 promoters in prokaryote, a promoter recognition model 70ProPred was established. In this work, two sequence-based features, including position-specific trinucleotide propensity based on single-stranded characteristic (PSTNPss) and electron-ion potential values for trinucleotides (PseEIIP), were assessed to build the best prediction model. It was found that 79 features of PSTNPSS combined with 64 features of PseEIIP obtained the best performance for sigma70 promoter identification, with a promising accuracy and the Matthews correlation coefficient (MCC) at 95.56% and 0.90, respectively.

CONCLUSION

The jackknife tests showed that 70ProPred outperforms the existing sigma70 promoter prediction approaches in terms of accuracy and stability. Additionally, this approach can also be extended to predict promoters of other species. In order to facilitate experimental biologists, an online web server for the proposed method was established, which is freely available at http://server.malab.cn/70ProPred/ .

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.

The sequence preference of DNA cleavage by T4 endonuclease VII

The enzyme T4 endonuclease VII is a resolvase that acts on branched DNA intermediates during genetic recombination, by cleaving DNA with staggered cuts approximately 3-6 bp apart. In this paper, we investigated the sequence preference of this cleavage reaction utilising two different DNA sequences. For the first time, the DNA sequence preference of T4 endonuclease VII cleavage sites has been examined without the presence of a known DNA substrate to mask any inherent nucleotide preference. The use of the ABI3730 platform enables the cleavage site to be determined at nucleotide resolution. We found that T4 endonuclease VII cleaves DNA with a sequence preference. We calculated the frequency of nucleotides surrounding the cleavage sites and found that following nucleotides had the highest incidence: AWTAN*STC, where N* indicates the cleavage site between positions 0 and 1, N is any base, W is A or T, and S is G or C. An A at position -1 and T at position +2 were the most predominant nucleotides at the cleavage site. Using a Sequence Logo method, the sequence TATTAN*CT was derived at the cleavage site. Note that A and T nucleotides were highly preferred 5' to the cleavage sites in both methods of analysis. It was proposed that the enzyme recognises the narrower minor groove of these consecutive AT base pairs and cleaves DNA 3' to this feature.

Measurements of translation initiation from all 64 codons in E. coli

Our understanding of translation underpins our capacity to engineer living systems. The canonical start codon (AUG) and a few near-cognates (GUG, UUG) are considered as the ‘start codons’ for translation initiation in Escherichia coli. Translation is typically not thought to initiate from the 61 remaining codons. Here, we quantified translation initiation of green fluorescent protein and nanoluciferase in E. coli from all 64 triplet codons and across a range of DNA copy number. We detected initiation of protein synthesis above measurement background for 47 codons. Translation from non-canonical start codons ranged from 0.007 to 3% relative to translation from AUG. Translation from 17 non-AUG codons exceeded the highest reported rates of non-cognate codon recognition. Translation initiation from non-canonical start codons may contribute to the synthesis of peptides in both natural and synthetic biological systems.

Mechanisms for Differential Protein Production in Toxin-Antitoxin Systems

Toxin-antitoxin (TA) systems are key regulators of bacterial persistence, a multidrug-tolerant state found in bacterial species that is a major contributing factor to the growing human health crisis of antibiotic resistance. Type II TA systems consist of two proteins, a toxin and an antitoxin; the toxin is neutralized when they form a complex. The ratio of antitoxin to toxin is significantly greater than 1.0 in the susceptible population (non-persister state), but this ratio is expected to become smaller during persistence. Analysis of multiple datasets (RNA-seq, ribosome profiling) and results from translation initiation rate calculators reveal multiple mechanisms that ensure a high antitoxin-to-toxin ratio in the non-persister state. The regulation mechanisms include both translational and transcriptional regulation. We classified E. coli type II TA systems into four distinct classes based on the mechanism of differential protein production between toxin and antitoxin. We find that the most common regulation mechanism is translational regulation. This classification scheme further refines our understanding of one of the fundamental mechanisms underlying bacterial persistence, especially regarding maintenance of the antitoxin-to-toxin ratio.

The role of mRNA structure in bacterial translational regulation

The characteristics of bacterial messenger RNAs (mRNAs) that influence translation efficiency provide many convenient handles for regulation of gene expression, especially when coupled with the processes of transcription termination and mRNA degradation. An mRNA's structure, especially near the site of initiation, has profound consequences for how readily it is translated. This property allows bacterial gene expression to be altered by changes to mRNA structure induced by temperature, or interactions with a wide variety of cellular components including small molecules, other RNAs (such as sRNAs and tRNAs), and RNA-binding proteins. This review discusses the links between mRNA structure and translation efficiency, and how mRNA structure is manipulated by conditions and signals within the cell to regulate gene expression. The range of RNA regulators discussed follows a continuum from very complex tertiary structures such as riboswitch aptamers and ribosomal protein-binding sites to thermosensors and mRNA:sRNA interactions that involve only base-pairing interactions. Furthermore, the high degrees of diversity observed for both mRNA structures and the mechanisms by which inhibition of translation occur have significant consequences for understanding the evolution of bacterial translational regulation. WIREs RNA 2017, 8:e1370. doi: 10.1002/wrna.1370 For further resources related to this article, please visit the WIREs website.

Sequence-based prediction of transcription upregulation by auxin in plants

Auxin is one of the main regulators of growth and development in plants. Prediction of auxin response based on gene sequence is of high importance. We found the TGTCNC consensus of 111 known natural and artificially mutated auxin response elements (AuxREs) with measured auxin-caused relative increase in genes' transcription levels, so-called either "a response to auxin" or "an auxin response." This consensus was identical to the most cited AuxRE motif. Also, we found several DNA sequence features that correlate with auxin-caused increase in genes' transcription levels, namely: number of matches with TGTCNC, homology score based on nucleotide frequencies at the consensus positions, abundances of five trinucleotides and five B-helical DNA features around these known AuxREs. We combined these correlations using a four-step empirical model of auxin response based on a gene's sequence with four steps, namely: (1) search for AuxREs with no auxin; (2) stop at the found AuxRE; (3) repression of the basal transcription of the gene having this AuxRE; and (4) manifold increase of this gene's transcription in response to auxin. Independently measured increases in transcription levels in response to auxin for 70 Arabidopsis genes were found to significantly correlate with predictions of this equation (r = 0.44, p < 0.001) as well as with TATA-binding protein (TBP)'s affinity to promoters of these genes and with nucleosome packing of these promoters (both, p < 0.025). Finally, we improved our equation for prediction of a gene's transcription increase in response to auxin by taking into account TBP-binding and nucleosome packing (r = 0.53, p < 10(-6)). Fisher's F-test validated the significant impact of both TBP/promoter-affinity and promoter nucleosome on auxin response in addition to those of AuxRE, F = 4.07, p < 0.025. It means that both TATA-box and nucleosome should be taken into account to recognize transcription factor binding sites upon DNA sequences: in the case of the TATA-less nucleosome-rich promoters, recognition scores must be higher than in the case of the TATA-containing nucleosome-free promoters at the same transcription activity.

BASIC: A New Biopart Assembly Standard for Idempotent Cloning Provides Accurate, Single-Tier DNA Assembly for Synthetic Biology

The ability to quickly and reliably assemble DNA constructs is one of the key enabling technologies for synthetic biology. Here we define a new Biopart Assembly Standard for Idempotent Cloning (BASIC), which exploits the principle of orthogonal linker based DNA assembly to define a new physical standard for DNA parts. Further, we demonstrate a new robust method for assembly, based on type IIs restriction enzyme cleavage and ligation of oligonucleotides with single stranded overhangs that determine the assembly order. It allows for efficient, parallel assembly with great accuracy: 4 part assemblies achieved 93% accuracy with single antibiotic selection and 99.7% accuracy with double antibiotic selection, while 7 part assemblies achieved 90% accuracy with double antibiotic selection. The linkers themselves may also be used as composable parts for RBS tuning or the creation of fusion proteins. The standard has one forbidden restriction site and provides for an idempotent, single tier organization, allowing all parts and composite constructs to be maintained in the same format. This makes the BASIC standard conceptually simple at both the design and experimental levels.

Quantitative analysis of mammalian translation initiation sites by FACS-seq

An approach combining fluorescence-activated cell sorting and high-throughput DNA sequencing (FACS-seq) was employed to determine the efficiency of start codon recognition for all possible translation initiation sites (TIS) utilizing AUG start codons. Using FACS-seq, we measured translation from a genetic reporter library representing all 65,536 possible TIS sequences spanning the −6 to +5 positions. We found that the motif RYMRMVAUGGC enhanced start codon recognition and translation efficiency. However, dinucleotide interactions, which cannot be conveyed by a single motif, were also important for modeling TIS efficiency. Our dataset combined with modeling allowed us to predict genome-wide translation initiation efficiency for all mRNA transcripts. Additionally, we screened somatic TIS mutations associated with tumorigenesis to identify candidate driver mutations consistent with known tumor expression patterns. Finally, we implemented a quantitative leaky scanning model to predict alternative initiation sites that produce truncated protein isoforms and compared predictions with ribosome footprint profiling data. The comprehensive analysis of the TIS sequence space enables quantitative predictions of translation initiation based on genome sequence.

Production of recombinant protein by a novel oxygen-induced system in Escherichia coli

Background

The SoxRS regulon of E. coli is activated in response to elevated dissolved oxygen concentration likely to protect the bacteria from possible oxygen damage. The soxS expression can be increased up to 16 fold, making it a possible candidate for recombinant protein expression. Compared with the existing induction approaches, oxygen induction is advantageous because it does not involve addition or depletion of growth factors or nutrients, addition of chemical inducers or temperature changes that can affect growth and metabolism of the producing bacteria. It also does not affect the composition of the growth medium simplifying the recovery and purification processes.

Results

The soxS promoter was cloned into the commercial pGFPmut3.1 plasmid creating pAB49, an expression vector that can be induced by increasing oxygen concentration. The efficiency and the regulatory properties of the soxS promoter were characterized by measuring the GFP expression when the culture dissolved oxygen concentration was increased from 30% to 300% air saturation. The expression level of recombinant GFP was proportional to the oxygen concentration, demonstrating that pAB49 is a controllable expression vector. A possible harmful effect of elevated oxygen concentration on the recombinant product was found to be negligible by determining the protein-carbonyl content and its specific fluorescence.

By performing high density growth in modified LB medium, the cells were induced by increasing the oxygen concentration. After 3 hours at 300% air saturation, GFP fluorescence reached 109000 FU (494 mg of GFP/L), representing 3.4% of total protein, and the cell concentration reached 29.1 g/L (DW).

Conclusions

Induction of recombinant protein expression by increasing the dissolved oxygen concentration was found to be a simple and efficient alternative expression strategy that excludes the use of chemical, nutrient or thermal inducers that have a potential negative effect on cell growth or the product recovery.

One step DNA assembly for combinatorial metabolic engineering

The rapid and efficient assembly of multi-step metabolic pathways for generating microbial strains with desirable phenotypes is a critical procedure for metabolic engineering, and remains a significant challenge in synthetic biology. Although several DNA assembly methods have been developed and applied for metabolic pathway engineering, many of them are limited by their suitability for combinatorial pathway assembly. The introduction of transcriptional (promoters), translational (ribosome binding site (RBS)) and enzyme (mutant genes) variability to modulate pathway expression levels is essential for generating balanced metabolic pathways and maximizing the productivity of a strain. We report a novel, highly reliable and rapid single strand assembly (SSA) method for pathway engineering. The method was successfully optimized and applied to create constructs containing promoter, RBS and/or mutant enzyme libraries. To demonstrate its efficiency and reliability, the method was applied to fine-tune multi-gene pathways. Two promoter libraries were simultaneously introduced in front of two target genes, enabling orthogonal expression as demonstrated by principal component analysis. This shows that SSA will increase our ability to tune multi-gene pathways at all control levels for the biotechnological production of complex metabolites, achievable through the combinatorial modulation of transcription, translation and enzyme activity.

Analysis of DevR regulated genes in Mycobacterium tuberculosis

The DevRS two component system of Mycobacterium tuberculosis is responsible for its dormancy in host and becomes operative under hypoxic condition. It is experimentally known that phosphorylated DevR controls the expression of several downstream genes in a complex manner. In the present work we propose a theoretical model to show role of binding sites in DevR mediated gene expression. Individual and collective role of binding sites in regulating DevR mediated gene expression has been shown via modeling. Objective of the present work is twofold. First, to describe qualitatively the temporal dynamics of wild type genes and their known mutants. Based on these results we propose that DevR controlled gene expression follows a specific pattern which is efficient in describing other DevR mediated gene expression. Second, to analyze behavior of the system from information theoretical point of view. Using the tools of information theory we have calculated molecular efficiency of the system and have shown that it is close to the maximum limit of isothermal efficiency.

Transcript level and sequence determinants of protein abundance and noise in Escherichia coli

The range over which a protein is expressed, and its cell-to-cell variability, is often thought to be linked to the demand for its activity. Steady-state protein level is determined by multiple mechanisms controlling transcription and translation, many of which are limited by DNA- and RNA-encoded signals that affect initiation, elongation and termination of polymerases and ribosomes. We performed a comprehensive analysis of >100 sequence features to derive a predictive model composed of a minimal non-redundant set of factors explaining 66% of the total variation of protein abundance observed in >800 genes in Escherichia coli. The model suggests that protein abundance is primarily determined by the transcript level (53%) and by effectors of translation elongation (12%), whereas only a small fraction of the variation is explained by translational initiation (1%). Our analyses uncover a new sequence determinant, not previously described, affecting translation initiation and suggest that elongation rate is affected by both codon biases and specific amino acid composition. We also show that transcription and translation efficiency may have an effect on expression noise, which is more similar than previously assumed.

D-Tailor: automated analysis and design of DNA sequences

Motivation: Current advances in DNA synthesis, cloning and sequencing technologies afford high-throughput implementation of artificial sequences into living cells. However, flexible computational tools for multi-objective sequence design are lacking, limiting the potential of these technologies.

Results: We developed DNA-Tailor (D-Tailor), a fully extendable software framework, for property-based design of synthetic DNA sequences. D-Tailor permits the seamless integration of multiple sequence analysis tools into a generic Monte Carlo simulation that evolves sequences toward any combination of rationally defined properties. As proof of principle, we show that D-Tailor is capable of designing sequence libraries comprising all possible combinations among three different sequence properties influencing translation efficiency in Escherichia coli. The capacity to design artificial sequences that systematically sample any given parameter space should support the implementation of more rigorous experimental designs.

Availability: Source code is available for download at https://sourceforge.net/projects/dtailor/

Contact:aparkin@lbl.gov or cambray.guillaume@gmail.com

Supplementary information:Supplementary data are available at Bioinformatics online (D-Tailor Tutorial).

Enhanced production of coenzyme Q10 by self-regulating the engineered MEP pathway in Rhodobacter sphaeroides

Fine-tuning the expression level of an engineered pathway is crucial for the metabolic engineering of a host toward a desired phenotype. However, most engineered hosts suffer from nonfunctional protein expression, metabolic imbalance, cellular burden or toxicity from intermediates when an engineered pathway is first introduced, which can decrease production of the desired product. To circumvent these obstacles, we developed a self-regulation system utilizing the trc/tac promoter, LacI(q) protein and ribosomal binding sites (RBS). With the purpose of improving coenzyme Q10 (CoQ10 ) production by increasing the decaprenyl diphosphate supplement, enzymes DXS, DXR, IDI, and IspD were constitutively overexpressed under the control of the trc promoter in Rhodobacter sphaeroides. Then, a self-regulation system combining a set of RBSs for adjusting the expression of the LacI(q) protein was applied to tune the expression of the four genes, resulting in improved CoQ10 production. Finally, another copy of the tac promoter with the UbiG gene (involved in the ubiquinone pathway of CoQ10 biosynthesis) was introduced into the engineered pathway. By optimizing the expression level of both the upstream and downstream pathway, CoQ10 production in the mutants was improved up to 93.34 mg/L (7.16 mg/g DCW), about twofold of the wild-type (48.25 mg/L, 3.24 mg/g DCW).

Comprehensive identification of mutations responsible for heterogeneous vancomycin-intermediate Staphylococcus aureus (hVISA)-to-VISA conversion in laboratory-generated VISA strains derived from hVISA clinical strain Mu3

Heterogeneous vancomycin-intermediate Staphylococcus aureus (hVISA) spontaneously produces VISA cells within its cell population at a frequency of 10(-6) or greater. We established a total of 45 VISA mutant strains independently obtained from hVISA Mu3 and its related strains by one-step vancomycin selection. We then performed high-throughput whole-genome sequencing of the 45 strains and their parent strains to identify the genes involved in the hVISA-to-VISA phenotypic conversion. A comparative genome study showed that all the VISA strains tested carried a unique set of mutations. All of the 45 VISA strains carried 1 to 4 mutations possibly affecting the expression of a total of 48 genes. Among them, 32 VISA strains carried only one gene affected by a single mutation. As many as 20 genes in more than eight functional categories were affected in the 32 VISA strains, which explained the extremely high rates of the hVISA-to-VISA phenotypic conversion. Five genes, rpoB, rpoC, walK, pbp4, and pp2c, were previously reported as being involved in vancomycin resistance. Fifteen remaining genes were newly identified as associated with vancomycin resistance in this study. The gene most frequently affected (6 out of 32 strains) was cmk, which encodes cytidylate kinase, followed closely by rpoB (5 out of 32), encoding the β subunit of RNA polymerase. A mutation prevalence study also revealed a sizable number of cmk mutants among clinical VISA strains (7 out of 38 [18%]). Reduced cytidylate kinase activity in cmk mutant strains is proposed to contribute to the hVISA-to-VISA phenotype conversion by thickening the cell wall and reducing the cell growth rate.

Design and optimization of short DNA sequences that can be used as 5' fusion partners for high-level expression of heterologous genes in Escherichia coli

The 5' terminal nucleotide sequence of a gene is often a bottleneck in recombinant protein production. The ifn-α2bS gene is poorly expressed in Escherichia coli unless a translocation signal sequence (pelB) is fused to the 5' end of the gene. A combined in silico and in vivo analysis reported here further indicates that the ifn-α2bS 5' coding sequence is suboptimal for efficient gene expression. ifn-α2bS therefore presents a suitable model gene for describing properties of 5' fusions promoting expression. We show that short DNA sequences corresponding to the 5' end of the highly expressed celB gene, whose protein product is cytosolic, can functionally replace pelB as a 5' fusion partner for efficient ifn-α2bS expression. celB fusions of various lengths (corresponding to a minimum of 8 codons) led to more than 7- and 60-fold stimulation of expression at the transcript and protein levels, respectively. Moreover, the presence of a celB-based fusion partner was found to moderately reduce the decay rate of the corresponding transcript. The 5' fusions thus appear to act by enhancing translation, and bound ribosomes may accordingly contribute to increased mRNA stability and reduced mRNA decay. However, other effects, such as altered protein stability, cannot be excluded. We also developed an experimental protocol that enabled us to identify improved variants of the celB fusion, and one of these (celBD11) could be used to additionally increase ifn-α2bS expression more than 4-fold at the protein level. Interestingly, celBD11 also stimulated greater protein production of three other medically important human genes than the wild-type celB fragment.

Composability of regulatory sequences controlling transcription and translation in Escherichia coli

The inability to predict heterologous gene expression levels precisely hinders our ability to engineer biological systems. Using well-characterized regulatory elements offers a potential solution only if such elements behave predictably when combined. We synthesized 12,563 combinations of common promoters and ribosome binding sites and simultaneously measured DNA, RNA, and protein levels from the entire library. Using a simple model, we found that RNA and protein expression were within twofold of expected levels 80% and 64% of the time, respectively. The large dataset allowed quantitation of global effects, such as translation rate on mRNA stability and mRNA secondary structure on translation rate. However, the worst 5% of constructs deviated from prediction by 13-fold on average, which could hinder large-scale genetic engineering projects. The ease and scale this of approach indicates that rather than relying on prediction or standardization, we can screen synthetic libraries for desired behavior.

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.

Determination of the ribosome-binding sequence and spacer length between binding site and initiation codon for efficient protein expression in Bifidobacterium longum 105-A

The frequency of the ribosome-binding site (RBS) and spacer length between binding site and initiation codon were analyzed using informatics, and compared the activities of various synthetic RBSs and spacers in Bifidobacterium longum 105-A. As the result, AAGGAG and a 5 nt spacer length produced the most efficient protein expression.

Translation Initiation

Selection of correct start codons on messenger RNAs is a key step required for faithful translation of the genetic message. Such a selection occurs in a complex process, during which a translation-competent ribosome assembles, eventually having in its P site a specialized methionyl-tRNAMet base-paired with the start codon on the mRNA. This chapter summarizes recent advances describing at the molecular level the successive steps involved in the process. Special emphasis is put on the roles of the three initiation factors and of the initiator tRNA, which are crucial for the efficiency and the specificity of the process. In particular, structural analyses concerning complexes containing ribosomal subunits, as well as detailed kinetic studies, have shed new light on the sequence of events leading to faithful initiation of protein synthesis in Bacteria.

An efficient design strategy for a whole-cell biosensor based on engineered ribosome binding sequences

In prokaryotes, the ribosome binding sequence (RBS), located in the 5' untranslated region (5' UTR) of an mRNA, plays a critical role in enhancing mRNA translation and stability. To evaluate the effect of the RBS on the sensitivity and signal intensity of an environmental whole-cell biosensor, three Escherichia coli-based biosensors that respond to benzene, toluene, ethylbenzene, and the xylenes (BTEX) were constructed; the three biosensors have the same Pu promoter and xylR regulator from the Pseudomonas putida TOL plasmid but differ in the engineered RBS in their reporter genes. The results from time and dose-dependent induction of luminescence activity by 2-chlorotoluene showed that the BTEX-SE and BTEX-SD biosensors with engineered RBS had signal intensities approximately 10-35 times higher than the primary BTEX-W biosensor. The limits of detection (LOD) of the BTEX-SE and BTEX-SD biosensors were also significantly lower than the LOD of the BTEX-W biosensor (20 ± 5 μmol L(-1) and 25 ± 5 μmol L(-1) vs. 120 ± 10 μmol L(-1)). Moreover, the BTEX-SE and BTEX-SD biosensors responded three times more rapidly to the analytes. These results suggest that rationally designed RBS in the 5' UTR of a reporter gene may be a promising strategy for increasing the sensitivity, signal intensity, and response speed of whole-cell biosensors.

A brief review of molecular information theory

The idea that we could build molecular communications systems can be advanced by investigating how actual molecules from living organisms function. Information theory provides tools for such an investigation. This review describes how we can compute the average information in the DNA binding sites of any genetic control protein and how this can be extended to analyze its individual sites. A formula equivalent to Claude Shannon's channel capacity can be applied to molecular systems and used to compute the efficiency of protein binding. This efficiency is often 70% and a brief explanation for that is given. The results imply that biological systems have evolved to function at channel capacity, which means that we should be able to build molecular communications that are just as robust as our macroscopic ones.

Activity-based selection of HIV-1 reverse transcriptase variants with decreased polymerization fidelity

HIV-1 reverse transcriptase (HIV-1 RT) copies the RNA genome of HIV-1 into DNA, thereby committing errors at an exceptionally high frequency. Viral offspring evolve rapidly and consequently are capable of evading the immune response as well as antiviral treatment. However, error-prone viral replication could drive HIV close to extinction owing to an intolerable load of deleterious mutations. We applied a genetic selection scheme to identify variants of HIV-1 RT with a further increased error rate to study the relationship between error rate and viral replication. Using this approach, we identified 16 mutator candidates, two of which were purified and further studiedin vitro. One of these variant enzymes showed a generally increased mutation frequency as compared with the reference enzyme. A single amino acid residue, R448, is probably responsible for the observed effect. Mutation of this residue, which is located within the RNase H domain of HIV-1 RT, seems to perturb the interaction with template RNA and consequently affects polymerase activity and fidelity.