Design parameters to control synthetic gene expression in Escherichia coli

Whole-cell biocatalysis for selective and productive C-O functional group introduction and modification

During the last decades, biocatalysis became of increasing importance for chemical and pharmaceutical industries. Regarding regio- and stereospecificity, enzymes have shown to be superior compared to traditional chemical synthesis approaches, especially in C-O functional group chemistry. Catalysts established on a process level are diverse and can be classified along a functional continuum starting with single-step biotransformations using isolated enzymes or microbial strains towards fermentative processes with recombinant microorganisms containing artificial synthetic pathways. The complex organization of respective enzymes combined with aspects such as cofactor dependency and low stability in isolated form often favors the use of whole cells over that of isolated enzymes. Based on an inventory of the large spectrum of biocatalytic C-O functional group chemistry, this review focuses on highlighting the potentials, limitations, and solutions offered by the application of self-regenerating microbial cells as biocatalysts. Different cellular functionalities are discussed in the light of their (possible) contribution to catalyst efficiency. The combined achievements in the areas of protein, genetic, metabolic, and reaction engineering enable the development of whole-cell biocatalysts as powerful tools in organic synthesis.

Modified 'one amino acid-one codon' engineering of high GC content TaqII-coding gene from thermophilic Thermus aquaticus results in radical expression increase

BACKGROUND

An industrial approach to protein production demands maximization of cloned gene expression, balanced with the recombinant host's viability. Expression of toxic genes from thermophiles poses particular difficulties due to high GC content, mRNA secondary structures, rare codon usage and impairing the host's coding plasmid replication.TaqII belongs to a family of bifunctional enzymes, which are a fusion of the restriction endonuclease (REase) and methyltransferase (MTase) activities in a single polypeptide. The family contains thermostable REases with distinct specificities: TspGWI, TaqII, Tth111II/TthHB27I, TspDTI and TsoI and a few enzymes found in mesophiles. While not being isoschizomers, the enzymes exhibit amino acid (aa) sequence homologies, having molecular sizes of ~120 kDa share common modular architecture, resemble Type-I enzymes, cleave DNA 11/9 nt from the recognition sites, their activity is affected by S-adenosylmethionine (SAM).

RESULTS

We describe the taqIIRM gene design, cloning and expression of the prototype TaqII. The enzyme amount in natural hosts is extremely low. To improve expression of the taqIIRM gene in Escherichia coli (E. coli), we designed and cloned a fully synthetic, low GC content, low mRNA secondary structure taqIIRM, codon-optimized gene under a bacteriophage lambda (λ) PR promoter. Codon usage based on a modified 'one amino acid-one codon' strategy, weighted towards low GC content codons, resulted in approximately 10-fold higher expression of the synthetic gene. 718 codons of total 1105 were changed, comprising 65% of the taqIIRM gene. The reason for we choose a less effective strategy rather than a resulting in high expression yields 'codon randomization' strategy, was intentional, sub-optimal TaqII in vivo production, in order to decrease the high 'toxicity' of the REase-MTase protein.

CONCLUSIONS

Recombinant wt and synthetic taqIIRM gene were cloned and expressed in E. coli. The modified 'one amino acid-one codon' method tuned for thermophile-coded genes was applied to obtain overexpression of the 'toxic' taqIIRM gene. The method appears suited for industrial production of thermostable 'toxic' enzymes in E. coli. This novel variant of the method biased toward increasing a gene's AT content may provide economic benefits for industrial applications.

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).

Advances in host and vector development for the production of plasmid DNA vaccines

Recent developments in DNA vaccine research provide a new momentum for this rather young and potentially disruptive technology. Gene-based vaccines are capable of eliciting protective immunity in humans to persistent intracellular pathogens, such as HIV, malaria, and tuberculosis, for which the conventional vaccine technologies have failed so far. The recent identification and characterization of genes coding for tumor antigens has stimulated the development of DNA-based antigen-specific cancer vaccines. Although most academic researchers consider the production of reasonable amounts of plasmid DNA (pDNA) for immunological studies relatively easy to solve, problems often arise during this first phase of production. In this chapter we review the current state of the art of pDNA production at small (shake flasks) and mid-scales (lab-scale bioreactor fermentations) and address new trends in vector design and strain engineering. We will guide the reader through the different stages of process design starting from choosing the most appropriate plasmid backbone, choosing the right Escherichia coli (E. coli) strain for production, and cultivation media and scale-up issues. In addition, we will address some points concerning the safety and potency of the produced plasmids, with special focus on producing antibiotic resistance-free plasmids. The main goal of this chapter is to make immunologists aware of the fact that production of the pDNA vaccine has to be performed with as much as attention and care as the rest of their research.

Ligand modifications to reduce the relative resistance of multi-drug resistant HIV-1 protease

