Multifactorial determinants of protein expression in prokaryotic open reading frames

Increasing recombinant protein production in E. coli via FACS-based selection of N-terminal coding DNA libraries

Here, we present a previously undescribed approach to modify N-terminal sequences of recombinant proteins to increase their production yield in Escherichia coli. Prior research has demonstrated that the nucleotides immediately following the start codon can significantly influence protein expression. However, the impact of these sequences is construct-specific and is not universally applicable to all proteins. Most of the previous research has been limited to selecting from a few rationally designed sequences. In contrast, we used a directed evolution-based methodology, screening large numbers of diversified sequences derived from DNA libraries coding for the N-termini of investigated proteins. To facilitate the identification of cells with increased expression of the target construct, we cloned a GFP gene at the C-terminus of the expressed genes and used fluorescent activated cell sorting (FACS) to separate cells based on their fluorescence. By following this systematic workflow, we successfully elevated the yield of soluble recombinant proteins of multiple constructs up to over 30-fold.

Accurate Identification of Periplasmic Urea-binding Proteins by Structure- and Genome Context-assisted Functional Analysis

ABC transporters are ancient and ubiquitous nutrient transport systems in bacteria and play a central role in defining lifestyles. Periplasmic solute-binding proteins (SBPs) are components that deliver ligands to their translocation machinery. SBPs have diversified to bind a wide range of ligands with high specificity and affinity. However, accurate assignment of cognate ligands remains a challenging problem in SBPs. Urea metabolism plays an important role in the nitrogen cycle; anthropogenic sources account for more than half of global nitrogen fertilizer. We report identification of urea-binding proteins within a large SBP sequence family that encodes diverse functions. By combining genetic linkage between SBPs, ABC transporter components, enzymes or transcription factors, we accurately identified cognate ligands, as we verified experimentally by biophysical characterization of ligand binding and crystallographic determination of the urea complex of a thermostable urea-binding homolog. Using three-dimensional structure information, these functional assignments were extrapolated to other members in the sequence family lacking genetic linkage information, which revealed that only a fraction bind urea. Using the same combined approaches, we also inferred that other family members bind various short-chain amides, aliphatic amino acids (leucine, isoleucine, valine), γ-aminobutyrate, and as yet unknown ligands. Comparative structural analysis revealed structural adaptations that encode diversification in these SBPs. Systematic assignment of ligands to SBP sequence families is key to understanding bacterial lifestyles, and also provides a rich source of biosensors for clinical and environmental analysis, such as the thermostable urea-binding protein identified here.

The expression of integron arrays is shaped by the translation rate of cassettes

Integrons are key elements in the rise and spread of multidrug resistance in Gram-negative bacteria. These genetic platforms capture cassettes containing promoterless genes and stockpile them in arrays of variable length. In the current integron model, expression of cassettes is granted by the Pc promoter in the platform and is assumed to decrease as a function of its distance. Here we explored this model using a large collection of 136 antibiotic resistance cassettes and show the effect of distance is in fact negligible. Instead, cassettes have a strong impact in the expression of downstream genes because their translation rate affects the stability of the whole polycistronic mRNA molecule. Hence, cassettes with reduced translation rates decrease the expression and resistance phenotype of cassettes downstream. Our data puts forward an integron model in which expression is contingent on the translation of cassettes upstream, rather than on the distance to the Pc.

Selection on synonymous sites: the unwanted transcript hypothesis

Although translational selection to favour codons that match the most abundant tRNAs is not readily observed in humans, there is nonetheless selection in humans on synonymous mutations. We hypothesize that much of this synonymous site selection can be explained in terms of protection against unwanted RNAs - spurious transcripts, mis-spliced forms or RNAs derived from transposable elements or viruses. We propose not only that selection on synonymous sites functions to reduce the rate of creation of unwanted transcripts (for example, through selection on exonic splice enhancers and cryptic splice sites) but also that high-GC content (but low-CpG content), together with intron presence and position, is both particular to functional native mRNAs and used to recognize transcripts as native. In support of this hypothesis, transcription, nuclear export, liquid phase condensation and RNA degradation have all recently been shown to promote GC-rich transcripts and suppress AU/CpG-rich ones. With such 'traps' being set against AU/CpG-rich transcripts, the codon usage of native genes has, in turn, evolved to avoid such suppression. That parallel filters against AU/CpG-rich transcripts also affect the endosomal import of RNAs further supports the unwanted transcript hypothesis of synonymous site selection and explains the similar design rules that have enabled the successful use of transgenes and RNA vaccines.

Nonspecific N-terminal tetrapeptide insertions disrupt the translation arrest induced by ribosome-arresting peptide sequences

The nascent polypeptide chains passing through the ribosome tunnel not only serve as an intermediate of protein synthesis but also, in some cases, act as dynamic genetic information, controlling translation through interaction with the ribosome. One notable example is Escherichia coli SecM, in which translation of the ribosome arresting peptide (RAP) sequence in SecM leads to robust elongation arrest. Translation regulations, including the SecM-induced translation arrest, play regulatory roles such as gene expression control. Recent investigations have indicated that the insertion of a peptide sequence, SKIK (or MSKIK), into the adjacent N-terminus of the RAP sequence of SecM behaves as an "arrest canceler". As the study did not provide a direct assessment of the strength of translation arrest, we conducted detailed biochemical analyses. The results revealed that the effect of SKIK insertion on weakening SecM-induced translation arrest was not specific to the SKIK sequence, that is, other tetrapeptide sequences inserted just before the RAP sequence also attenuated the arrest. Our data suggest that SKIK or other tetrapeptide insertions disrupt the context of the RAP sequence rather than canceling or preventing the translation arrest.

A new look at the fluorescent protein-based approach for identifying optimal coding sequence for recombinant protein expression in E. coli

