- Open Access
Comparison of two codon optimization strategies to enhance recombinant protein production in Escherichia coli
Microbial Cell Factoriesvolume 10, Article number: 15 (2011)
Variations in codon usage between species are one of the major causes affecting recombinant protein expression levels, with a significant impact on the economy of industrial enzyme production processes. The use of codon-optimized genes may overcome this problem. However, designing a gene for optimal expression requires choosing from a vast number of possible DNA sequences and different codon optimization methods have been used in the past decade. Here, a comparative study of the two most common methods is presented using calf prochymosin as a model.
Seven sequences encoding calf prochymosin have been designed, two using the "one amino acid-one codon" method and five using a "codon randomization" strategy. When expressed in Escherichia coli, the variants optimized by the codon randomization approach produced significantly more proteins than the native sequence including one gene that produced an increase of 70% in the amount of prochymosin accumulated. On the other hand, no significant improvement in protein expression was observed for the variants designed with the one amino acid-one codon method. The use of codon-optimized sequences did not affect the quality of the recovered inclusion bodies.
The results obtained in this study indicate that the codon randomization method is a superior strategy for codon optimization. A significant improvement in protein expression was obtained for the largely established process of chymosin production, showing the power of this strategy to reduce production costs of industrial enzymes in microbial hosts.
Industrial enzymes, included those used in food industry, are now traded as commodity products and there is a continuing need to reduce manufacturing costs in order to remain competitive in the global markets. Escherichia coli is a preferred host for the production of recombinant proteins because it combines fast growth rate, inexpensive fermentation media and well understood genetics; and the cost of production in this microorganism depends in large part upon the protein expression levels [1, 2].
Species-specific variations in codon usage are often cited as one of the major causes impacting protein expression levels [3, 4]. The presence of rare codons, which are correlated with low levels of their cognate tRNAs species in the cell, can reduce the translation rate and induce translation errors with a significant impact on the economy of the production process [5, 6]. In the past decade, a high number of genes have been re-designed to increase their expression level [1–3, 7–10]. However, designing a gene for optimal expression requires choosing from a large space containing a vast number of possible DNA sequences. Typically, two strategies have been used for codon optimization. The first one, known as "one amino acid-one codon", assigns the most abundant codon of the host or a set of selected genes to all instances of a given amino acid in the target sequence [4, 8, 11–15]. The second one, designed here "codon randomization", uses translation tables, based on the frequency distribution of the codons in an entire genome or a subset of highly expressed genes, to attach weights to each codon. In this case, codons are assigned randomly with a probability given by the weights [2, 7, 8, 16, 17].
Many transgenic proteins expressed in E. coli are recovered as insoluble aggregates in the form of inclusion bodies. The formation of these aggregates seems to be independent of the type of protein, and this drawback has been proven difficult to overcome [18, 19]. Nevertheless, the fact that inclusion bodies are easy to isolate and mainly composed by the over-expressed protein facilitates product recovery at industrial scale [20, 21]. Thus, the production of recombinant proteins as inclusion bodies represent a cost-effective alternative for enzyme manufacturing, provided that efficient large scale refolding methods are available. This is the case for calf prochymosin, the precursor of chymosin, which is widely used in cheese making and mostly obtained from recombinant microorganisms [18, 22].
It has been recently shown that proteins contained in inclusion bodies possess a degree of secondary structure and exhibit biological activity in many cases [23–28]. Moreover, the quality of the inclusion bodies, determined by the degree of folding of the aggregated proteins, depends on different factors; many of them related to the translation rate of the corresponding mRNA [23, 24]. In addition, it has been recently demonstrated that synonymous codon replacement is not always silent . Codon optimization affects translation rate which, in turn, may alter protein structure and function and thus the efficiency of refolding of proteins recovered from inclusion bodies [29–31]. Therefore, codon-optimized genes may lead to the formation of inclusion bodies from which the recovered proteins could be difficult to refold.
In this study, calf prochymosin was used as a model to perform a quantitative and qualitative evaluation of the inclusion bodies obtained from the expression of a set of synthetic genes encoding this protein. Seven genes were designed and synthesized using two different codon optimization strategies, and the amount of recombinant protein and the refolding yield of the inclusion bodies obtained with each gene were compared.
Gene design, synthesis and expression vector construction
Analysis of the native calf prochymosin gene revealed that almost 59% of codons are not the preferred for E. coli (Table 1). For example, seven out of the eight codons encoding for arginine found in the native sequence are represented with a frequency below 6% in the E. coli genome. Six of these codons are AGA or AGG, which have been shown to cause low levels of expression and mistranslational errors [1, 32]. Thus, it was surmised that codon optimization of the calf prochymosin gene might result in an increase in protein expression.
Seven variants of the calf prochymosin gene were designed using two different codon optimization strategies. In the first approach, only one codon was assigned for each amino acid to create two sequences named V0 and V1. For the V0 gene, the preferred codon found in the entire genome of E. coli W3110 was assigned to each amino acid. For the design of the V1 gene, a similar strategy was used; but in this case the favorite codon found in a set of highly expressed genes was employed to encode for each amino acid.
The second strategy used for codon optimization consisted on randomly assigning a triplet for each amino acid using a preference table http://www.kazusa.or.jp/codon/cgi-bin/showcodon.cgi?species=316407, with a probability based on the weight of each codon within the set encoding a given amino acid. Using this algorithm, five sequences were independently designed using the GeMS software package , and named V2-V6. The codon distribution for sequences V0, V1 and V2 is shown in Table 1 and the sequences of the seven genes and their codon usage is shown in Additional file 1.
The codon-optimized synthetic genes were created by using single strand 5´phosphorylated complementary primers. In all the cases 27 primers with a length ranging from 38 to 42 bases were used to create the leading strand and 27 primers with a length between from 38 to 43 bases were used to create the lagging strand. For all the genes the designed single stranded oligonucleotides overlapped each other by a minimum of 18 bases to ensure annealing. Two TAA stop codons in tandem were added at the end of each coding region followed by an Eco RI site. Additionally, an Nde I site overlapping the initial ATG was included in all the genes to mobilize the ORFs. Finally, all the synthetic genes were inserted into the expression vector pBru , where gene expression is driven by the PBAD promoter, inducible by the addition of L-arabinose.
Protein expression and yield
Recombinant plasmids were transformed into the E. coli W3110 strain for expression tests. Cell cultures were induced by the addition of L-arabinose when OD600 reached 0.5 and cells were harvested after 5 h. As previously reported, analysis of soluble and insoluble fractions of cell lysates by SDS-PAGE showed that all detectable prochymosin was located in the insoluble fraction in the form of inclusion bodies [18, 22]. The amount of prochymosin produced by the expression of each gene variant was quantified by densitometry of the stained gels and is shown in Figure 1 and Table 2. The variants V2-V6 optimized by the codon randomization approach produced significantly more proteins than the native sequence. The best result was obtained for the V2 sequence with an increase of 70% in the amount of prochymosin accumulated. On the other hand, no significant improvement in protein expression was observed for the V0 and V1 variants. In both cases the amount of prochymosin produced was similar to that obtained upon the expression of the wild type gene. Inclusion bodies were washed and isolated and the yield of recombinant protein was also measured by the Lowry method . Consistent with the results obtained from gel scanning quantification, the yields for V2 and the wild type sequence were 490 mg/L and 282 mg/L respectively. As reported by other authors, no correlation was found between the codon bias, measured using the codon adaptation index , and the quantity of recombinant protein produced by the codon-optimized genes [2, 35].
Cell growth rates were determined for E. coli cells harboring each expression plasmid. No significant difference was observed among the growth curves of the recombinant strains indicating that none of the genes had a toxic effect. As expected, the growth rate of the strain harboring the empty vector was higher after the addition of the inducer (data not shown).
It has been reported that the 5'coding region is particularly important in modulating translation initiation [2, 35–37]. Predicted mRNA secondary structures did not correlate strongly with expression when the sequences were analyzed using the online RNA folding program "mfold" http://mfold.rna.albany.edu/. In order to investigate whether the reason for the low expression of the V0 variant was due to the presence of a 5' local element, the first 10 codons of this gene were replaced by those of the V2 gene; the variant showing the greatest level of prochymosin production. Expression of the chimeric gene V2/0 showed a modest increase in prochymosin production (Table 2), indicating that downstream elements account for the lower production of prochymosin observed for the V0 variant.
In order to investigate whether a further increase in protein expression level could be obtained by using a stronger promoter, an additional experiment was carried out. For this, the V2 sequence was cloned into the pET 24b expression vector, where the expression is driven by the T7 promoter. The resulting plasmid was transformed into the BL21(DE3) strain and three colonies grown on LB medium. In all the cases a fall in the optical density with a concomitant increase in the viscosity of the cultures was observed within 6 h after the inoculation, indicating cell lysis. This observation suggests a toxic effect of the V2 gene when controlled by the strong T7 promoter.
In vitro refolding efficiency
Codon optimization affects translation rate which, in turn, may alter protein structure and function. It has been described that inclusion bodies formed in E. coli under different expression conditions may differ in their quality and, therefore, their ability to yield active proteins [25, 39, 40]. To the best of my knowledge, the impact of codon optimization on the quality of inclusion bodies has not been previously studied. Thus, I decided to investigate the ability of inclusion bodies obtained from the expression of different gene variants to yield functional prochymosin.
Inclusion bodies were prepared from cultures of E. coli W3110 strain harboring expression plasmids for the seven prochymosin gene variants described above. In all cases, inclusion bodies were washed after cell lysis and recovered as a white paste. The paste was dissolved in 8 M urea and the total protein concentration adjusted to 20 g/L. For all the preparations, PAGE analysis showed that more than 90% of the protein contained in the urea solution corresponds to calf prochymosin (data not shown). The urea solution was rapidly diluted into refolding buffer supplemented with 0.5 M arginine and 10 μM Cu++ since these additives were previously shown to increase the refolding efficiency of prochymosin .
Figure 2 shows the refolding efficiency obtained for the different inclusion body preparations. No significant differences were found in the recovery of calf prochymosin for inclusion bodies prepared using the different gene variants. In the refolding method employed, air oxidation was used to promote the formation of disulphide bonds and the activity recovered was similar to that previously reported for the native sequence [18, 41]
Codon optimization of the calf prochymosin gene was chosen due to the commercial value of improving its expression and to study the impact of codon optimization in an established production process. Even when the wild type sequence has been reported to express well in E. coli[18, 22], the presence of some rare codons led me to investigate codon optimization strategies in order to increase the expression of this protein. In the present study, seven gene variants were designed and synthesized to evaluate the effect of the two most common gene design strategies on the production of calf prochymosin in E. coli. The five sequences designed using the codon randomization strategy yielded higher protein quantities than those designed with the one amino acid-one codon method. These results suggest that the former method is a superior strategy for codon optimization. In addition, codon randomization permits flexibility in codon selection to facilitate gene design by avoiding: (i) repetitive elements that may lead to gene deletions, (ii) internal Shine-Dalgarno sequences, (iii) secondary mRNA structures and (iv) unwanted restriction sites. Some of these advantages of the codon randomization over the one amino acid-one codon method have been previously highlighted by Villalobos and coworkers . However, no studies have been conducted comparing these methods side by side and many authors still propose synthetic gene-based production improvements using the one amino acid-one codon method [4, 8, 11–15].
Sequences V2-V6 were designed based on a codon usage table obtained from the entire genome of E. coli W3110. The analysis of the codon distribution of these sequences shows that the differences in expression among these genes cannot be explained by the random assignation of rare codons when using the codon randomization method (see Additional file 2). Recently, Welch and co-workers have shown that most favorable codons are those read by tRNAs that are most highly charged during amino acid starvation rather than those that are most abundant in the genome . Using a codon table created based on these findings may provide genes to further increase calf prochymosin expression in E. coli.
No differences in protein production were found between V0 and V1 sequences, where the one amino acid-one codon algorithm using different codon tables was used for the design. This result suggests that, when employing this method, the translation efficiency may be limited by other constraints rather than the choice of the favorite codon to encode a given amino acid in the designed gene. The deleterious effect on gene expression of an imbalanced tRNA pool, previously proposed by several authors, is a likely explanation [4, 42, 43]. All the experiments described in this study were conducted using the PBAD promoter to drive the expression of the synthetic genes. Attempts to increase the level of recombinant protein using the stronger T7 promoter resulted in early lysis of the cells. A likely explanation for this observation is that the higher translation rate of the redesigned genes, associated with the leaky repression of the lac based T7 system in the abscense of inducer, may result in early accumulation of recombinant protein which prevents the healthy growth of the cultures. The finding of significant levels of prochymosin production in the absence of L-arabinose for E. coli harboring the V2 sequence supports this hypothesis.
Synonymous codon replacement may influence protein structure and function indicating that protein folding is DNA sequence dependent [24, 29, 44]. Polypeptides entrapped in E. coli inclusion bodies exhibit a variable degree of folding organization under different production conditions [19, 24]. Such degree of folding is frequently correlated with the "quality" of inclusion bodies because it has an impact on the yield of refolding and, therefore, the overall economy of the production process [45, 46]. This led me to explore the impact of using codon-optimized sequences on the ability of the resulting inclusion bodies to yield active chymosin. The refolding efficiency obtained for inclusion bodies recovered from recombinant clones expressing the seven individual variants were very similar, suggesting that the tested DNA sequences did not alter the conformational quality of protein contained in the inclusion bodies.
Variations in the rate of mRNA translation may influence the formation of secondary structures in the nascent polypeptide [47, 48], and analysis of gene sequences and the structure of their encoded proteins show that frequently used codons are associated with structural elements, while strings of less used codons tend to be present in boundaries separating such elements. Thus, a redesigned gene where most abundant codons are placed to encode secondary structures (like alpha helices) and rare codons are placed to encode linkers, may lead to the formation of inclusion bodies of superior quality. In the redesigned genes tested in this work, codon assignment frequency was equally distributed all along the entire gene. A calf prochymosin gene designed taking into account the sequence/structure relationship may provide insights into the influence of codon optimization on the refolding efficiency. This work is currently in progress in our laboratory.
Two alternative strategies for codon optimization have been evaluated in E. coli using calf prochymosin as a model. In all the cases the sequences created using the codon randomization method provided significantly more protein than their counterparts designed with one amino acid-one codon strategy, suggesting that this is a superior method for codon optimization. One of the obtained sequences produced more than 70% prochymosin than the native sequence, showing the potential of the approach to considerably reduce the production cost of well established production processes like in the case of chymosin.
Enzymes were obtained from New England Biolabs (USA) and used as recommended. DH5α, BL21(DE3) and W3110 E. coli strains were made chemically competent with a kit from Zymo Research (USA). Oligonucleotides were from Operon (USA). NTPs were PCR-grade from Roche Applied Sciences. DNA sequencing was performed on an ABI 3730 DNA analyzer (Applied Biosystems, USA) according to the manufacturer's recommended protocol. All other reagents were obtained from Sigma (USA)
Codon optimization, gene synthesis and cloning
Seven versions of the calf prochymosin A gene were designed and constructed. The variant V0 was designed using the software GeMS  and a codon table containing only the most abundant codon found in the entire genome of E. coli W3110 for each amino acid. The variant V1 was designed using the one amino acid-one codon algorithm from the Optimizer software . In this case, the favorite codon found in a set of highly expressed genes was used to encode each amino acid. Variants V2-V6 were designed using a codon randomization algorithm with the GeMS software and a codon table containing a fractional preference for each codon equal to that found in the genome of E. coli W3110. DNA sequences were synthesized using the method described by Reisinger and co-workers , digested with Nde I and Eco RI, inserted into the expression vector pBru  or pET 24b and verified by sequencing. In all the cases, E. coli DH5α was used for cloning.
Culture growth, calf prochymosin expression
E. coli W3110 cells harboring the expression vectors were grown with agitation in 1 l erlenmeyer flaks containing 100 ml of Luria-Bertani medium supplemented with glycerol (10 g/l) kanamycin (50 mg/l) at 30°C. Protein expression was induced when OD600 reached 0.5 units by adding 2 g/l L-arabinose and incubation continued for an additional 5 h period. Final OD600 was typically between 8-10 units. Cells were then harvested by centrifugation at 10,000g for 15 min at 4°C.
Protein analysis and in vitro refolding of calf prochymosin
Harvested cells (1 g wet weight) were resuspended in 40 ml of Tris-HCl 50 mM (pH 8.0) and incubated at 37°C for 30 minutes in the presence of lysozyme (0.2 mg/ml final concentration). Then, the mixture was sonicated on ice for 5 min. with 5 s pulses. Total extracts and proteins in different fractions were separated by SDS-PAGE on 10% gels, stained with Sypro-red and quantified by densitometry using a Typhoon scanner and BSA as a standard.
The inclusion bodies were isolated from the lysates by centrifugation at 10,000 g for 20 min at 20°C, washed twice with 50 ml of 10 mM EDTA (pH 8.0), 0.5% (v/v) Triton X-100 and once with 20 mM KH2PO4 (pH 7.5). Washed inclusion bodies were dissolved in deionized 8 M urea in 50 mM KH2PO4 (pH 10.5), rendering a final protein concentration of 20 mg/ml. The resulting solution was incubated with agitation for 2 h at 30°C, and centrifuged at 10,000 g for 10 min at 20°C. The preparation contained more than 95% prochymosin, as judged by SDS-PAGE. Refolding was carried out by rapid dilution of 1 ml of unfolded protein solution (20 mg/ml) in 20 ml of 50 mM KH2PO4, 0.5 M arginine, 10 μM CuSO4 (pH 10.5). The refolding solution was incubated for 12 h at 4°C. The renatured prochymosin was acidified to pH 2.0 with 2 M HCl and incubated 15 min at 20°C. Finally, samples were brought to pH 6.3 by the addition of 1 N NaOH and chymosin acivity measured using a milk clotting assay as previously described using authentic calf prochymosin (Sigma) as standard .
Burgess-Brown NA, Sharma S, Sobott F, Loenarz C, Oppermann U, Gileadi O: Codon optimization can improve expression of human genes in Escherichia coli: A multi-gene study. Protein Expr Purif. 2008, 59 (1): 94-102. 10.1016/j.pep.2008.01.008.
Welch M, Govindarajan S, Ness JE, Villalobos A, Gurney A, Minshull J, Gustafsson C: Design parameters to control synthetic gene expression in Escherichia coli. PLoS One. 2009, 4 (9): e7002-10.1371/journal.pone.0007002.
Zhou Z, Schnake P, Xiao L, Lal AA: Enhanced expression of a recombinant malaria candidate vaccine in Escherichia coli by codon optimization. Protein Expr Purif. 2004, 34 (1): 87-94. 10.1016/j.pep.2003.11.006.
Villalobos A, Ness JE, Gustafsson C, Minshull J, Govindarajan S: Gene Designer: a synthetic biology tool for constructing artificial DNA segments. BMC Bioinformatics. 2006, 7: 285-10.1186/1471-2105-7-285.
Ikemura T: Correlation between the abundance of Escherichia coli transfer RNAs and the occurrence of the respective codons in its protein genes. J Mol Biol. 1981, 146 (1): 1-21. 10.1016/0022-2836(81)90363-6.
Gustafsson C, Govindarajan S, Minshull J: Codon bias and heterologous protein expression. Trends Biotechnol. 2004, 22 (7): 346-353. 10.1016/j.tibtech.2004.04.006.
Kodumal SJ, Patel KG, Reid R, Menzella HG, Welch M, Santi DV: Total synthesis of long DNA sequences: synthesis of a contiguous 32-kb polyketide synthase gene cluster. Proc Natl Acad Sci USA. 2004, 101 (44): 15573-15578. 10.1073/pnas.0406911101.
Wang X, Li X, Zhang Z, Shen X, Zhong F: Codon optimization enhances secretory expression of Pseudomonas aeruginosa exotoxin A in E. coli. Protein Expr Purif. 2010, 72 (1): 101-106. 10.1016/j.pep.2010.02.011.
Menzella HG, Carney JR, Santi DV: Rational design and assembly of synthetic trimodular polyketide synthases. Chem Biol. 2007, 14 (2): 143-151. 10.1016/j.chembiol.2006.12.002.
Menzella HG, Reisinger SJ, Welch M, Kealey JT, Kennedy J, Reid R, Tran CQ, Santi DV: Redesign, synthesis and functional expression of the 6-deoxyerythronolide B polyketide synthase gene cluster. J Ind Microbiol Biotechnol. 2006, 33 (1): 22-28. 10.1007/s10295-005-0038-3.
Marlatt NM, Spratt DE, Shaw GS: Codon optimization for enhanced Escherichia coli expression of human S100A11 and S100A1 proteins. Protein Expr Purif. 2010, 73 (1): 58-64. 10.1016/j.pep.2010.03.015.
Zhen Feng LZ, Han Xue, Zhang Yanhe: Codon optimization of the calf prochymosin gene and its expression in Kluyveromyces lactis. World Journal of Microbiology and Biotechnology. 2010, 26 (5): 895-901. 10.1007/s11274-009-0249-2.
Gao W, Rzewski A, Sun H, Robbins PD, Gambotto A: UpGene: Application of a web-based DNA codon optimization algorithm. Biotechnol Prog. 2004, 20 (2): 443-448. 10.1021/bp0300467.
Fuglsang A: Codon optimizer: a freeware tool for codon optimization. Protein Expr Purif. 2003, 31 (2): 247-249. 10.1016/S1046-5928(03)00213-4.
Supek F, Vlahovicek K: INCA: synonymous codon usage analysis and clustering by means of self-organizing map. Bioinformatics. 2004, 20 (14): 2329-2330. 10.1093/bioinformatics/bth238.
Jayaraj S, Reid R, Santi DV: GeMS: an advanced software package for designing synthetic genes. Nucleic Acids Res. 2005, 33 (9): 3011-3016. 10.1093/nar/gki614.
Menzella HG, Reid R, Carney JR, Chandran SS, Reisinger SJ, Patel KG, Hopwood DA, Santi DV: Combinatorial polyketide biosynthesis by de novo design and rearrangement of modular polyketide synthase genes. Nat Biotechnol. 2005, 23 (9): 1171-1176. 10.1038/nbt1128.
Menzella HG, Gramajo HC, Ceccarelli EA: High recovery of prochymosin from inclusion bodies using controlled air oxidation. Protein Expr Purif. 2002, 25 (2): 248-255. 10.1016/S1046-5928(02)00006-2.
Parrilli E, Giuliani M, Marino G, Tutino ML: Influence of production process design on inclusion bodies protein: the case of an Antarctic flavohemoglobin. Microb Cell Fact. 2010, 9: 19-10.1186/1475-2859-9-19.
Burgess RR: Refolding solubilized inclusion body proteins. Methods Enzymol. 2009, 463: 259-282. full_text.
Singh SM, Panda AK: Solubilization and refolding of bacterial inclusion body proteins. J Biosci Bioeng. 2005, 99 (4): 303-310. 10.1263/jbb.99.303.
Wei C, Tang B, Zhang Y, Yang K: Oxidative refolding of recombinant prochymosin. Biochem J. 1999, 340 (Pt 1): 345-351. 10.1042/0264-6021:3400345.
Martinez-Alonso M, Gonzalez-Montalban N, Garcia-Fruitos E, Villaverde A: Learning about protein solubility from bacterial inclusion bodies. Microb Cell Fact. 2009, 8: 4-10.1186/1475-2859-8-4.
Martinez-Alonso M, Garcia-Fruitos E, Villaverde A: Yield, solubility and conformational quality of soluble proteins are not simultaneously favored in recombinant Escherichia coli. Biotechnol Bioeng. 2008, 101 (6): 1353-1358. 10.1002/bit.21996.
Ventura S, Villaverde A: Protein quality in bacterial inclusion bodies. Trends Biotechnol. 2006, 24 (4): 179-185. 10.1016/j.tibtech.2006.02.007.
Garcia-Fruitos E: Inclusion bodies: a new concept. Microb Cell Fact. 2010, 9: 80-10.1186/1475-2859-9-80.
Rodriguez-Carmona E, Cano-Garrido O, Seras-Franzoso J, Villaverde A, Garcia-Fruitos E: Isolation of cell-free bacterial inclusion bodies. Microb Cell Fact. 2010, 9: 71-10.1186/1475-2859-9-71.
Peternel S, Komel R: Isolation of biologically active nanomaterial (inclusion bodies) from bacterial cells. Microb Cell Fact. 2010, 9: 66-10.1186/1475-2859-9-66.
Kimchi-Sarfaty C, Oh JM, Kim IW, Sauna ZE, Calcagno AM, Ambudkar SV, Gottesman MM: A "silent" polymorphism in the MDR1 gene changes substrate specificity. Science. 2007, 315 (5811): 525-528. 10.1126/science.1135308.
Angov E, Hillier CJ, Kincaid RL, Lyon JA: Heterologous protein expression is enhanced by harmonizing the codon usage frequencies of the target gene with those of the expression host. PLoS One. 2008, 3 (5): e2189-10.1371/journal.pone.0002189.
Marin M: Folding at the rhythm of the rare codon beat. Biotechnol J. 2008, 3 (8): 1047-1057. 10.1002/biot.200800089.
Calderone TL, Stevens RD, Oas TG: High-level misincorporation of lysine for arginine at AGA codons in a fusion protein expressed in Escherichia coli. J Mol Biol. 1996, 262 (4): 407-412. 10.1006/jmbi.1996.0524.
Lowry OH, Rosebrough NJ, Farr AL, Randall RJ: Protein measurement with the Folin phenol reagent. J Biol Chem. 1951, 193 (1): 265-275.
Carbone A, Zinovyev A, Kepes F: Codon adaptation index as a measure of dominating codon bias. Bioinformatics. 2003, 19 (16): 2005-2015. 10.1093/bioinformatics/btg272.
Welch M, Villalobos A, Gustafsson C, Minshull J: You're one in a googol: optimizing genes for protein expression. J R Soc Interface. 2009, 6 (Suppl 4): S467-476. 10.1098/rsif.2008.0520.focus.
Kudla G, Murray AW, Tollervey D, Plotkin JB: Coding-sequence determinants of gene expression in Escherichia coli. Science. 2009, 324 (5924): 255-258. 10.1126/science.1170160.
Stenstrom CM, Isaksson LA: Influences on translation initiation and early elongation by the messenger RNA region flanking the initiation codon at the 3' side. Gene. 2002, 288 (1-2): 1-8. 10.1016/S0378-1119(02)00501-2.
Zuker M: Mfold web server for nucleic acid folding and hybridization prediction. Nucleic Acids Res. 2003, 31 (13): 3406-3415. 10.1093/nar/gkg595.
Martinez-Alonso M, Gonzalez-Montalban N, Garcia-Fruitos E, Villaverde A: The Functional quality of soluble recombinant polypeptides produced in Escherichia coli is defined by a wide conformational spectrum. Appl Environ Microbiol. 2008, 74 (23): 7431-7433. 10.1128/AEM.01446-08.
Rosano GL, Ceccarelli EA: Rare codon content affects the solubility of recombinant proteins in a codon bias-adjusted Escherichia coli strain. Microb Cell Fact. 2009, 8: 41-10.1186/1475-2859-8-41.
Tichy PJ, Kapralek F, Jecmen P: Improved procedure for a high-yield recovery of enzymatically active recombinant calf chymosin from Escherichia coli inclusion bodies. Protein Expr Purif. 1993, 4 (1): 59-63. 10.1006/prep.1993.1009.
Puigbo P, Guzman E, Romeu A, Garcia-Vallve S: OPTIMIZER: a web server for optimizing the codon usage of DNA sequences. Nucleic Acids Res. 2007, W126-131. 10.1093/nar/gkm219. 35 Web Server
Kurland C, Gallant J: Errors of heterologous protein expression. Curr Opin Biotechnol. 1996, 7 (5): 489-493. 10.1016/S0958-1669(96)80050-4.
Komar AA, Lesnik T, Reiss C: Synonymous codon substitutions affect ribosome traffic and protein folding during in vitro translation. FEBS Lett. 1999, 462 (3): 387-391. 10.1016/S0014-5793(99)01566-5.
Margreiter G, Messner P, Caldwell KD, Bayer K: Size characterization of inclusion bodies by sedimentation field-flow fractionation. J Biotechnol. 2008, 138 (3-4): 67-73. 10.1016/j.jbiotec.2008.07.1995.
Kischnick S, Weber B, Verdino P, Keller W, Sanders EA, Anspach FB, Fiebig H, Cromwell O, Suck R: Bacterial fermentation of recombinant major wasp allergen Antigen 5 using oxygen limiting growth conditions improves yield and quality of inclusion bodies. Protein Expr Purif. 2006, 47 (2): 621-628. 10.1016/j.pep.2006.01.009.
Purvis IJ, Bettany AJ, Santiago TC, Coggins JR, Duncan K, Eason R, Brown AJ: The efficiency of folding of some proteins is increased by controlled rates of translation in vivo. A hypothesis. J Mol Biol. 1987, 193 (2): 413-417. 10.1016/0022-2836(87)90230-0.
Marin M: Folding at the rhythm of the rare codon beat. Biotechnology Juornal. 2008, 3 (8): 1047-1057. 10.1002/biot.200800089.
Reisinger SJ, Patel KG, Santi DV: Total synthesis of multi-kilobase DNA sequences from oligonucleotides. Nat Protoc. 2006, 1 (6): 2596-2603. 10.1038/nprot.2006.426.
The author is the inventor of a patent application that includes part of the work described in this paper.