Proper proteolytic processing of the HIV-1 Gag/Pol polyprotein is required for HIV infection and viral replication. This feature has made HIV-1 protease an attractive target for antiretroviral drug design for the treatment of HIV-1 infected patients. To examine the role of the P1 and P1'positions of the substrate in inhibitory efficacy of multi-drug resistant HIV-1 protease 769 (MDR 769), we performed a series of structure-function studies. Using the original CA/p2 cleavage site sequence, we generated heptapeptides containing one reduced peptide bond with an L to F and A to F double mutation at P1 and P1' (F-r-F), and an A to F at P1' (L-r-F) resulting in P1/P1' modified ligands. Here, we present an analysis of co-crystal structures of CA/p2 F-r-F, and CA/p2 L-r-F in complex with MDR 769. To examine conformational changes in the complex structure, molecular dynamic (MD) simulations were performed with MDR769-ligand complexes. MD trajectories show the isobutyl group of both the lopinavir analog and the CA/p2 L-r-F substrate cause a conformational change of in the active site of MDR 769. IC50 measurements suggest the non identical P1/P1' ligands (CA/p2 L-r-F and lopinavir analog) are more effective against MDR proteases as opposed to identical P1/P1'ligands. Our results suggest that a non identical P1/P1'composition may be more favorable for the inhibition of MDR 769 as they induce conformational changes in the active site of the enzyme resulting in disruption of the two-fold symmetry of the protease, thus, stabilizing the inhibitor in the active site.

One-pot DNA construction for synthetic biology: the Modular Overlap-Directed Assembly with Linkers (MODAL) strategy

Overlap-directed DNA assembly methods allow multiple DNA parts to be assembled together in one reaction. These methods, which rely on sequence homology between the ends of DNA parts, have become widely adopted in synthetic biology, despite being incompatible with a key principle of engineering: modularity. To answer this, we present MODAL: a Modular Overlap-Directed Assembly with Linkers strategy that brings modularity to overlap-directed methods, allowing assembly of an initial set of DNA parts into a variety of arrangements in one-pot reactions. MODAL is accompanied by a custom software tool that designs overlap linkers to guide assembly, allowing parts to be assembled in any specified order and orientation. The in silico design of synthetic orthogonal overlapping junctions allows for much greater efficiency in DNA assembly for a variety of different methods compared with using non-designed sequence. In tests with three different assembly technologies, the MODAL strategy gives assembly of both yeast and bacterial plasmids, composed of up to five DNA parts in the kilobase range with efficiencies of between 75 and 100%. It also seamlessly allows mutagenesis to be performed on any specified DNA parts during the process, allowing the one-step creation of construct libraries valuable for synthetic biology applications.

Tuning response curves for synthetic biology

Synthetic biology may be viewed as an effort to establish, formalize, and develop an engineering discipline in the context of biological systems. The ability to tune the properties of individual components is central to the process of system design in all fields of engineering, and synthetic biology is no exception. A large and growing number of approaches have been developed for tuning the responses of cellular systems, and here we address specifically the issue of tuning the rate of response of a system: given a system where an input affects the rate of change of an output, how can the shape of the response curve be altered experimentally? This affects a system's dynamics as well as its steady-state properties, both of which are critical in the design of systems in synthetic biology, particularly those with multiple components. We begin by reviewing a mathematical formulation that captures a broad class of biological response curves and use this to define a standard set of varieties of tuning: vertical shifting, horizontal scaling, and the like. We then survey the experimental literature, classifying the results into our defined categories, and organizing them by regulatory level: transcriptional, post-transcriptional, and post-translational.

Causes and effects of N-terminal codon bias in bacterial genes

Most amino acids are encoded by multiple codons, and codon choice has strong effects on protein expression. Rare codons are enriched at the N terminus of genes in most organisms, although the causes and effects of this bias are unclear. Here, we measure expression from >14,000 synthetic reporters in Escherichia coli and show that using N-terminal rare codons instead of common ones increases expression by ~14-fold (median 4-fold). We quantify how individual N-terminal codons affect expression and show that these effects shape the sequence of natural genes. Finally, we demonstrate that reduced RNA structure and not codon rarity itself is responsible for expression increases. Our observations resolve controversies over the roles of N-terminal codon bias and suggest a straightforward method for optimizing heterologous gene expression in bacteria.

Building better drugs: developing and regulating engineered therapeutic proteins

Most native proteins do not make optimal drugs and thus a second- and third-generation of therapeutic proteins, which have been engineered to improve product attributes or to enhance process characteristics, are rapidly becoming the norm. There has been unprecedented progress, during the past decade, in the development of platform technologies that further these ends. Although the advantages of engineered therapeutic proteins are considerable, the alterations can affect the safety and efficacy of the drugs. We discuss both the key technological innovations with respect to engineered therapeutic proteins and advancements in the underlying basic science. The latter would permit the design of science-based criteria for the prediction and assessment of potential risks and the development of appropriate risk management plans. This in turn holds promise for more predictable criteria for the licensure of a class of products that are extremely challenging to develop but represent an increasingly important component of modern medical practice.

Genome-scale engineering for systems and synthetic biology

This review provides an overview of methodologies and technologies enabling genome-scale engineering, focusing on the design, construction, and testing of modified genomes in a variety of organisms. Future applications for systems and synthetic biology are discussed.

Genome-modification technologies enable the rational engineering and perturbation of biological systems. Historically, these methods have been limited to gene insertions or mutations at random or at a few pre-defined locations across the genome. The handful of methods capable of targetedgene editing suffered from low efficiencies, significant labor costs, or both. Recent advances have dramatically expanded our ability to engineer cells in a directed and combinatorial manner. Here, we review current technologies and methodologies for genome-scale engineering, discuss the prospects for extending efficient genome modification to new hosts, and explore the implications of continued advances toward the development of flexibly programmable chasses, novel biochemistries, and safer organismal and ecological engineering.