Due to the degeneracy of the genetic code, most amino acids are encoded by several codons. The choice among synonymous codons at the N-terminus of genes has a profound effect on protein expression in Escherichia coli. This is often explained by the different contributions of synonymous codons to mRNA secondary structure formation. Strong secondary structures at the 5'-end of mRNA interfere with ribosome binding and affect the process of translation initiation. In silico optimization of the gene 5'-end can significantly increase the level of protein expression; however, this method is not always effective due to the uncertainty of the exact mechanism by which synonymous substitutions affect expression; thus, it may produce nonoptimal variants as well as miss some of the best producers. In this paper, an alternative approach is proposed based on screening a partially randomized library of expression constructs comprising hundreds of selected synonymous variants. The effect of such substitutions was evaluated using the gene of interest fused to the reporter gene of the fluorescent protein with subsequent screening for the most promising candidates according to the reporter's signal intensity. The power of the approach is demonstrated by a significant increase in the prokaryotic expression of three proteins: canine cystatin C, human BCL2-associated athanogene 3 and human cardiac troponin I. This simple approach was suggested which may provide an efficient, easy, and inexpensive optimization method for poorly expressed proteins in bacteria.

Genes for highly abundant proteins in Escherichia coli avoid 5' codons that promote ribosomal initiation

In many species highly expressed genes (HEGs) over-employ the synonymous codons that match the more abundant iso-acceptor tRNAs. Bacterial transgene codon randomization experiments report, however, that enrichment with such "translationally optimal" codons has little to no effect on the resultant protein level. By contrast, consistent with the view that ribosomal initiation is rate limiting, synonymous codon usage following the 5' ATG greatly influences protein levels, at least in part by modifying RNA stability. For the design of bacterial transgenes, for simple codon based in silico inference of protein levels and for understanding selection on synonymous mutations, it would be valuable to computationally determine initiation optimality (IO) scores for codons for any given species. One attractive approach is to characterize the 5' codon enrichment of HEGs compared with the most lowly expressed genes, just as translational optimality scores of codons have been similarly defined employing the full gene body. Here we determine the viability of this approach employing a unique opportunity: for Escherichia coli there is both the most extensive protein abundance data for native genes and a unique large-scale transgene codon randomization experiment enabling objective definition of the 5' codons that cause, rather than just correlate with, high protein abundance (that we equate with initiation optimality, broadly defined). Surprisingly, the 5' ends of native genes that specify highly abundant proteins avoid such initiation optimal codons. We find that this is probably owing to conflicting selection pressures particular to native HEGs, including selection favouring low initiation rates, this potentially enabling high efficiency of ribosomal usage and low noise. While the classical HEG enrichment approach does not work, rendering simple prediction of native protein abundance from 5' codon content futile, we report evidence that initiation optimality scores derived from the transgene experiment may hold relevance for in silico transgene design for a broad spectrum of bacteria.

Large-scale analysis of mRNA sequences localized near the start and amber codons and their impact on the diversity of mRNA display libraries

Extremely diverse libraries are essential for effectively selecting functional peptides or proteins, and mRNA display technology is a powerful tool for generating such libraries with over 1012-1013 diversity. Particularly, the protein-puromycin linker (PuL)/mRNA complex formation yield is determining for preparing the libraries. However, how mRNA sequences affect the complex formation yield remains unclear. To study the effects of N-terminal and C-terminal coding sequences on the complex formation yield, puromycin-attached mRNAs containing three random codons after the start codon (32768 sequences) or seven random bases next to the amber codon (6480 sequences) were translated. Enrichment scores were calculated by dividing the appearance rate of every sequence in protein-PuL/mRNA complexes by that in total mRNAs. The wide range of enrichment scores (0.09-2.10 for N-terminal and 0.30-4.23 for C-terminal coding sequences) indicated that the N-terminal and C-terminal coding sequences strongly affected the complex formation yield. Using C-terminal GGC-CGA-UAG-U sequences, which resulted in the highest enrichment scores, we constructed highly diverse libraries of monobodies and macrocyclic peptides. The present study provides insights into how mRNA sequences affect the protein/mRNA complex formation yield and will accelerate the identification of functional peptides and proteins involved in various biological processes and having therapeutic applications.

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.

Flow-Seq Method: Features and Application in Bacterial Translation Studies

The Flow-seq method is based on using reporter construct libraries, where a certain element regulating the gene expression of fluorescent reporter proteins is represented in many thousands of variants. Reporter construct libraries are introduced into cells, sorted according to their fluorescence level, and then subjected to next-generation sequencing. Therefore, it turns out to be possible to identify patterns that determine the expression efficiency, based on tens and hundreds of thousands of reporter constructs in one experiment. This method has become common in evaluating the efficiency of protein synthesis simultaneously by multiple mRNA variants. However, its potential is not confined to this area. In the presented review, a comparative analysis of the Flow-seq method and other alternative approaches used for translation efficiency evaluation of mRNA was carried out; the features of its application and the results obtained by Flow-seq were also considered.

Unexpected myriad of co-occurring viral strains and species in one of the most abundant and microdiverse viruses on Earth

Viral genetic microdiversity drives adaptation, pathogenicity, and speciation and has critical consequences for the viral-host arms race occurring at the strain and species levels, which ultimately impact microbial community structure and biogeochemical cycles. Despite the fact that most efforts have focused on viral macrodiversity, little is known about the microdiversity of ecologically important viruses on Earth. Recently, single-virus genomics discovered the putatively most abundant ocean virus in temperate and tropical waters: the uncultured dsDNA virus vSAG 37-F6 infecting Pelagibacter, the most abundant marine bacteria. In this study, we report the cooccurrence of up to ≈1,500 different viral strains (>95% nucleotide identity) and ≈30 related species (80-95% nucleotide identity) in a single oceanic sample. Viral microdiversity was maintained over space and time, and most alleles were the result of synonymous mutations without any apparent adaptive benefits to cope with host translation codon bias and efficiency. Gene flow analysis used to delimitate species according to the biological species concept (BSC) revealed the impact of recombination in shaping vSAG 37-F6 virus and Pelagibacter speciation. Data demonstrated that this large viral microdiversity somehow mirrors the host species diversity since ≈50% of the 926 analyzed Pelagibacter genomes were found to belong to independent BSC species that do not significantly engage in gene flow with one another. The host range of this evolutionarily successful virus revealed that a single viral species can infect multiple Pelagibacter BSC species, indicating that this virus crosses not only formal BSC barriers but also biomes since viral ancestors are found in freshwater.

