Formal <i>vs</i>. experimental diversity.

<p>The number of randomized positions in each fragment library is indicated in brackets. The theoretical diversity represents the number of all possible combinations of 19 different codons by the given number of randomized positions. Fragment library B2(6) was used for the synthesis of L26; fragment library B2*(4)–for the synthesis of L24.</p><p>Formal <i>vs</i>. experimental diversity.</p

Summary of genetic complementation results.

<p>Amino acid substitutions in library members that complemented the HisF deficiency. The growth rate of <i>ΔhisF</i> cells, transformed with the corresponding plasmid variants was determined. The doubling time is a mean of five independent experiments.</p><p>Summary of genetic complementation results.</p

Observed average substitution rates.

<p>Codon exchange rates per randomized position in C- and B-fragment libraries and in the final gene libraries (L = L24+L26).</p><p>Observed average substitution rates.</p

A Robust and Versatile Method of Combinatorial Chemical Synthesis of Gene Libraries <i>via</i> Hierarchical Assembly of Partially Randomized Modules

<div><p>A major challenge in gene library generation is to guarantee a large functional size and diversity that significantly increases the chances of selecting different functional protein variants. The use of trinucleotides mixtures for controlled randomization results in superior library diversity and offers the ability to specify the type and distribution of the amino acids at each position. Here we describe the generation of a high diversity gene library using tHisF of the hyperthermophile <i>Thermotoga maritima</i> as a scaffold. Combining various rational criteria with contingency, we targeted 26 selected codons of the <i>thisF</i> gene sequence for randomization at a controlled level. We have developed a novel method of creating full-length gene libraries by combinatorial assembly of smaller sub-libraries. Full-length libraries of high diversity can easily be assembled on demand from smaller and much less diverse sub-libraries, which circumvent the notoriously troublesome long-term archivation and repeated proliferation of high diversity ensembles of phages or plasmids. We developed a generally applicable software tool for sequence analysis of mutated gene sequences that provides efficient assistance for analysis of library diversity. Finally, practical utility of the library was demonstrated in principle by assessment of the conformational stability of library members and isolating protein variants with HisF activity from it. Our approach integrates a number of features of nucleic acids synthetic chemistry, biochemistry and molecular genetics to a coherent, flexible and robust method of combinatorial gene synthesis.</p></div

Library analysis.

<p>(A) Exchange rate per type of resident amino acid. Comparison of the expected exchange probabilities (the fraction of foreign trimers in each trinucleotide mixture) with the observed exchange rate per type of resident codon (cumulative, sample size: 2762 exchanges). The number of resident amino acids among the 26 selected positions is indicated in parentheses, below is the respective codon. The differences in the expected exchange probabilities per type of residue are due to the biased composition of the trimer mix. The thick horizontal line indicates the average exchange probability (0.28), the dashed line–the average observed exchange rate (0.27). (B) Pertaining to replacing residues. Comparison of the expected frequency distribution of the amino acid exchanges (the molar ratio of each trimer in the trinucleotide mixture) with the observed frequency distribution for a sample of 2762 trinucleotide exchanges. (C) Binomial distribution of the number of mutations in L24 and L26. The observed distribution of the number of mutations per molecule (mutant class) is compared with the expected one (number of representatives in each mutant class for the given sample size = occurrences). The expected distribution of the mutant classes (L24: k = 1–24; L26: k = 1–26) is calculated according to the binomial distribution for a sample size of 239 sequences (L24) and 76 sequences (L26).</p

Assembly of the gene library.

<p>(A) Gel-eluted fragments from C-fragment libraries after restriction digestion (10% PAGE). The randomized fragments are marked in red; the brackets indicate the C-fragments, participating in one ligation assembly. Marker: 50 bp ladder (Thermo Scientific) (B) Ligation assembly of B1-fragment prior to gel elution. The full-length ligation product is indicated with an arrow. The length of the assembled B1 fragment with flanking BsaI restriction sites comprises 208 bp. Marker: 50 bp ladder (Thermo Scientific) (C) Gel-eluted fragments from B-fragment libraries. The length of the B1 fragment, liberated from the plasmid vector after restriction digestion comprises 182 bp. (D) Gene library ligation assembly. The full-length ligation product is indicated with an arrow. Marker: 1 kb ladder (Thermo Scientific).</p