Differential translation tunes uneven production of operon-encoded proteins

Clustering of functionally related genes in operons allows for coregulated gene expression in prokaryotes. This is advantageous when equal amounts of gene products are required. Production of protein complexes with an uneven stoichiometry, however, requires tuning mechanisms to generate subunits in appropriate relative quantities. Using comparative genomic analysis, we show that differential translation is a key determinant of modulated expression of genes clustered in operons and that codon bias generally is the best in silico indicator of unequal protein production. Variable ribosome density profiles of polycistronic transcripts correlate strongly with differential translation patterns. In addition, we provide experimental evidence that de novo initiation of translation can occur at intercistronic sites, allowing for differential translation of any gene irrespective of its position on a polycistronic messenger. Thus, modulation of translation efficiency appears to be a universal mode of control in bacteria and archaea that allows for differential production of operon-encoded proteins.

Genetic code-guided protein synthesis and folding in Escherichia coli

Universal genetic codes are degenerated with 61 codons specifying 20 amino acids, thus creating synonymous codons for a single amino acid. Synonymous codons have been shown to affect protein properties in a given organism. To address this issue and explore how Escherichia coli selects its "codon-preferred" DNA template(s) for synthesis of proteins with required properties, we have designed synonymous codon libraries based on an antibody (scFv) sequence and carried out bacterial expression and screening for variants with altered properties. As a result, 342 codon variants have been identified, differing significantly in protein solubility and functionality while retaining the identical original amino acid sequence. The soluble expression level varied from completely insoluble aggregates to a soluble yield of ~2.5 mg/liter, whereas the antigen-binding activity changed from no binding at all to a binding affinity of > 10(-8) m. Not only does our work demonstrate the involvement of genetic codes in regulating protein synthesis and folding but it also provides a novel screening strategy for producing improved proteins without the need to substitute amino acids.

Rate-limiting steps in yeast protein translation

Deep sequencing now provides detailed snapshots of ribosome occupancy on mRNAs. We leverage these data to parameterize a computational model of translation, keeping track of every ribosome, tRNA, and mRNA molecule in a yeast cell. We determine the parameter regimes in which fast initiation or high codon bias in a transgene increases protein yield and infer the initiation rates of endogenous Saccharomyces cerevisiae genes, which vary by several orders of magnitude and correlate with 5′ mRNA folding energies. Our model recapitulates the previously reported 5′-to-3′ ramp of decreasing ribosome densities, although our analysis shows that this ramp is caused by rapid initiation of short genes rather than slow codons at the start of transcripts. We conclude that protein production in healthy yeast cells is typically limited by the availability of free ribosomes, whereas protein production under periods of stress can sometimes be rescued by reducing initiation or elongation rates.

•

Computational model of translation tracks all ribosomes, tRNAs, and mRNAs in a cell

•

Translation is generally limited by initiation, not elongation

•

Model allows inference of initiation rates for all yeast genes

•

Ramp of 5′ ribosomes is caused by rapid initiation of short genes

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.

Physical constraints on biological integral control design for homeostasis and sensory adaptation

Synthetic biology includes an effort to use design-based approaches to create novel controllers, biological systems aimed at regulating the output of other biological processes. The design of such controllers can be guided by results from control theory, including the strategy of integral feedback control, which is central to regulation, sensory adaptation, and long-term robustness. Realization of integral control in a synthetic network is an attractive prospect, but the nature of biochemical networks can make the implementation of even basic control structures challenging. Here we present a study of the general challenges and important constraints that will arise in efforts to engineer biological integral feedback controllers or to analyze existing natural systems. Constraints arise from the need to identify target output values that the combined process-plus-controller system can reach, and to ensure that the controller implements a good approximation of integral feedback control. These constraints depend on mild assumptions about the shape of input-output relationships in the biological components, and thus will apply to a variety of biochemical systems. We summarize our results as a set of variable constraints intended to provide guidance for the design or analysis of a working biological integral feedback controller.

Enhanced expression of codon optimized interferon gamma in CHO cells

The human interferon-gamma (IFN-γ) is a potential drug candidate for treating various diseases due to its immunomodulatory properties. The efficient production of this protein can be achieved through a popular industrial host, Chinese hamster ovary (CHO) cells. However, recombinant expression of foreign proteins is typically suboptimal possibly due to the usage of non-native codon patterns within the coding sequence. Therefore, we demonstrated the application of a recently developed codon optimization approach to design synthetic IFN-γ coding sequences for enhanced heterologous expression in CHO cells. For codon optimization, earlier studies suggested to establish the target usage distribution pattern in terms of selected design parameters such as individual codon usage (ICU) and codon context (CC), mainly based on the host's highly expressed genes. However, our RNA-Seq based transcriptome profiling indicated that the ICU and CC distribution patterns of different gene expression classes in CHO cell are relatively similar, unlike other microbial expression hosts, Escherichia coli and Saccharomyces cerevisiae. This finding was further corroborated through the in vivo expression of various ICU and CC optimized IFN-γ in CHO cells. Interestingly, the CC-optimized genes exhibited at least 13-fold increase in expression level compared to the wild-type IFN-γ while a maximum of 10-fold increase was observed for the ICU-optimized genes. Although design criteria based on individual codons, such as ICU, have been widely used for gene optimization, our experimental results suggested that codon context is relatively more effective parameter for improving recombinant IFN-γ expression in CHO cells.