Discovery of Thermostable, Fluorescently Responsive Glucose Biosensors by Structure-Assisted Function Extrapolation

Accurate assignment of protein function from sequence remains a fascinating and difficult challenge. The periplasmic-binding protein (PBP) superfamily present an interesting case of function prediction because they are both ubiquitous in prokaryotes and tend to diversify through gene duplication "explosions" that can lead to large numbers of paralogs in a genome. An engineered version of the moderately thermostable glucose-binding PBP from Escherichia coli has been used successfully as a reagentless fluorescent biosensor both in vitro and in vivo. To develop more robust sensors that meet the challenges of real-world applications, we report the discovery of thermostable homologues that retain a glucose-mediated conformationally coupled fluorescence response. Accurately identifying a glucose-binding PBP homologue among closely related paralogs is challenging. We demonstrate that a structure-based method that filters sequences by residues that bind glucose in an archetype structure is highly effective. Using fully sequenced bacterial genomes, we found that this filter reduced high paralog numbers to single hits in a genome, consistent with the accurate separation of glucose binding from other functions. We expressed engineered proteins for eight homologues, chosen to represent different degrees of sequence identity, and tested their glucose-mediated fluorescence responses. We accurately predicted the presence of glucose binding down to 31% sequence identity. We have also successfully identified suitable candidates for next-generation robust, fluorescent glucose sensors.

Nascent polypeptide within the exit tunnel stabilizes the ribosome to counteract risky translation

Continuous translation elongation, irrespective of amino acid sequences, is a prerequisite for living organisms to produce their proteomes. However, nascent polypeptide products bear an inherent risk of elongation abortion. For example, negatively charged sequences with occasional intermittent prolines, termed intrinsic ribosome destabilization (IRD) sequences, weaken the translating ribosomal complex, causing certain nascent chain sequences to prematurely terminate translation. Here, we show that most potential IRD sequences in the middle of open reading frames remain cryptic and do not interrupt translation, due to two features of the nascent polypeptide. Firstly, the nascent polypeptide itself spans the exit tunnel, and secondly, its bulky amino acid residues occupy the tunnel entrance region, thereby serving as a bridge and protecting the large and small ribosomal subunits from dissociation. Thus, nascent polypeptide products have an inbuilt ability to ensure elongation continuity.

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.

A 68-codon genetic code to incorporate four distinct non-canonical amino acids enabled by automated orthogonal mRNA design

Orthogonal (O) ribosome-mediated translation of O-mRNAs enables the incorporation of up to three distinct non-canonical amino acids (ncAAs) into proteins in Escherichia coli (E. coli). However, the general and efficient incorporation of multiple distinct ncAAs by O-ribosomes requires scalable strategies for both creating efficiently and specifically translated O-mRNAs, and the compact expression of multiple O-aminoacyl-tRNA synthetase (O-aaRS)/O-tRNA pairs. We automate the discovery of O-mRNAs that lead to up to 40 times more protein, and are up to 50-fold more orthogonal, than previous O-mRNAs; protein yields from our O-mRNAs match or exceed those from wild-type mRNAs. These advances enable a 33-fold increase in yield for incorporating three distinct ncAAs. We automate the creation of operons for O-tRNA genes, and develop operons for O-aaRS genes. Combining our advances creates a 68-codon, 24-amino-acid genetic code to efficiently incorporate four distinct ncAAs into a single protein in response to four distinct quadruplet codons.

Inferring Adaptive Codon Preference to Understand Sources of Selection Shaping Codon Usage Bias

Alternative synonymous codons are often used at unequal frequencies. Classically, studies of such codon usage bias (CUB) attempted to separate the impact of neutral from selective forces by assuming that deviations from a predicted neutral equilibrium capture selection. However, GC-biased gene conversion (gBGC) can also cause deviation from a neutral null. Alternatively, selection has been inferred from CUB in highly expressed genes, but the accuracy of this approach has not been extensively tested, and gBGC can interfere with such extrapolations (e.g., if expression and gene conversion rates covary). It is therefore critical to examine deviations from a mutational null in a species with no gBGC. To achieve this goal, we implement such an analysis in the highly AT rich genome of Dictyostelium discoideum, where we find no evidence of gBGC. We infer neutral CUB under mutational equilibrium to quantify "adaptive codon preference," a nontautologous genome wide quantitative measure of the relative selection strength driving CUB. We observe signatures of purifying selection consistent with selection favoring adaptive codon preference. Preferred codons are not GC rich, underscoring the independence from gBGC. Expression-associated "preference" largely matches adaptive codon preference but does not wholly capture the influence of selection shaping patterns across all genes, suggesting selective constraints associated specifically with high expression. We observe patterns consistent with effects on mRNA translation and stability shaping adaptive codon preference. Thus, our approach to quantifying adaptive codon preference provides a framework for inferring the sources of selection that shape CUB across different contexts within the genome.

Rational Design of the N-Terminal Coding Sequence for Regulating Enzyme Expression in Bacillus subtilis

Synonymous mutation of the N-terminal coding sequence (NCS) has been used to regulate gene expression. We here developed a statistical model to predict the effect of the NCSs on protein expression in Bacillus subtilis WB600. First, a synonymous mutation was performed within the first 10 residues of a superfolder green fluorescent protein to generate a library of 172 NCS synonymous mutants with different expression levels. A prediction model was then developed, which adopted G/C frequency at the third position of each codon and minimum free energy of mRNA as the independent variables, using multiple regression analysis between the 11 sequence parameters of the NCS and their fluorescence intensities. By designing the NCS of the 10 signal peptides de novo according to the model, the extracellular yield of B. subtilis pullulanase fused to each signal peptide was up-regulated by up to 515% or down-regulated by at most 79%. This work provided a candidate tool for fine-tuning gene expression or enzyme production in B. subtilis.