Summary workflow diagram of the method for combinatorial gene synthesis.

<p><i>Library design</i>: the method imposes no restrictions to the number of residues selected for mutagenesis or to their location. The composition of the trinucleotide mixture and the exchange rate can be freely chosen without any limitations. Contact the oligonucleotide synthesis company of your choice well in advance and coordinate the purchase of trimer phosphoramidites. <i>Assembly design</i>: The fractionation of the gene sequence into modules depends on the locations of the residues for randomization. Divide the gene sequence into modules with a length between 40 bp and 90 bp, containing the mutagenized codons. The fragment borders of the modules should create unique overhands for ligation assembly after restriction digestion (see <a href="http://www.plosone.org/article/info:doi/10.1371/journal.pone.0136778#pone.0136778.t001" target="_blank">Table 1</a> and main text for details). Use PCR for longer stretches of wild type gene sequence. <i>Chemical synthesis</i>: Invest efforts and resources into high quality oligonucleotide synthesis as well as cloning and analysis of the C-libraries. The chemically synthesized diversity is stored in them and they are the starting point for all future gene libraries. Clone a wild type sequence, corresponding to each mutagenized module. <i>Library synthesis</i>: Generate B-fragment libraries by ligation from C-fragment libraries. Different B-fragment libraries (1, …, n) can be obtained in parallel by exchanging or mixing mutagenized modules with the corresponding wild type modules. Archive the B-libraries. Generate gene libraries by ligation assembly of B-libraries. Use the combinatorics approach for generation of multiple unique gene libraries.</p

Molar fractions of 19 codons in standard mixture and their derivation from observed occurrences on the surface of proteins and in catalytic sites of enzymes.

<p>The frequency of the residues in the trinucleotide mixture is calculated from the surface residue frequencies (mesophiles) [<a href="http://www.plosone.org/article/info:doi/10.1371/journal.pone.0136778#pone.0136778.ref038" target="_blank">38</a>], multiplied with a factor of 0.75, plus catalytic residue frequencies [<a href="http://www.plosone.org/article/info:doi/10.1371/journal.pone.0136778#pone.0136778.ref039" target="_blank">39</a>], multiplied with a factor of 0.25. The frequency for cystein was set to zero and all other values were normalized to 1.</p

Proteins produced from gene variants chosen at random from library L24.

<p>Soluble crude protein extracts of <i>E</i>.<i>coli</i> DH5α after 4 h gene expression were analysed by Western blotting using <i>anti-</i>His-tag antibody. C- negative control: cells, transformed with empty vector; WT: wild type <i>thisF (</i>MW_tHisF-His6 = 29 kDa); M–marker. Mutations in each clone are listed in <a href="http://www.plosone.org/article/info:doi/10.1371/journal.pone.0136778#pone.0136778.s002" target="_blank">S1 Table</a>.</p

Amino acid sequence of tHisF and randomized positions.

<p>(A) Nucleotide sequence of synthetic <i>thisF</i> gene. The protein comprises 253 amino acid residues, corresponding to 759 nucleotide of its structural gene. His6-tag was fused C-terminally to the coding sequence (depicted in blue together with stop codon). Red font indicates randomized positions. (B) Three dimensional location of randomized chain positions. The residues selected for randomization are presented as surfaces and colored according to their distance from the top (C-terminal ends of the β-strands). Blue: I52, A54, S55, S144, I173, D176, G177, K179; green: D11, L50, G80, H84, N103, T104, A105, D130, T171, S201, A204, A224, S225: orange: A8, S101, A128, L222, A223.</p