Translational sensitivity of the Escherichia coli genome to fluctuating tRNA availability

The synthesis of protein from messenger RNA during translation is a highly dynamic process that plays a key role in controlling the efficiency and fidelity of genome-wide protein expression. The availability of aminoacylated transfer RNA (tRNA) is a major factor influencing the speed of ribosomal movement, which depending on codon choices, varies considerably along a transcript. Furthermore, it has been shown experimentally that tRNA availability can vary significantly under different growth and stress conditions, offering the cell a way to adapt translational dynamics across the genome. Existing models of translation have neglected fluctuations of tRNA pools, instead assuming fixed tRNA availabilities over time. This has lead to an incomplete understanding of this process. Here, we show for the entire Escherichia coli genome how and to what extent translational speed profiles, which capture local aspects of translational elongation, respond to measured shifts in tRNA availability. We find that translational profiles across the genome are affected to differing degrees, with genes that are essential or related to fundamental processes such as translation, being more robust than those linked to regulation. Furthermore, we reveal how fluctuating tRNA availability influences profiles of specific sequences known to play a significant role in translational control of gene expression.

Improved blue, green, and red fluorescent protein tagging vectors for S. cerevisiae

Fluorescent protein fusions are a powerful tool to monitor the localization and trafficking of proteins. Such studies are particularly easy to carry out in the budding yeast Saccharomyces cerevisiae due to the ease with which tags can be introduced into the genome by homologous recombination. However, the available yeast tagging plasmids have not kept pace with the development of new and improved fluorescent proteins. Here, we have constructed yeast optimized versions of 19 different fluorescent proteins and tested them for use as fusion tags in yeast. These include two blue, seven green, and seven red fluorescent proteins, which we have assessed for brightness, photostability and perturbation of tagged proteins. We find that EGFP remains the best performing green fluorescent protein, that TagRFP-T and mRuby2 outperform mCherry as red fluorescent proteins, and that mTagBFP2 can be used as a blue fluorescent protein tag. Together, the new tagging vectors we have constructed provide improved blue and red fluorescent proteins for yeast tagging and three color imaging.

Non-Invasive Analysis of Recombinant mRNA Stability in Escherichia coli by a Combination of Transcriptional Inducer Wash-Out and qRT-PCR

mRNA stability is one among many parameters that can potentially affect the level of recombinant gene expression in bacteria. Blocking of the entire prokaryotic transcription machinery by addition of rifampicin is commonly used in protocols for analysis of mRNA stability. Here we show that such treatment can be effectively replaced by a simple, non-invasive method based on removal of the relevant transcriptional inducers and that the mRNA decay can then be followed by qRT-PCR. To establish the methodology we first used the m-toluate-inducible XylS/Pm expression cassette as a model system and analyzed several examples of DNA modifications causing gene expression stimulation in Escherichia coli. The new method allowed us to clearly discriminate whether an improvement in mRNA stability contributes to observed increases in transcript amounts for each individual case. To support the experimental data a simple mathematical fitting model was developed to calculate relative decay rates. We extended the relevance of the method by demonstrating its application also for an IPTG-inducible expression cassette (LacI/P_tac) and by analyzing features of the bacteriophage T7-based expression system. The results suggest that the methodology is useful in elucidating factors controlling mRNA stability as well as other specific features of inducible expression systems. Moreover, as expression systems based on diffusible inducers are almost universally available, the concept can be most likely used to measure mRNA decay for any gene in any cell type that is heavily used in molecular biology research.

Directed evolution study unveiling key sequence factors that affect translation efficiency in Escherichia coli

Synonymous mutations in protein coding genes significantly impact translation efficiency. We synthesized a pair of genes encoding green fluorescent protein that were separated by 160 synonymous mutations to investigate key factors that affect translation efficiency. One sequence was optimized for Escherichia coli (GFP(Eco)) and the other for Bacillus subtilis (GFP(Bsu)). When the genes were expressed in E. coli, GFP(Eco) fluoresced 12-fold stronger than GFP(Bsu), confirming the suboptimal nature of the GFP(Bsu) gene. We then employed directed evolution to improve the expression of GFP(Bsu). Random mutagenesis and DNA shuffling was used to generate mutant libraries, which were screened for fluorescence. A variant showing 6-fold fluorescence enhancement was identified, which contained a single mutation (G10A) in a rare codon for Gly-4. However, the substitution generated another type of rare codon, AGA, for Arg, suggesting that the improvement was caused by a factor other than the rare codon. We next applied saturation mutagenesis to Gly-4. The darkest variant contained a GGG codon (GFP(Bsu)-G) for Gly-4. Taking the location of the mutation into account, we hypothesized that destabilization of the mRNA secondary structure around the initiation codon improved the expression. We then randomized the nucleotide triplet in 5'-untranslated region (5'UTR) of GFP(Bsu), which is complementary to the Gly-4 codon. A variant showing 6-fold fluorescence enhancement was identified, which exhibited a destabilized secondary structure. When this 5'UTR sequence was combined with GFP(Bsu)-G, 22-fold fluorescent improvement was achieved. Collectively, the stability of the mRNA secondary structure around the initiation codon predominantly affected the translation efficiency.