Translational Control by Ribosome Pausing in Bacteria: How a Non-uniform Pace of Translation Affects Protein Production and Folding

Protein homeostasis of bacterial cells is maintained by coordinated processes of protein production, folding, and degradation. Translational efficiency of a given mRNA depends on how often the ribosomes initiate synthesis of a new polypeptide and how quickly they read the coding sequence to produce a full-length protein. The pace of ribosomes along the mRNA is not uniform: periods of rapid synthesis are separated by pauses. Here, we summarize recent evidence on how ribosome pausing affects translational efficiency and protein folding. We discuss the factors that slow down translation elongation and affect the quality of the newly synthesized protein. Ribosome pausing emerges as important factor contributing to the regulatory programs that ensure the quality of the proteome and integrate the cellular and environmental cues into regulatory circuits of the cell.

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.

Improving the Expression of Human Granulocyte Colony Stimulating Factor in Escherichia coli by Reducing the GC-content and Increasing mRNA Folding Free Energy at 5'-Terminal End

Purpose: Strategy for improving the production of biopharmaceutical protein continues to develop due to increasing market demand. Human granulocyte colony stimulating factor (hG-CSF) is one of biopharmaceutical proteins that has many applications, and easily produced in Escherichia coli expression system. Previous studies reported that codon usage, rare codon, mRNA folding and GC-content at 5’-terminal end were crucial for protein production in E. coli. In the present study, the effect of reducing the GC-content and increasing the mRNA folding free energy at the 5’-terminal end on the expression level of hG-CSF proteins was investigated.

Methods: Synonymous codon substitutions were performed to generate mutant variants of open reading frame (ORF) with lower GC-content at 5’-terminal ends. Oligoanalyzer tool was used to calculate the GC content of eight codons sequence after ATG. Whereas, mRNA folding free energy was predicted using KineFold and RNAfold tools. The template DNA was amplified using three variant forward primers and one same reverse primer. Those DNA fragments were individually cloned into pJexpress414 expression vector and were confirmed using restriction and DNA sequencing analyses. The confirmed constructs were transformed into E. coli NiCo21(DE3) host cells and the recombinant protein was expressed using IPTG-induction. Total protein obtained were characterized using SDS-PAGE, Western blot and ImageJ software analyses.

Results: The result showed that the mutant variant with lower GC-content and higher mRNA folding free energy near the translation initiation region (TIR) could produce a higher amount of hG-CSF proteins compared to the original gene sequence.

Conclusion: This study emphasized the important role of the nucleotide composition immediately downstream the start codon to achieve high-yield protein product on heterologous expression in E. coli.

Translation at first sight: the influence of leading codons

First triplets of mRNA coding region affect the yield of translation. We have applied the flowseq method to analyze >30 000 variants of the codons 2-11 of the fluorescent protein reporter to identify factors affecting the protein synthesis. While the negative influence of mRNA secondary structure on translation has been confirmed, a positive role of rare codons at the beginning of a coding sequence for gene expression has not been observed. The identity of triplets proximal to the start codon contributes more to the protein yield then more distant ones. Additional in-frame start codons enhance translation, while Shine-Dalgarno-like motifs downstream the initiation codon are inhibitory. The metabolic cost of amino acids affects the yield of protein in the poor medium. The most efficient translation was observed for variants with features resembling those of native Escherichia coli genes.

Determinants of translation efficiency in the evolutionarily-divergent protist Trichomonas vaginalis

BACKGROUND

Trichomonas vaginalis, the causative agent of a prevalent urogenital infection in humans, is an evolutionarily divergent protozoan. Protein-coding genes in T. vaginalis are largely controlled by two core promoter elements, producing mRNAs with short 5' UTRs. The specific mechanisms adopted by T. vaginalis to fine-tune the translation efficiency (TE) of mRNAs remain largely unknown.

RESULTS

Using both computational and experimental approaches, this study investigated two key factors influencing TE in T. vaginalis: codon usage and mRNA secondary structure. Statistical dependence between TE and codon adaptation index (CAI) highlighted the impact of codon usage on mRNA translation in T. vaginalis. A genome-wide interrogation revealed that low structural complexity at the 5' end of mRNA followed closely by a highly structured downstream region correlates with TE variation in this organism. To validate these findings, a synthetic library of 15 synonymous iLOV genes was created, representing five mRNA folding profiles and three codon usage profiles. Fluorescence signals produced by the expression of these synonymous iLOV genes in T. vaginalis were consistent with and validated our in silico predictions.

CONCLUSIONS

This study demonstrates the role of codon usage bias and mRNA secondary structure in TE of T. vaginalis mRNAs, contributing to a better understanding of the factors that influence, and possibly regulate, gene expression in this human pathogen.

Widespread non-modular overlapping codes in the coding regions

Messenger RNAs (mRNAs) consist of a coding region (open reading frame (ORF)) and two untranslated regions (UTRs), 5'UTR and 3'UTR. Ribosomes travel along the coding region, translating nucleotide triplets (called codons) to a chain of amino acids. The coding region was long believed to mainly encode the amino acid content of proteins, whereas regulatory signals reside in the UTRs and in other genomic regions. However, in recent years we have learned that the ORF is expansively populated with various regulatory signals, or codes, which are related to all gene expression steps and additional intracellular aspects. In this paper, we review the current knowledge related to overlapping codes inside the coding regions, such as the influence of synonymous codon usage on translation speed (and, in turn, the effect of translation speed on protein folding), ribosomal frameshifting, mRNA stability, methylation, splicing, transcription and more. All these codes come together and overlap in the ORF sequence, ensuring production of the right protein at the right time.

A short translational ramp determines the efficiency of protein synthesis

