"

From 2011.igem.org

---

Software - Joint Codon Optimisation Algorithm

We wanted the option of expressing the genes responsible for auxin production in both B.subtilis and E.coli. To achieve this, we decided to joint codon optimise the IaaM and IaaH coding sequences. We could not find any software to do this so wrote our own.

---

1. Background

The genetic code is redundant which means that multiple codons can encode the same amino acid. Synonymous codons are circumstantially decoded at different speeds. Codon usage also varies between different species. This phenomenon means that sequence optimisation can be used to tune protein expression levels. It is tempting to think that one could codon optimise a sequence by selectively using an organism’s preferred codons. This is commonly referred to as the "one amino acid-one codon" method. Unfortunately it does not work. Recent optimisation studies^[1] have highlighted the importance of maintaining a diverse codon population in a given coding sequence. This said, the inclusion of ‘rare codons’ has also been shown to dramatically reduce protein expression E.coli.

2. Solution

In our approach to codon optimisation, we attempted to maintain codon diversity while simultaneously limiting rare codon inclusion. This was achieved by weighting codon selection using bias tables obtained from the Codon Usage Database^[2]. Joint optimisation was facilitated by combining the bias tables of E.coli and B.subtilis. Following the generation of a seed-sequence, the stochastic pruning of rare codons was used to iteratively optimise the sequence.

3. In-silico Testing

To test our codon optimisation algorithm, we ran the protein coding sequence for Dendra2 (BBa_K515007) through our software. The resultant DNA sequences were then fed into Genscript’s Codon Adaptation Index (CAI) analyser. This online tool measures the suitability of a sequence for expression in E.coli. Genscript’s own codon optimisation claims to be able to generate sequences with a CAI > 0.8. We were able to match this.

4. Future Work

Recent studies^[3] have suggested that rather than favouring highly used codons, it is preferable to select a subset for which the accompanying tRNAs are the most frequently charged during amino acid starvation. Once these codons have been tabulated for both E.coli and B.subtilis this data could be incorporated into our program.

5. R code

6. References

[1] Menzella, H.G. Comparison of two codon optimization strategies to enhance recombinant protein production in Escherichia coli. Microb Cell Fact 10, 15 (2011).
[2] Nakamura, Y., Gojobori, T. & Ikemura, T. Codon usage tabulated from international DNA sequence databases: status for the year 2000. Nucleic Acids Res 28, 292 (2000).
[3] Welch M, Govindarajan S, Ness JE, Villalobos A, Gurney A, et al. 2009 Design Parameters to Control Synthetic Gene Expression in Escherichia coli. PLoS ONE 4(9): e7002. doi:10.1371/journal.pone.0007002