Antigenic differences between the EG95-related proteins from Echinococcus granulosus G1 and G6 genotypes: implications for vaccination

Cystic echinococcosis caused by Echinococcus granulosus remains an important and neglected issue in public health. The study of the likely efficacy of the currently available EG95 vaccine against other genotypes of the parasite is important to improve the vaccine as a potential tool to be used in control programmes. The recombinant vaccine EG95-1G1 was developed based on the G1 genotype of E. granulosus. Characterization of the eg95 gene family in the G6 genotype by genomic DNA cloning previously produced the first unequivocal information about the composition of the gene family in a different genotype. The information was used in this study to predict and express two EG95-related proteins from the G6 genotype as recombinants, for assessment of their capacity to bind antibodies raised in sheep vaccinated with the EG95-1G1 vaccine. The proteins (EG95-1G6 and EG95-5G6) from the G6 genotype of E. granulosus were unable to bind all the antibodies raised by sheep vaccinated with EG95-1G1. Differences in the amino acid sequence of EG95-related proteins from G6 and likely the differences in the encoded FnIII domain may be responsible for changes in the conformation of these epitopes.

Towards a whole-cell modeling approach for synthetic biology

Despite rapid advances over the last decade, synthetic biology lacks the predictive tools needed to enable rational design. Unlike established engineering disciplines, the engineering of synthetic gene circuits still relies heavily on experimental trial-and-error, a time-consuming and inefficient process that slows down the biological design cycle. This reliance on experimental tuning is because current modeling approaches are unable to make reliable predictions about the in vivo behavior of synthetic circuits. A major reason for this lack of predictability is that current models view circuits in isolation, ignoring the vast number of complex cellular processes that impinge on the dynamics of the synthetic circuit and vice versa. To address this problem, we present a modeling approach for the design of synthetic circuits in the context of cellular networks. Using the recently published whole-cell model of Mycoplasma genitalium, we examined the effect of adding genes into the host genome. We also investigated how codon usage correlates with gene expression and find agreement with existing experimental results. Finally, we successfully implemented a synthetic Goodwin oscillator in the whole-cell model. We provide an updated software framework for the whole-cell model that lays the foundation for the integration of whole-cell models with synthetic gene circuit models. This software framework is made freely available to the community to enable future extensions. We envision that this approach will be critical to transforming the field of synthetic biology into a rational and predictive engineering discipline.

Tuning the dials of Synthetic Biology

Synthetic Biology is the ‘Engineering of Biology’ – it aims to use a forward-engineering design cycle based on specifications, modelling, analysis, experimental implementation, testing and validation to modify natural or design new, synthetic biology systems so that they behave in a predictable fashion. Motivated by the need for truly plug-and-play synthetic biological components, we present a comprehensive review of ways in which the various parts of a biological system can be modified systematically. In particular, we review the list of ‘dials’ that are available to the designer and discuss how they can be modelled, tuned and implemented. The dials are categorized according to whether they operate at the global, transcriptional, translational or post-translational level and the resolution that they operate at. We end this review with a discussion on the relative advantages and disadvantages of some dials over others.

Rational design of orthogonal libraries of protein coding genes

Array-based oligonucleotide synthesis technologies provide access to thousands of custom-designed sequence variants at low cost. Large-scale synthesis and high-throughput assays have become valuable experimental tools to study in detail the interplay between sequence and function. We have developed a methodology and corresponding algorithms for the design of diverse protein coding gene libraries, to exploit the potential of multiplex synthesis and help elucidate the effects of codon utilization and other factors in gene expression. Using our algorithm, we have computationally designed gene libraries with hundreds to thousands of orthogonal codon usage variants, uniformly exploring the design space of codon utilization, while demanding only a small fraction of the synthesis cost that would be required if these variants were synthesized independently.

Fitness costs of synonymous mutations in the rpsT gene can be compensated by restoring mRNA base pairing

We previously reported that the distribution of fitness effects for non-synonymous and synonymous mutations in Salmonella typhimurium ribosomal proteins S20 and L1 are similar, suggesting that fitness constraints are present at the level of mRNA. Here we explore the hypothesis that synonymous mutations confer their fitness-reducing effect by alterating the secondary structure of the mRNA. To this end, we constructed a set of synonymous substitutions in the rpsT gene, encoding ribosomal protein S20, that are located in predicted paired regions in the mRNA and measured their effect on bacterial fitness. Our results show that for 3/9 cases tested, the reduced fitness conferred by a synonymous mutation could be fully or partly restored by introducing a second synonymous substitution that restore base pairing in a mRNA stem. In addition, random mutations in predicted paired regions had larger fitness effects than those in unpaired regions. Finally, we did not observe any correlation between fitness effects of the synonymous mutations and their rarity. These results suggest that for ribosomal protein S20, the deleterious effects of synonymous mutations are not generally due to codon usage effects, but that mRNA secondary structure is a major fitness constraint.