Translation initiation is a major rate-limiting step for protein synthesis. However, recent studies strongly suggest that the efficiency of protein synthesis is additionally regulated by multiple factors that impact the elongation phase. To assess the influence of early elongation on protein synthesis, we employed a library of more than 250,000 reporters combined with in vitro and in vivo protein expression assays. Here we report that the identity of the amino acids encoded by codons 3 to 5 impact protein yield. This effect is independent of tRNA abundance, translation initiation efficiency, or overall mRNA structure. Single-molecule measurements of translation kinetics revealed pausing of the ribosome and aborted protein synthesis on codons 4 and 5 of distinct amino acid and nucleotide compositions. Finally, introduction of preferred sequence motifs only at specific codon positions improves protein synthesis efficiency for recombinant proteins. Collectively, our data underscore the critical role of early elongation events in translational control of gene expression.

Several factors contribute to the efficiency of protein expression. Here the authors show that the identity of amino acids encoded by codons at position 3–5 significantly impact translation efficiency and protein expression levels.

Constructive approach for synthesis of a functional IgG using a reconstituted cell-free protein synthesis system

IgG is an indispensable biological experimental tool as well as a widely-used therapeutic protein. However, cell culture-based expression of monoclonal IgG is costly and time-consuming, making this process difficult to use for high-throughput screening in early-stage evaluation of biologics. With the goal of establishing a fast, simple, and robust high-throughput expression system for IgG, we implemented the synthesis of functional aglycosylated IgG by constructive approach based on a reconstituted prokaryotic cell-free protein synthesis system (PURE system). Optimization of the PURE system revealed that the following factors and reaction conditions were needed for IgG synthesis: (1) inclusion of the disulfide bond isomerase DsbC, (2) adjustment of the GSH/GSSG ratio, (3) inclusion of the molecular chaperone DnaK and its cofactors, and (4) use of an extended incubation time. Synthesis temperature and template DNA ratio (light chain-/heavy chain-encoding) also had been optimized for each IgG. Under optimal conditions, peak production of the anti-HER2 antibody trastuzumab reached 124 µg/mL. Furthermore, the active forms of other IgGs, including IgG1, IgG2, and IgG4 subclasses, also were synthesized. These results provide basic information for the development of novel high-throughput expression and functional screening systems for IgG, as well as useful information for understanding the IgG synthesis process.

Transcription and translation of the sigG gene is tuned for proper execution of the switch from early to late gene expression in the developing Bacillus subtilis spore

A cascade of alternative sigma factors directs developmental gene expression during spore formation by the bacterium Bacillus subtilis. As the spore develops, a tightly regulated switch occurs in which the early-acting sigma factor σ^F is replaced by the late-acting sigma factor σ^G. The gene encoding σ^G (sigG) is transcribed by σ^F and by σ^G itself in an autoregulatory loop; yet σ^G activity is not detected until σ^F-dependent gene expression is complete. This separation in σ^F and σ^G activities has been suggested to be due at least in part to a poorly understood intercellular checkpoint pathway that delays sigG expression by σ^F. Here we report the results of a careful examination of sigG expression during sporulation. Unexpectedly, our findings argue against the existence of a regulatory mechanism to delay sigG transcription by σ^F and instead support a model in which sigG is transcribed by σ^F with normal timing, but at levels that are very low. This low-level expression of sigG is the consequence of several intrinsic features of the sigG regulatory and coding sequence—promoter spacing, secondary structure potential of the mRNA, and start codon identity—that dampen its transcription and translation. Especially notable is the presence of a conserved hairpin in the 5’ leader sequence of the sigG mRNA that occludes the ribosome-binding site, reducing translation by up to 4-fold. Finally, we demonstrate that misexpression of sigG from regulatory and coding sequences lacking these features triggers premature σ^G activity in the forespore during sporulation, as well as inappropriate σ^G activity during vegetative growth. Altogether, these data indicate that transcription and translation of the sigG gene is tuned to prevent vegetative expression of σ^G and to ensure the precise timing of the switch from σ^F to σ^G in the developing spore.

Global changes in gene expression occur during normal cellular growth and development, as well as during cancer cell transformation and bacterial pathogenesis. In this study we have investigated the molecular mechanisms that drive the switch from early to late developmental gene expression during spore formation by the model bacterium Bacillus subtilis. At early times, gene expression in the developing spore is directed by the transcription factor σ^F; at later times σ^F is replaced by σ^G. An important, yet poorly understood aspect of this σ^F-to-σ^G transition is how σ^G activation is delayed until the early, σ^F-directed phase of gene expression is complete. Here we have carefully examined expression of the gene encoding σ^G, sigG, and found that its transcription and translation are ordinarily dampened by several features of its regulatory and coding sequences. Moreover, we have found that this “tuning” of sigG expression is required for proper timing of the switch to σ^G. These results reframe our understanding of how sigG is regulated during B. subtilis sporulation and, more broadly, advance our understanding of how global changes in gene expression can be precisely executed at the molecular/genetic level.

Pathway Design, Engineering, and Optimization

The microbial metabolic versatility found in nature has inspired scientists to create microorganisms capable of producing value-added compounds. Many endeavors have been made to transfer and/or combine pathways, existing or even engineered enzymes with new function to tractable microorganisms to generate new metabolic routes for drug, biofuel, and specialty chemical production. However, the success of these pathways can be impeded by different complications from an inherent failure of the pathway to cell perturbations. Pursuing ways to overcome these shortcomings, a wide variety of strategies have been developed. This chapter will review the computational algorithms and experimental tools used to design efficient metabolic routes, and construct and optimize biochemical pathways to produce chemicals of high interest.

In vitro gene expression-coupled bacterial cell chip for screening species-specific antimicrobial enzymes

Targeting infectious bacterial pathogens is important for reducing the evolution of antibiotic-resistant bacteria and preserving the endogenous human microbiome. Cell lytic enzymes including bacteriophage endolysins, bacterial autolysins, and other bacteriolysins are useful antibiotic alternatives due to their exceptional target selectivity, which may be used to lysins rapidly kill target bacteria and their high specificity permit the normal commensal microflora to be left undisturbed. Genetic information of numerous lysins is currently available, but the identification of their antimicrobial function and specificity has been limited because most lysins are often poorly expressed and exhibit low solubilities. Here, we report the development of bacterial cell chip for rapidly accessing the function of diverse genes that are suggestive of encoding lysins. This approach can be used to evaluate rapidly the species-specific antimicrobial activity of diverse lysins synthesized from in vitro transcription and translation (TNT) of plasmid DNA. In addition, new potent lysins can be assessed that are not expressed in hosts and display low solubility. As a result of evaluating the species-specific antimicrobial function of 11 (un)known lysins with an in vitro TNT-coupled bacterial cell chip, a potent recombinant lysin against Staphylococcus strains, SA1, was identified. The SA1 was highly potent against not only S. aureus, but also both lysostaphin-resistant S. simulans and S. epidermidis cells. To this end, the SA1 may be applicable to treat both methicillin-resistant S. aureus (MRSA) and lysostaphin-resistant MRSA mutants. Biotechnol. Bioeng. 2017;114: 1648-1657. © 2017 Wiley Periodicals, Inc.

A bicistronic vector with destabilized mRNA secondary structure yields scalable higher titer expression of human neurturin in E. coli

Human neurturin (NTN) is a cystine knot growth factor with potential therapeutic use in diseases such as Parkinson's and diabetes. Scalable high titer production of native NTN is particularly challenging because of the cystine knot structure which consists of an embedded ring comprised of at least three disulfide bonds. We sought to pursue enhanced scalable production of NTN in Escherichia coli. Our initial efforts focused on codon optimization of the first two codons following AUG, but these studies resulted in only a marginal increase in NTN expression. Therefore, we pursued an alternative strategy of using a bicistronic vector for NTN expression designed to reduce mRNA secondary structure to achieve increased ribosome binding and re-initiation. The first cistron was designed to prevent sequestration of the translation initiation region in a secondary conformation. The second cistron, which contained the NTN coding sequence itself, was engineered to disrupt double bonded base pairs and destabilize the secondary structure for ribosome re-initiation. The ensemble approach of reducing NTN's mRNA secondary structure and using the bicistronic vector had an additive effect resulting in significantly increased NTN expression. Our strain selection studies were conducted in a miniaturized bioreactor. An optimized strain was selected and scaled up to a 100 L fermentor, which yielded an inclusion body titer of 2 g/L. The inclusion bodies were refolded to yield active NTN. We believe that our strategy is applicable to other candidate proteins that are difficult-to-express due to stable mRNA secondary structures. Biotechnol. Bioeng. 2017;114: 1753-1761. © 2017 Wiley Periodicals, Inc.

The Frequency of Internal Shine-Dalgarno-like Motifs in Prokaryotes

In prokaryotes, translation initiation typically depends on complementary binding between a G-rich Shine–Dalgarno (SD) motif in the 5′ untranslated region of mRNAs, and the 3′ tail of the 16S ribosomal RNA (the anti-SD sequence). In some cases, internal SD-like motifs in the coding region generate “programmed” ribosomal pauses that are beneficial for protein folding or accurate targeting. On the other hand, such pauses can also reduce protein production, generating purifying selection against internal SD-like motifs. This selection should be stronger in GC-rich genomes that are more likely to harbor the G-rich SD motif. However, the nature and consequences of selection acting on internal SD-like motifs within genomes and across species remains unclear. We analyzed the frequency of SD-like hexamers in the coding regions of 284 prokaryotes (277 with known anti-SD sequences and 7 without a typical SD mechanism). After accounting for GC content, we found that internal SD-like hexamers are avoided in 230 species, including three without a typical SD mechanism. The degree of avoidance was higher in GC-rich genomes, mesophiles, and N-terminal regions of genes. In contrast, 54 species either showed no signature of avoidance or were enriched in internal SD-like motifs. C-terminal gene regions were relatively enriched in SD-like hexamers, particularly for genes in operons or those followed closely by downstream genes. Together, our results suggest that the frequency of internal SD-like hexamers is governed by multiple factors including GC content and genome organization, and further empirical work is necessary to understand the evolution and functional roles of these motifs.

The ribosome in action: Tuning of translational efficiency and protein folding

The cellular proteome is shaped by the combined activities of the gene expression and quality control machineries. While transcription plays an undoubtedly important role, in recent years also translation emerged as a key step that defines the composition and quality of the proteome and the functional activity of proteins in the cell. Among the different post-transcriptional control mechanisms, translation initiation and elongation provide multiple checkpoints that can affect translational efficiency. A multitude of specific signals in mRNAs can determine the frequency of translation initiation, choice of the open reading frame, global and local elongation velocities, and the folding of the emerging protein. In addition to specific signatures in the mRNAs, also variations in the global pools of translation components, including ribosomes, tRNAs, mRNAs, and translation factors can alter translational efficiencies. The cellular outcomes of phenomena such as mRNA codon bias are sometimes difficult to understand due to the staggering complexity of covariates that affect codon usage, translation, and protein folding. Here we summarize the experimental evidence on how the ribosome-together with the other components of the translational machinery-can alter translational efficiencies of mRNA at the initiation and elongation stages and how translation velocity affects protein folding. We seek to explain these findings in the context of mechanistic work on the ribosome. The results argue in favour of a new understanding of translation control as a hub that links mRNA homeostasis to production and quality control of proteins in the cell.

Cooperation between the chloroplast psbA 5'-untranslated region and coding region is important for translational initiation: the chloroplast translation machinery cannot read a human viral gene coding region