Current state and recent advances in biopharmaceutical production in Escherichia coli, yeasts and mammalian cells

Almost all of the 200 or so approved biopharmaceuticals have been produced in one of three host systems: the bacterium Escherichia coli, yeasts (Saccharomyces cerevisiae, Pichia pastoris) and mammalian cells. We describe the most widely used methods for the expression of recombinant proteins in the cytoplasm or periplasm of E. coli, as well as strategies for secreting the product to the growth medium. Recombinant expression in E. coli influences the cell physiology and triggers a stress response, which has to be considered in process development. Increased expression of a functional protein can be achieved by optimizing the gene, plasmid, host cell, and fermentation process. Relevant properties of two yeast expression systems, S. cerevisiae and P. pastoris, are summarized. Optimization of expression in S. cerevisiae has focused mainly on increasing the secretion, which is otherwise limiting. P. pastoris was recently approved as a host for biopharmaceutical production for the first time. It enables high-level protein production and secretion. Additionally, genetic engineering has resulted in its ability to produce recombinant proteins with humanized glycosylation patterns. Several mammalian cell lines of either rodent or human origin are also used in biopharmaceutical production. Optimization of their expression has focused on clonal selection, interference with epigenetic factors and genetic engineering. Systemic optimization approaches are applied to all cell expression systems. They feature parallel high-throughput techniques, such as DNA microarray, next-generation sequencing and proteomics, and enable simultaneous monitoring of multiple parameters. Systemic approaches, together with technological advances such as disposable bioreactors and microbioreactors, are expected to lead to increased quality and quantity of biopharmaceuticals, as well as to reduced product development times.

Metabolic reconstruction and flux analysis of industrial Pichia yeasts

Pichia yeasts have been recognized as important microbial cell factories in the biotechnological industry. Notably, the Pichia pastoris and Pichia stipitis species have attracted much research interest due to their unique cellular physiology and metabolic capability: P. pastoris has the ability to utilize methanol for cell growth and recombinant protein production, while P. stipitis is capable of assimilating xylose to produce ethanol under oxygen-limited conditions. To harness these characteristics for biotechnological applications, it is highly required to characterize their metabolic behavior. Recently, following the genome sequencing of these two Pichia species, genome-scale metabolic networks have been reconstructed to model the yeasts' metabolism from a systems perspective. To date, there are three genome-scale models available for each of P. pastoris and P. stipitis. In this mini-review, we provide an overview of the models, discuss certain limitations of previous studies, and propose potential future works that can be conducted to better understand and engineer Pichia yeasts for industrial applications.

Good codons, bad transcript: large reductions in gene expression and fitness arising from synonymous mutations in a key enzyme

Biased codon usage in protein-coding genes is pervasive, whereby amino acids are largely encoded by a specific subset of possible codons. Within individual genes, codon bias is stronger at evolutionarily conserved residues, favoring codons recognized by abundant tRNAs. Although this observation suggests an overall pattern of selection for translation speed and/or accuracy, other work indicates that transcript structure or binding motifs drive codon usage. However, our understanding of codon bias evolution is constrained by limited experimental data on the fitness effects of altering codons in functional genes. To bridge this gap, we generated synonymous variants of a key enzyme-coding gene in Methylobacterium extorquens. We found that mutant gene expression, enzyme production, enzyme activity, and fitness were all significantly lower than wild-type. Surprisingly, encoding the gene using only rare codons decreased fitness by 40%, whereas an allele coded entirely by frequent codons decreased fitness by more than 90%. Increasing gene expression restored mutant fitness to varying degrees, demonstrating that the fitness disadvantage of synonymous mutants arose from a lack of beneficial protein rather than costs of protein production. Protein production was negatively correlated with the frequency of motifs with high affinity for the anti-Shine-Dalgarno sequence, suggesting ribosome pausing as the dominant cause of low mutant fitness. Together, our data support the idea that, although a particular set of codons are favored on average across a genome, in an individual gene selection can either act for or against codons depending on their local context.

Design and assembly of large synthetic DNA constructs

The availability of custom synthetic gene-length DNA products removes numerous bottlenecks in research efforts, making gene synthesis an increasingly common commercial service. However, the assembly of synthetic oligonucleotides into large, custom DNA constructs is not especially difficult, and performing "in-house" gene synthesis has time and cost advantages. This unit will treat both the concerns of design and physical assembly in gene synthesis, including how to design DNA sequences for synthesis and the design of overlapping oligonucleotide schemes to ensure facile assembly into the final product. Assembly is accomplished using a reliable series of PCR reactions, with a troubleshooting assembly protocol included, which not only assembles difficult sequences but allows identification of the source of a failure down to a pair of oligonucleotides.

On the steady-state distribution in the homogeneous ribosome flow model