Chloroplast mRNA translation is regulated by the 5'-untranslated region (5'-UTR). Chloroplast 5'-UTRs also support translation of the coding regions of heterologous genes. Using an in vitro translation system from tobacco chloroplasts, we detected no translation from a human immunodeficiency virus tat coding region fused directly to the tobacco chloroplast psbA 5'-UTR. This lack of apparent translation could have been due to rapid degradation of mRNA templates or synthesized protein products. Replacing the psbA 5'-UTR with the E. coli phage T7 gene 10 5'-UTR, a highly active 5'-UTR, and substituting synonymous codons led to some translation of the tat coding region. The Tat protein thus synthesized was stable during translation reactions. No significant degradation of the added tat mRNAs was observed after translation reactions. These results excluded the above two possibilities and confirmed that the tat coding region prevented its own translation. The tat coding region was then fused to the psbA 5'-UTR with a cognate 5'-coding segment. Significant translation was detected from the tat coding region when fused after 10 or more codons. That is, translation could be initiated from the tat coding region once translation had started, indicating that the tat coding region inhibits translational initiation but not elongation. Hence, cooperation/compatibility between the 5'-UTR and its coding region is important for translational initiation.

Large-Effect Beneficial Synonymous Mutations Mediate Rapid and Parallel Adaptation in a Bacterium

Contrary to previous understanding, recent evidence indicates that synonymous codon changes may sometimes face strong selection. However, it remains difficult to generalize the nature, strength, and mechanism(s) of such selection. Previously, we showed that synonymous variants of a key enzyme-coding gene (fae) of Methylobacterium extorquens AM1 decreased enzyme production and reduced fitness dramatically. We now show that during laboratory evolution, these variants rapidly regained fitness via parallel yet variant-specific, highly beneficial point mutations in the N-terminal region of fae. These mutations (including four synonymous mutations) had weak but consistently positive impacts on transcript levels, enzyme production, or enzyme activity. However, none of the proposed mechanisms (including internal ribosome pause sites or mRNA structure) predicted the fitness impact of evolved or additional, engineered point mutations. This study shows that synonymous mutations can be fixed through strong positive selection, but the mechanism for their benefit varies depending on the local sequence context.

Enhanced Protein Production in Escherichia coli by Optimization of Cloning Scars at the Vector-Coding Sequence Junction

Protein production in Escherichia coli is a fundamental activity for a large fraction of academic, pharmaceutical, and industrial research laboratories. Maximum production is usually sought, as this reduces costs and facilitates downstream purification steps. Frustratingly, many coding sequences are poorly expressed even when they are codon-optimized and expressed from vectors with powerful genetic elements. In this study, we show that poor expression can be caused by certain nucleotide sequences (e.g., cloning scars) at the junction between the vector and the coding sequence. Since these sequences lie between the Shine-Dalgarno sequence and the start codon, they are an integral part of the translation initiation region. To identify the most optimal sequences, we devised a simple and inexpensive PCR-based step that generates sequence variants at the vector-coding sequence junction. These sequence variants modulated expression by up to 1000-fold. FACS-seq analyses indicated that low GC content and relaxed mRNA stability (ΔG) in this region were important, but not the only, determinants for high expression.

How do bacteria tune translation efficiency?

Bacterial proteins are translated with precisely determined rates to meet cellular demand. In contrast, efforts to express recombinant proteins in bacteria are often met with large unpredictability in their levels of translation. The disconnect between translation of natural and synthetic mRNA stems from the lack of understanding of the strategy used by bacteria to tune translation efficiency (TE). The development of array-based oligonucleotide synthesis and ribosome profiling provides new approaches to address this issue. Although the major determinant for TE is still unknown, these high-throughput studies point out a statistically significant but mild contribution from the mRNA secondary structure around the start codon. Here I summarize those findings and provide a theoretical framework for measuring TE.

Rationally designed, heterologous S. cerevisiae transcripts expose novel expression determinants

Deducing generic causal relations between RNA transcript features and protein expression profiles from endogenous gene expression data remains a major unsolved problem in biology. The analysis of gene expression from heterologous genes contributes significantly to solving this problem, but has been heavily biased toward the study of the effect of 5' transcript regions and to prokaryotes. Here, we employ a synthetic biology driven approach that systematically differentiates the effect of different regions of the transcript on gene expression up to 240 nucleotides into the ORF. This enabled us to discover new causal effects between features in previously unexplored regions of transcripts, and gene expression in natural regimes. We rationally designed, constructed, and analyzed 383 gene variants of the viral HRSVgp04 gene ORF, with multiple synonymous mutations at key positions along the transcript in the eukaryote S. cerevisiae. Our results show that a few silent mutations at the 5'UTR can have a dramatic effect of up to 15 fold change on protein levels, and that even synonymous mutations in positions more than 120 nucleotides downstream from the ORF 5'end can modulate protein levels up to 160%-300%. We demonstrate that the correlation between protein levels and folding energy increases with the significance of the level of selection of the latter in endogenous genes, reinforcing the notion that selection for folding strength in different parts of the ORF is related to translation regulation. Our measured protein abundance correlates notably(correlation up to r = 0.62 (p=0.0013)) with mean relative codon decoding times, based on ribosomal densities (Ribo-Seq) in endogenous genes, supporting the conjecture that translation elongation and adaptation to the tRNA pool can modify protein levels in a causal/direct manner. This report provides an improved understanding of transcript evolution, design principles of gene expression regulation, and suggests simple rules for engineering synthetic gene expression in eukaryotes.

Multiple roles of the coding sequence 5' end in gene expression regulation

The codon composition of the coding sequence's (ORF) 5′ end first few dozen codons is known to be distinct to that of the rest of the ORF. Various explanations for the unusual codon distribution in this region have been proposed in recent years, and include, among others, novel regulatory mechanisms of translation initiation and elongation. However, due to the fact that many overlapping regulatory signals are suggested to be associated with this relatively short region, its research is challenging. Here, we review the currently known signals that appear in this region, the theories related to the way they regulate translation and affect the organismal fitness, and the debates they provoke.

Modular construction of a functional artificial epothilone polyketide pathway