A central biological process in all living organisms is gene translation. Developing a deeper understanding of this complex process may have ramifications to almost every biomedical discipline. Reuveni et al. recently proposed a new computational model of gene translation called the Ribosome Flow Model (RFM). In this paper, we consider a particular case of this model, called the Homogeneous Ribosome Flow Model (HRFM). From a biological viewpoint, this corresponds to the case where the transition rates of all the coding sequence codons are identical. This regime has been suggested recently based on experiments in mouse embryonic cells. We consider the steady-state distribution of the HRFM. We provide formulas that relate the different parameters of the model in steady state. We prove the following properties: 1) the ribosomal density profile is monotonically decreasing along the coding sequence; 2) the ribosomal density at each codon monotonically increases with the initiation rate; and 3) for a constant initiation rate, the translation rate monotonically decreases with the length of the coding sequence. In addition, we analyze the translation rate of the HRFM at the limit of very high and very low initiation rate, and provide explicit formulas for the translation rate in these two cases. We discuss the relationship between these theoretical results and biological findings on the translation process.

Codon preference optimization increases prokaryotic cystatin C expression

Gene expression is closely related to optimal vector-host system pairing in many prokaryotes. Redesign of the human cystatin C (cysC) gene using the preferred codons of the prokaryotic system may significantly increase cysC expression in Escherichia coli (E. coli). Specifically, cysC expression may be increased by removing unstable sequences and optimizing GC content. According to E. coli expression system codon preferences, the gene sequence was optimized while the amino acid sequence was maintained. The codon-optimized cysC (co-cysC) and wild-type cysC (wt-cysC) were expressed by cloning the genes into a pET-30a plasmid, thus transforming the recombinant plasmid into E. coli BL21. Before and after the optimization process, the prokaryotic expression vector and host bacteria were examined for protein expression and biological activation of CysC. The recombinant proteins in the lysate of the transformed bacteria were purified using Ni(2+)-NTA resin. Recombinant protein expression increased from 10% to 46% based on total protein expression after codon optimization. Recombinant CysC purity was above 95%. The significant increase in cysC expression in E. coli expression produced by codon optimization techniques may be applicable to commercial production systems.

Pseudomonas syringae self-protection from tabtoxinine-β-lactam by ligase TblF and acetylase Ttr

Plant pathogenic Pseudomonas syringae produce the hydroxy-β-lactam antimetabolite tabtoxinine-β-lactam (TβL) as a time-dependent inactivating glutamine analogue of plant glutamine synthetases. The producing pseudomonads use multiple modes of self-protection, two of which are characterized in this study. The first is the dipeptide ligase TblF which converts tabtoxinine-β-lactam to the TβL-Thr dipeptide known as tabtoxin. The dipeptide is not recognized by glutamine synthetase. This represents a Trojan Horse strategy: the dipeptide is secreted, taken up by dipeptide permeases in neighboring cells, and TβL is released by peptidase action. The second self-protection mode is elaboration by the acetyltransferase Ttr, which acetylates the α-amino group of the proximal inactivator TβL, but not the tabtoxin dipeptide.

Effects of codon sequence on the dynamics of genetic networks

In prokaryotes, the rate at which codons are translated varies from one codon to the next. Using a stochastic model of transcription and translation at the nucleotide and codon levels, we investigate the effects of the codon sequence on the dynamics of protein numbers. For sequences generated according to the codon frequencies in Escherichia coli, we find that mean protein numbers at near equilibrium differ with the codon sequence, due to the mean codon translation efficiencies, in particular of the codons at the ribosome binding site region. We find close agreement between these predictions and measurements of protein expression levels as a function of the codon sequence. Next, we investigate the effects of short codon sequences at the start/end of the RNA sequence with linearly increasing/decreasing translation efficiencies, known as slow ramps. The ramps affect the mean, but not the fluctuations, in proteins numbers by affecting the rate of translation initiation. Finally, we show that slow ramps affect the dynamics of small genetic circuits, namely, switches and clocks. In switches, ramps affect the frequency of switching and bias the robustness of the noisy attractors. In repressilators, ramps alter the robustness of periodicity. We conclude that codon sequences affect the dynamics of gene expression and genetic circuits and, thus, are likely to be under selection regarding both mean codon frequency as well as spatial arrangement along the sequence.

A selective force favoring increased G+C content in bacterial genes

Bacteria display considerable variation in their overall base compositions, which range from 13% to over 75% G+C. This variation in genomic base compositions has long been considered to be a strictly neutral character, due solely to differences in the mutational process; however, recent sequence comparisons indicate that mutational input alone cannot produce the observed base compositions, implying a role for natural selection. Because bacterial genomes have high gene content, forces that operate on the base composition of individual genes could help shape the overall genomic base composition. To explore this possibility, we tested whether genes that encode the same protein but vary only in their base compositions at synonymous sites have effects on bacterial fitness. Escherichia coli strains harboring G+C-rich versions of genes display higher growth rates, indicating that despite a pervasive mutational bias toward A+T, a selective force, independent of adaptive codon use, is driving genes toward higher G+C contents.

The PLOS ONE synthetic biology collection: six years and counting

Since it was launched in 2006, PLOS ONE has published over fifty articles illustrating the many facets of the emerging field of synthetic biology. This article reviews these publications by organizing them into broad categories focused on DNA synthesis and assembly techniques, the development of libraries of biological parts, the use of synthetic biology in protein engineering applications, and the engineering of gene regulatory networks and metabolic pathways. Finally, we review articles that describe enabling technologies such as software and modeling, along with new instrumentation. In order to increase the visibility of this body of work, the papers have been assembled into the PLOS ONE Synthetic Biology Collection (www.ploscollections.org/synbio). Many of the innovative features of the PLOS ONE web site will help make this collection a resource that will support a lively dialogue between readers and authors of PLOS ONE synthetic biology papers. The content of the collection will be updated periodically by including relevant articles as they are published by the journal. Thus, we hope that this collection will continue to meet the publishing needs of the synthetic biology community.

EuGene: maximizing synthetic gene design for heterologous expression

Numerous software applications exist to deal with synthetic gene design, granting the field of heterologous expression a significant support. However, their dispersion requires the access to different tools and online services in order to complete one single project. Analyzing codon usage, calculating codon adaptation index (CAI), aligning orthologs and optimizing genes are just a few examples. A software application, EuGene, was developed for the optimization of multiple gene synthetic design algorithms. In a seamless automatic form, EuGene calculates or retrieves genome data on codon usage (relative synonymous codon usage and CAI), codon context (CPS and codon pair bias), GC content, hidden stop codons, repetitions, deleterious sites, protein primary, secondary and tertiary structures, gene orthologs, species housekeeping genes, performs alignments and identifies genes and genomes. The main function of EuGene is analyzing and redesigning gene sequences using multi-objective optimization techniques that maximize the coding features of the resulting sequence.

AVAILABILITY

EuGene is freely available for non-commercial use, at http://bioinformatics.ua.pt/eugene.

Design and biophysical characterization of a monomeric four-alpha-helix bundle protein Aα₄ with affinity for the volatile anesthetic halothane

A monomeric four-α-helix bundle protein Aα₄ was designed as a step towards investigating the interaction of volatile general anesthetics with their putative membrane protein targets. The alpha helices, connected by glycine loops, have the sequence A, B, B', A'. The DNA sequence was designed to make the helices with the same amino acid sequences (helix A and A', B and B', respectively) as different as possible, while using codons which are favorable for expression in E. coli. The protein was bacterially expressed and purified to homogeneity using reversed-phase HPLC. Protein identity was verified using MALDI-TOF mass spectrometry. Far-UV circular dichroism spectroscopy confirmed the predominantly alpha-helical nature of the protein Aα₄. Guanidinium chloride induced denaturation showed that the monomeric four-α-helix bundle protein Aα₄ is considerably more stable compared to the dimeric di-α-helical protein (Aα₂-L38M)₂. The sigmoidal character of the unfolding reaction is conserved while the sharpness of the transition is increased 1.8-fold. The monomeric four-α-helix bundle protein Aα₄ bound halothane with a dissociation constant (K(d)) of 0.93 ± 0.02mM, as shown by both tryptophan fluorescence quenching and isothermal titration calorimetry. This monomeric four-α-helix bundle protein can now be used as a scaffold to incorporate natural central nervous system membrane protein sequences in order to examine general anesthetic interactions with putative targets in detail.

Efficient manipulations of synonymous mutations for controlling translation rate: an analytical approach

Gene translation is a central process in all living organism with important ramifications to almost every biomedical field. Previous systems evolutionary studies in the field have demonstrated that in many organisms coding sequence features undergo selection to optimize this process. In the current study, we report for the first time analytical proofs related to the various aspects of this process and its optimality. Among our results we show that coding sequences with mono- tonic increasing profiles of translation efficiency (i.e., with slower codons near the 5'UTR), mathematically optimize ribosomal allocation by minimizing the number of ribosomes needed for translating a codon per time unit. Thus, the genomic translation efficiency profile reported in previous studies for many organisms is optimal in this sense. In addition, we show that improving translation efficiency of a codon in a gene may result in a decrease in the translation rate of other genes, demonstrating that the relation between codon bias and protein translation rate is less trivial than was assumed before. Based on these observations we describe an efficient heuristic for designing coding sequences with specific translation efficiency and minimal ribosomal allocation for heterologous gene expression. We demonstrate how this heuristic can be used in biotechnology for engineering a heterologous gene before expressing it in a new host.

Expression of lignocellulolytic enzymes in Pichia pastoris

BACKGROUND

Sustainable utilization of plant biomass as renewable source for fuels and chemical building blocks requires a complex mixture of diverse enzymes, including hydrolases which comprise the largest class of lignocellulolytic enzymes. These enzymes need to be available in large amounts at a low price to allow sustainable and economic biotechnological processes.Over the past years Pichia pastoris has become an attractive host for the cost-efficient production and engineering of heterologous (eukaryotic) proteins due to several advantages.

RESULTS

In this paper codon optimized genes and synthetic alcohol oxidase 1 promoter variants were used to generate Pichia pastoris strains which individually expressed cellobiohydrolase 1, cellobiohydrolase 2 and beta-mannanase from Trichoderma reesei and xylanase A from Thermomyces lanuginosus. For three of these enzymes we could develop strains capable of secreting gram quantities of enzyme per liter in fed-batch cultivations. Additionally, we compared our achieved yields of secreted enzymes and the corresponding activities to literature data.

CONCLUSION

In our experiments we could clearly show the importance of gene optimization and strain characterization for successfully improving secretion levels. We also present a basic guideline how to correctly interpret the interplay of promoter strength and gene dosage for a successful improvement of the secretory production of lignocellulolytic enzymes in Pichia pastoris.