Natural products of microbial origin continue to be an important source of pharmaceuticals and agrochemicals exhibiting potent activities and often novel modes of action. Due to their inherent structural complexity chemical synthesis is often hardly possible, leaving fermentation as the only viable production route. In addition, the pharmaceutical properties of natural products often need to be optimized for application by sophisticated medicinal chemistry and/or biosynthetic engineering. The latter requires a detailed understanding of the biosynthetic process and genetic tools to modify the producing organism that are often unavailable. Consequently, heterologous expression of complex natural product pathways has been in the focus of development over recent years. However, piecing together existing DNA cloned from natural sources and achieving efficient expression in heterologous circuits represent several limitations that can be addressed by synthetic biology. In this work we have redesigned and reassembled the 56 kb epothilone biosynthetic gene cluster from Sorangium cellulosum for expression in the high GC host Myxococcus xanthus. The codon composition was adapted to a modified codon table for M. xanthus, and unique restriction sites were simultaneously introduced and others eliminated from the sequence in order to permit pathway assembly and future interchangeability of modular building blocks from the epothilone megasynthetase. The functionality of the artificial pathway was demonstrated by successful heterologous epothilone production in M. xanthus at significant yields that have to be improved in upcoming work. Our study sets the stage for future engineering of epothilone biosynthesis and production optimization using a highly flexible assembly strategy.

A novel methodology for enhanced and consistent heterologous expression of unmodified human cytochrome P450 1B1 (CYP1B1)

Cytochrome P450 1B1 (CYP1B1) is a universal cancer marker and is implicated in many other disorders. Mutations in CYP1B1 are also associated with childhood blindness due to primary congenital glaucoma (PCG). To understand the CYP1B1 mediated etiopathology of PCG and pathomechanism of various cancers, it is important to carry out its functional studies. Heterologous expression of CYP1B1 in prokaryotes is imperative because bacteria yield a higher amount of heterologous proteins in lesser time and so the expressed protein is ideal for functional studies. In such expression system there is no interference by other eukaryotic proteins. But the story is not that simple as expression of heterologous CYP1B1 poses many technical difficulties. Investigators have employed various modifications/deletions of CYP N-terminus to improve CYP1B1 expression. However, the drawback of these studies is that it changes the original protein and, as a result, invalidates functional studies. The present study examines the role of various conditions and reagents in successful and consistent expression of sufficient quantities of unmodified/native human CYP1B1 in E. coli. We aimed at expressing CYP1B1 in various strains of E. coli and in the course developed a protocol that results in high expression of unmodified protein sufficient for functional/biophysical studies. We examined CYP1B1 expression with respect to different expression vectors, bacterial strains, types of culture media, time, Isopropyl β-D-1-thiogalactopyranoside concentrations, temperatures, rotations per minute, conditioning reagents and the efficacy of a newly described technique called double colony selection. We report a protocol that is simple, easy and can be carried out in any laboratory without the requirement of a fermentor. Though employed for CYP1B1 expression, this protocol can ideally be used to express any eukaryotic membrane protein.

Design rules for efficient transgene expression in plants

Sustained expression of transgenes in specified developmental patterns is commonly needed in plant biotechnology, but obstructed by transgene silencing. Here, we present a set of gene design rules, tested on the silencing-susceptible beetle luc and bacterial ims genes, expressed in sugarcane. Designs tested independently or in combination included removal of rare codons, removal of RNA instability sequences, blocking of likely endogenous sRNA binding sites and randomization of non-rare codons. Stable transgene expression analyses, on multiple independent lines per construct, showed greatest improvement from the removal of RNA instability sequences, accompanied by greatly reduced transcript degradation evident in northern blot analysis. We provide a set of motifs that readily can be eliminated concurrently with rare codons and undesired structural features such as repeat sequences, using Gene Designer 2.0 software. These design rules yielded 935- and 5-fold increased expression in transgenic callus, relative to the native luc and ims sequences; and gave sustained expression under the control of sugarcane and heterologous promoters over several years in greenhouse and field trials. The rules can be applied easily with codon usage tables from any plant species, providing a simple and effective means to achieve sustained expression of otherwise silencing-prone transgenes in plants.

Large-scale de novo DNA synthesis: technologies and applications

For over 60 years, the synthetic production of new DNA sequences has helped researchers understand and engineer biology. Here we summarize methods and caveats for the de novo synthesis of DNA, with particular emphasis on recent technologies that allow for large-scale and low-cost production. In addition, we discuss emerging applications enabled by large-scale de novo DNA constructs, as well as the challenges and opportunities that lie ahead.

Enhanced expression of GCRV VP6 in CIK cells by relative sequence optimization

Efficient expression of target protein is one of strategies for gene therapy or vaccine design. Many studies showed that codon optimization could enhance the expression of target proteins. In this paper, a target sequence of about 1.26 kb encoding the major capsid protein VP6 of grass carp reovirus (GCRV) and an optimized counterpart were synthesized and inserted into vectors for expressing VP6. The final constructs (named as pcDV6G and pcDV6YG) were transfected in Ctenopharyngodon idellus kidney (CIK) cells. The fluorescence analysis and the Western blot results showed that the gene fragment was transfected and expressed in CIK cells successfully. Although the qRT-PCR results showed no difference at the messenger RNA (mRNA) levels between the different versions of vp6 in the indicated stages, the enzyme-linked immunosorbent assay (ELISA) results showed that the protein level of VP6 expressed by pcDV6YG was higher than that by pcDV6G in the indicated hours. Taken together, these results suggest that the enhanced expression of GCRV VP6 in CIK cells by relative sequence optimization may be a good choice for making DNA vaccine against GCRV.

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.

Relative specificity: all substrates are not created equal

A biological molecule, e.g., an enzyme, tends to interact with its many cognate substrates, targets, or partners differentially. Such a property is termed relative specificity and has been proposed to regulate important physiological functions, even though it has not been examined explicitly in most complex biochemical systems. This essay reviews several recent large-scale studies that investigate protein folding, signal transduction, RNA binding, translation and transcription in the context of relative specificity. These results and others support a pervasive role of relative specificity in diverse biological processes. It is becoming clear that relative specificity contributes fundamentally to the diversity and complexity of biological systems, which has significant implications in disease processes as well.