Difference between revisions of "Team:Heidelberg/Notebook"

Line 5: Line 5:
 
In our Notebooks you will find all information on how specifically we developed the different parts of our project. Have a look at our weekly reports below to find descriptions of what we worked on chronologically. Information on <a href="https://2017.igem.org/Team:Heidelberg/Materials">Materials</a> and <a href="https://2017.igem.org/Team:Heidelberg/Experiments">Methods</a> is provided on a separate pages. Finally you will find a detailled documentation of all our labwork in the <a href="https://2017.igem.org/Team:Heidelberg/Database">Notebook Database</a>.  
 
In our Notebooks you will find all information on how specifically we developed the different parts of our project. Have a look at our weekly reports below to find descriptions of what we worked on chronologically. Information on <a href="https://2017.igem.org/Team:Heidelberg/Materials">Materials</a> and <a href="https://2017.igem.org/Team:Heidelberg/Experiments">Methods</a> is provided on a separate pages. Finally you will find a detailled documentation of all our labwork in the <a href="https://2017.igem.org/Team:Heidelberg/Database">Notebook Database</a>.  
 
}}}}
 
}}}}
{{Heidelberg/templateus/Contentsection|{{#tag:html|{{Heidelberg/accpanel|#F8991D|{{Heidelberg/accord|Week  
+
{{Heidelberg/templateus/Contentsection|{{#tag:html|{{Heidelberg/accpanel|#F8991D|{{Heidelberg/accord|Calendar Week
 
27|{{#tag:html|
 
27|{{#tag:html|
 
<h2>Split_Polymerase</h2>
 
<h2>Split_Polymerase</h2>
Line 150: Line 150:
 
Carlson, J.C., Badran, A.H., Guggiana-Nilo, D.A., and Liu, D.R. (2014). Negative selection and stringency modulation in phage-assisted continuous evolution. Nat Chem Biol 10, 216-222.<br>Pu, J., Zinkus-Boltz, J., and Dickinson, B.C. (2017). Evolution of a split RNA polymerase as a versatile biosensor platform. Nat Chem Biol 13, 432-438.
 
Carlson, J.C., Badran, A.H., Guggiana-Nilo, D.A., and Liu, D.R. (2014). Negative selection and stringency modulation in phage-assisted continuous evolution. Nat Chem Biol 10, 216-222.<br>Pu, J., Zinkus-Boltz, J., and Dickinson, B.C. (2017). Evolution of a split RNA polymerase as a versatile biosensor platform. Nat Chem Biol 13, 432-438.
 
}}|week_id-0}}
 
}}|week_id-0}}
{{Heidelberg/accord|Week  
+
{{Heidelberg/accord|Calendar Week
 
28|{{#tag:html|
 
28|{{#tag:html|
 
<h2>Split_Polymerase</h2>
 
<h2>Split_Polymerase</h2>
Line 503: Line 503:
 
Esvelt, K.M., Carlson, J.C., and Liu, D.R. (2011). A system for the continuous directed evolution of biomolecules. Nature 472, 499-503.
 
Esvelt, K.M., Carlson, J.C., and Liu, D.R. (2011). A system for the continuous directed evolution of biomolecules. Nature 472, 499-503.
 
}}|week_id-1}}
 
}}|week_id-1}}
{{Heidelberg/accord|Week  
+
{{Heidelberg/accord|Calendar Week
 
30|{{#tag:html|
 
30|{{#tag:html|
 
<h2>Modeling</h2>
 
<h2>Modeling</h2>
Line 681: Line 681:
 
</ul>
 
</ul>
 
}}|week_id-2}}
 
}}|week_id-2}}
{{Heidelberg/accord|Week  
+
{{Heidelberg/accord|Calendar Week
 
31|{{#tag:html|
 
31|{{#tag:html|
 
<h2>Modeling</h2>
 
<h2>Modeling</h2>
Line 1,107: Line 1,107:
 
For this reason, a set of test reactions were performed with the commercially available Dimethylphenylsilan as a replacement of the 4-(dimethylsilyl)-anilin, always requiring custom synthesis. Hereby, we referred to the procedure of >“Lysate-small scale biocatalytic reactions“, mainly differing from the previous experiment by the addition of Dithionite solution as a reducing agent and GC Septum vials as reaction vials. Besides four samples with purified cyt c under various conditions, two samples still contained the whole cell solution.
 
For this reason, a set of test reactions were performed with the commercially available Dimethylphenylsilan as a replacement of the 4-(dimethylsilyl)-anilin, always requiring custom synthesis. Hereby, we referred to the procedure of >“Lysate-small scale biocatalytic reactions“, mainly differing from the previous experiment by the addition of Dithionite solution as a reducing agent and GC Septum vials as reaction vials. Besides four samples with purified cyt c under various conditions, two samples still contained the whole cell solution.
 
}}|week_id-3}}
 
}}|week_id-3}}
{{Heidelberg/accord|Week  
+
{{Heidelberg/accord|Calendar Week
 
32|{{#tag:html|
 
32|{{#tag:html|
 
<h2>Modeling</h2>
 
<h2>Modeling</h2>
Line 1,369: Line 1,369:
 
The Golden Gate product is purified using the ZymoResearch Kit and then 1 µl is<br>electroporated into 25 µl of S1059 E. coli. After 4 hours incubation in 3 ml SOC medium the cell culture is used for an plaque assay. A plaque PCR showed that the Tet(X)-insert is not included in the Golden Gate product and<br>therefore new Golden Gate overhangs are designed.
 
The Golden Gate product is purified using the ZymoResearch Kit and then 1 µl is<br>electroporated into 25 µl of S1059 E. coli. After 4 hours incubation in 3 ml SOC medium the cell culture is used for an plaque assay. A plaque PCR showed that the Tet(X)-insert is not included in the Golden Gate product and<br>therefore new Golden Gate overhangs are designed.
 
}}|week_id-4}}
 
}}|week_id-4}}
{{Heidelberg/accord|Week  
+
{{Heidelberg/accord|Calendar Week
 
33|{{#tag:html|
 
33|{{#tag:html|
 
<h2>Cytochrome_Engineering</h2>
 
<h2>Cytochrome_Engineering</h2>
Line 2,090: Line 2,090:
 
PCR proofed that the insertes are present in AP_destructase and a glycerol stock is created.
 
PCR proofed that the insertes are present in AP_destructase and a glycerol stock is created.
 
}}|week_id-5}}
 
}}|week_id-5}}
{{Heidelberg/accord|Week  
+
{{Heidelberg/accord|Calendar Week
 
34|{{#tag:html|
 
34|{{#tag:html|
 
<h2>Cytochrome_Engineering</h2>
 
<h2>Cytochrome_Engineering</h2>
Line 2,442: Line 2,442:
 
To test under which conditions the selection phage propagates best, 12 variations of the destrucatase accessory plasmid are assembled by gibson assembly.<br>The variations include different RBS and origins of replication. Because time was running out we decided to end this project at that point and concentrated on other PACE and PREDCEL experiments.
 
To test under which conditions the selection phage propagates best, 12 variations of the destrucatase accessory plasmid are assembled by gibson assembly.<br>The variations include different RBS and origins of replication. Because time was running out we decided to end this project at that point and concentrated on other PACE and PREDCEL experiments.
 
}}|week_id-6}}
 
}}|week_id-6}}
{{Heidelberg/accord|Week  
+
{{Heidelberg/accord|Calendar Week
 
35|{{#tag:html|
 
35|{{#tag:html|
 
<h2>Cytochrome_Engineering</h2>
 
<h2>Cytochrome_Engineering</h2>
Line 2,567: Line 2,567:
 
In analogy to the Bt toxin paper, it should be possible to implement a sequencing workflow, enabling a reliable basecalling accuracy. A scheme of this workflow is provided as a picture.
 
In analogy to the Bt toxin paper, it should be possible to implement a sequencing workflow, enabling a reliable basecalling accuracy. A scheme of this workflow is provided as a picture.
 
}}|week_id-7}}
 
}}|week_id-7}}
{{Heidelberg/accord|Week  
+
{{Heidelberg/accord|Calendar Week
 
36|{{#tag:html|
 
36|{{#tag:html|
 
<h2>Cytochrome_Engineering</h2>
 
<h2>Cytochrome_Engineering</h2>
Line 2,823: Line 2,823:
 
The embedding script did not really support restoring, that was fixed. The first embedding on uniprot that reached the second epoch was trained. Finding the best hyperparameters, especially the embedding dimension, and the length of the k-mers the embedding is based on, still remains. Checkpoints are saved as picklefiles now, making them accessible for applications without tensorflow.
 
The embedding script did not really support restoring, that was fixed. The first embedding on uniprot that reached the second epoch was trained. Finding the best hyperparameters, especially the embedding dimension, and the length of the k-mers the embedding is based on, still remains. Checkpoints are saved as picklefiles now, making them accessible for applications without tensorflow.
 
}}|week_id-8}}
 
}}|week_id-8}}
{{Heidelberg/accord|Week  
+
{{Heidelberg/accord|Calendar Week
 
37|{{#tag:html|
 
37|{{#tag:html|
 
<h2>Cytochrome_Engineering</h2>
 
<h2>Cytochrome_Engineering</h2>
Line 2,908: Line 2,908:
 
Training of the embedding on Uniprot continues. First approaches to optimize the current multilabel classifier by adding parameters to it, either in depth or in width started.<br>The concept of a generator that models natural evolution was designed. There, a starting sequence is exchanged by randomly chosen point mutations, which are then scored by the classifier with a sigmoid output for each label. Sequences that score higher are more likely to be chosen to continue working with. Those will then be changed by point mutations again. In this scenario the selection is modeled by a function that uses the classifier. However it remains unclear, if the classifier supports gradual improvements in a direction that improves the function in reality and in the classifier.
 
Training of the embedding on Uniprot continues. First approaches to optimize the current multilabel classifier by adding parameters to it, either in depth or in width started.<br>The concept of a generator that models natural evolution was designed. There, a starting sequence is exchanged by randomly chosen point mutations, which are then scored by the classifier with a sigmoid output for each label. Sequences that score higher are more likely to be chosen to continue working with. Those will then be changed by point mutations again. In this scenario the selection is modeled by a function that uses the classifier. However it remains unclear, if the classifier supports gradual improvements in a direction that improves the function in reality and in the classifier.
 
}}|week_id-9}}
 
}}|week_id-9}}
{{Heidelberg/accord|Week  
+
{{Heidelberg/accord|Calendar Week
 
38|{{#tag:html|
 
38|{{#tag:html|
 
<h2>Cytochrome_Engineering</h2>
 
<h2>Cytochrome_Engineering</h2>
Line 3,060: Line 3,060:
 
The concept of a genetical artificial intelligent algorithm (GAIA)  that is based on the classifier was developed. Starting from a given sequence it maximizes given labels and minimizes other given labels. At the same time the garbage score of the sequence is minimized. The model starts with one sequence but always works with B sequences. Therefore the starting sequence is in silico mutated B times to yield B different sequences. Those are then classified by DeeProtein and the outputs are used to calculate a sequence score. A given number of best sequences is selected and again mutated in silico to yield B sequences. This cycle of mutation and selection repeats and optimizes the starting sequence.<br>Additionally it is possible to also use recombination after the mutation to combine good subsequences
 
The concept of a genetical artificial intelligent algorithm (GAIA)  that is based on the classifier was developed. Starting from a given sequence it maximizes given labels and minimizes other given labels. At the same time the garbage score of the sequence is minimized. The model starts with one sequence but always works with B sequences. Therefore the starting sequence is in silico mutated B times to yield B different sequences. Those are then classified by DeeProtein and the outputs are used to calculate a sequence score. A given number of best sequences is selected and again mutated in silico to yield B sequences. This cycle of mutation and selection repeats and optimizes the starting sequence.<br>Additionally it is possible to also use recombination after the mutation to combine good subsequences
 
}}|week_id-10}}
 
}}|week_id-10}}
{{Heidelberg/accord|Week  
+
{{Heidelberg/accord|Calendar Week
 
39|{{#tag:html|
 
39|{{#tag:html|
 
<h2>Cytochrome_Engineering</h2>
 
<h2>Cytochrome_Engineering</h2>
Line 3,082: Line 3,082:
 
Embedding on Uniprot is calculated over 15 epochs, vectors that are close to each other in the high dimensional space seem to be similar on sequence level on first sight.<br>Genetic artificial intelligent algorithm is implemented. Features are mutation of a given maximum of residues, mutation of specific residues. Those are inserted in the model in small letters while the rest of the sequence is written in capital letters. Furthermore weighting arbitrary numbers of goal- and avoid- GO-terms are possible. The variance of the given classes, a measure for the authenticity of a protein, is substracted from the overall score. Recombination of sequences was not implemented as it corrupts maintaining less then the maximum allowed mutations in a sequence.
 
Embedding on Uniprot is calculated over 15 epochs, vectors that are close to each other in the high dimensional space seem to be similar on sequence level on first sight.<br>Genetic artificial intelligent algorithm is implemented. Features are mutation of a given maximum of residues, mutation of specific residues. Those are inserted in the model in small letters while the rest of the sequence is written in capital letters. Furthermore weighting arbitrary numbers of goal- and avoid- GO-terms are possible. The variance of the given classes, a measure for the authenticity of a protein, is substracted from the overall score. Recombination of sequences was not implemented as it corrupts maintaining less then the maximum allowed mutations in a sequence.
 
}}|week_id-11}}
 
}}|week_id-11}}
{{Heidelberg/accord|Week  
+
{{Heidelberg/accord|Calendar Week
 
40|{{#tag:html|
 
40|{{#tag:html|
 
<h2>Cytochrome_Engineering</h2>
 
<h2>Cytochrome_Engineering</h2>
Line 3,105: Line 3,105:
 
GAIA was expanded by allowing multiple mutations before selection takes place, a linear decay of the number of mutations before a round of selection is applied. Results look slightly better compared to one mutation before selection. With this mode epistatic mutations can be found, that means mutations that are beneficial together but decrease the score of a protein when occuring alone.
 
GAIA was expanded by allowing multiple mutations before selection takes place, a linear decay of the number of mutations before a round of selection is applied. Results look slightly better compared to one mutation before selection. With this mode epistatic mutations can be found, that means mutations that are beneficial together but decrease the score of a protein when occuring alone.
 
}}|week_id-12}}
 
}}|week_id-12}}
{{Heidelberg/accord|Week  
+
{{Heidelberg/accord|Calendar Week
 
41|{{#tag:html|
 
41|{{#tag:html|
 
<h2>Cytochrome_Engineering</h2>
 
<h2>Cytochrome_Engineering</h2>
Line 3,280: Line 3,280:
 
GAIA was used to generate mutations in a glucoronidase sequence, beneficial for galactosidase activity. A subset of the generated sequences was ordered as oligonucleotides for kinetic assays. In order to show that the score from the used classifier correlates with the real-world activity, a set of less active beta-lactamase sequences was predicted by gaia. Those were also bought as oligonucleotides and will be tested for their activity.<br>A randomization function was implemented to examine how the different scores develop when the number of mutations increases. The results look sigmoid, showing almost no change in the score for the first mutations, then rapidly dropping below the threshold for positive classification, after which they only decrease slowly.<br>A combination function to examine all possible combinations of a set of mutations was implemented in order to aid the decisions on which sequences to test in the lab.
 
GAIA was used to generate mutations in a glucoronidase sequence, beneficial for galactosidase activity. A subset of the generated sequences was ordered as oligonucleotides for kinetic assays. In order to show that the score from the used classifier correlates with the real-world activity, a set of less active beta-lactamase sequences was predicted by gaia. Those were also bought as oligonucleotides and will be tested for their activity.<br>A randomization function was implemented to examine how the different scores develop when the number of mutations increases. The results look sigmoid, showing almost no change in the score for the first mutations, then rapidly dropping below the threshold for positive classification, after which they only decrease slowly.<br>A combination function to examine all possible combinations of a set of mutations was implemented in order to aid the decisions on which sequences to test in the lab.
 
}}|week_id-13}}
 
}}|week_id-13}}
{{Heidelberg/accord|Week  
+
{{Heidelberg/accord|Calendar Week
 
42|{{#tag:html|
 
42|{{#tag:html|
 
<h2>Modeling</h2>
 
<h2>Modeling</h2>
Line 3,412: Line 3,412:
 
The infected overnight culture was centrifuged at 8000 xg for 3 min and phage supernatant was transferred into a fresh tube stored at 4 °C.
 
The infected overnight culture was centrifuged at 8000 xg for 3 min and phage supernatant was transferred into a fresh tube stored at 4 °C.
 
<h2>Protein_Interaction</h2>
 
<h2>Protein_Interaction</h2>
<h2 id="-week-42-"><strong>Week 42:</strong></h1>
+
<h2 id="-week-42-"><strong>Week 42:</strong></h2>
 
By plaque assay performance phage propagation on Jin69 was proved to work, so gene circuit can be used in PACE and/or PREDCEL. After several attempts of positive AP sequencing AP were finally proven right. Related stocks were proven not to be contaminated by contamination detection PCR. <em>S1030 Cells</em> containing approved APs were made electrocompetent and MP1 and MP4 were transformed into these cells separately by electroporation. For proper and complete sequencing of T7 RNAP to detect mutations new sequencing primers were designed and ordered. Furthermore to exchange SD8 RBS sequence with weaker sd8 RBS sequence GoldenGate primers were bought and used to perform GoldenGate cloning.
 
By plaque assay performance phage propagation on Jin69 was proved to work, so gene circuit can be used in PACE and/or PREDCEL. After several attempts of positive AP sequencing AP were finally proven right. Related stocks were proven not to be contaminated by contamination detection PCR. <em>S1030 Cells</em> containing approved APs were made electrocompetent and MP1 and MP4 were transformed into these cells separately by electroporation. For proper and complete sequencing of T7 RNAP to detect mutations new sequencing primers were designed and ordered. Furthermore to exchange SD8 RBS sequence with weaker sd8 RBS sequence GoldenGate primers were bought and used to perform GoldenGate cloning.
 
<h2>Software</h2>
 
<h2>Software</h2>
Line 3,418: Line 3,418:
 
For further insight in properties of a sequence, a function was implemented, that scores every possible single mutant and plots changes in the score. Potentially active site and other crucial sequence regions can be predicted by this.
 
For further insight in properties of a sequence, a function was implemented, that scores every possible single mutant and plots changes in the score. Potentially active site and other crucial sequence regions can be predicted by this.
 
}}|week_id-14}}
 
}}|week_id-14}}
{{Heidelberg/accord|Week  
+
{{Heidelberg/accord|Calendar Week
 
43|{{#tag:html|
 
43|{{#tag:html|
 
<h2>Optogenetics</h2>
 
<h2>Optogenetics</h2>

Revision as of 23:06, 1 November 2017

Notebook
Weekly Reports
In our Notebooks you will find all information on how specifically we developed the different parts of our project. Have a look at our weekly reports below to find descriptions of what we worked on chronologically. Information on Materials and Methods is provided on a separate pages. Finally you will find a detailled documentation of all our labwork in the Notebook Database.
Calendar Week 27

Split_Polymerase

Documentation week 07/03 - 07/09/2017a

Wednesday 07/05/2017

Transformation of electrocompetent cells for reproduction of Dickinson-Experiments

In order to gain the three strains which are required for the reproduction of the first seven days of the “Evolution of a split RNA polymerase as a versatile biosensor platform” (Pu et al., 2017), transformation of electrocompetent cells with different strains was performed. S1030 electrocompetent cells were thawed and DNA of the appropriate plasmids were added. The amount of DNA of each plasmid was adjusted at equivalent amounts. First, the recovery medium was cooled on ice and thawed. The electroporation cuvettes were also refrigerated on ice. Following this, the 40 µl of S1030 cells were added to the DNA in a separate tube and mixed carefully. Next, the mixture was transferred to the cuvettes and electroporation was performed using 2500 V. Cells were immediately recovered using 1 ml SOC medium and grown at 37 °C for 1 h and 275 rpm. The resulting cell culture was centrifuged at 6000 rpm for 2 min and the cell pellet was resuspended using 200 µl SOC medium. After this, 200 µl of the cell culture were applied on the related LB-Agar plates, displaying the corresponding antibiotic resistance and 10 mM Glucose concentration. At last, the plates were incubated at 37 °C over night. For further information about the detailed transformation: awesome_sheet/Trafos/ ID 66-69

Thursday 07/06/2017

Transformation of electrocompetent cells for reproduction of Dickinson-Experiments

Due to the non-functionality of the first transformations on the day before, the transformations were repeated. The general procedure was implemented following the instructions of Wednesday 07/05/2017. Several changes were made in order to complete the transformation successfully. All DNA amounts were adjusted at 100 ng to improve the efficiency. Beyond that, 25 mM Glucose instead of 10 mM Glucose were added to the LB-agar plates, following the instructions of the information supplied by addgene. For further information about the detailed transformation: awesome_sheet/Trafos/ ID 70 -73

Friday 07/07/2017

Transformation of electrocompetent cells for reproduction of Dickinson-Experiments

Since the second transformations did not work as well, several changes were made for enable a successfull transformation. The general protocol of Wednesday 07/05/2017 was implemented like before. Glucose concentration in LB-agar plates as well as in SOC medium was adjusted at 100 mM Glucose for prevention of the activation of the MP1 plasmid. In comparison, another transformation was implemented, using TOP10 chemo competent cells to investigate the effect of MP1 as a non-transformable plasmid. Furthermore, pSB1C3 plasmid was also used for transformation to test the general competence of the electrocompetent S1030 cells. A last, a transformation using only 1 ng of MP1 DNA was implemented under equivalent conditions with S1030 cells.

Reagent Stock concentration [mg/ml]
Ampicillin 50
Chloramphenicol 34
Kanamycin 10
Tetracycline 5
Streptomycin ???
Glucose 1.1 M

ID Plasmids Purification concentration [ng/µl]
1 MP1 83,8
41 pJin 70 N/A
60 pJin 178 28
62 2-22: pJC173b 29
63 pJin 69 (JM) 140
76 pJin 69 (MW) 373
77 pJin 177 167
79 pJin 172 514
80 pJin 173 232
124 pSB1C3 232
Plasmids with the IDs 41, 60, 62 and 63 were used incorrectly. The correct plasmids for further transformations have the following IDs: 76,77,79 and 80.

Preparation of Davis Rich Medium for transformation of MP1 in S1030

Davis Rich Medium (DRM) was used in all PACE papers for transformations as well as the PACE device itself. Therefore, using DRM should be of any great advantage against using normal media like LB, 2xYT or simple Esvelt medium. For preparation of DRM, the amounts of different chemicals was mixed according to the methods of the “Negative selection and stringency modulation in phage-assisted continuous evolution” (Carlson et al., 2014). The amounts were adjusted to 1 l medium instead of 20 l. 10 mM of 50 ml Borsäure and Coppersulfat were prepared using the appropriate amounts calculated using the molecular weights.

Reagent Molecular weight [g/mol]
Guanin 151,13
H3BO3 61,83
CuSO4 159,61

Gibson Assembly preparation

For gibson assembly, reaction mix was prepared according to the conditions of NEB. The overall reaction volume was 20 µl. 10 µl Gibson MM were mixed with 0.02 – 0.5 pmol of the DNA fragments. The reaction mixture was filled up with ddH2O. For calculation of the molecular weight of the fragments, NEBioCalculator was used. The respective length of the fragments was extracted from the awesome_sheet. If the insert contains under 200 bp, the respective amount was multiplicated by five. The reaction mix was incubated at 50 °C for 15 min (2-3 fragments) or for 60 min (4-6 fragments) and the resulting sample was stored at -20 °C. For the following transformation, 2 µl of the assembly mix were used. For further information, please visit Gibson Assembly - NEB.

Saturday 07/08/2017

Preparation of Davis Rich Medium-like medium for transformation of MP1 in S1030

Sunday 07/09/2017

Transformation of electrocompetent cells for reproduction of Dickinson-Experiments

References:

Carlson, J.C., Badran, A.H., Guggiana-Nilo, D.A., and Liu, D.R. (2014). Negative selection and stringency modulation in phage-assisted continuous evolution. Nat Chem Biol 10, 216-222.
Pu, J., Zinkus-Boltz, J., and Dickinson, B.C. (2017). Evolution of a split RNA polymerase as a versatile biosensor platform. Nat Chem Biol 13, 432-438.
Calendar Week 28

Split_Polymerase

Documentation week 07/10 – 07/16/2017

Monday 07/10/2017

Preparation of Davis Rich medium for triple and double transformation of Dickinson strains

Equation for calculation of the required mass for preparation of a 10 mM stock concentration. m = c x M x V

Reagent Molecular weight [g/mol]
Guanin 151,13
H3BO3 61,83
CuSO4 159,61
CoCl2 129.839
MnCl2 125.844
ZnSO4 161,47
(NH4)6Mo7 1235,86

Davis Rich Medium 20 liters 1 liter unit
K2HPO4 140 7 g
KH2PO4 40 2 g
(NH4)2SO4 20 1 g
NaCl 58 2,9 g
casamino acids 100 5 g
Tween-80 20ml 1ml
l-cysteine 1 0,05 g
l-tryptophan 0,5 0.025 g
adenine 0,5 0.025 g
guanine 0,5 0.025 g
cytosine 0,5 0.025 g
uracil 0,5 0.025 g
CaCl2 0.5 µM
After autoclaving
NaHCO3 16,8 0,84 g
glucose 90 4,5 g
sodium citrate 10 0,5 g
MgSO4 1 0,05 g
FeSO4 56 mg 2,8 mg
thiamine 134 mg 6,7 mg
calcium pantothenate 94 mg 4,7 mg
para-aminobenzoic acid 54 mg 2,7 mg
para-hydroxybenzoic acid 54 mg 2,7 mg
2,3-dihydroxybenzoic acid 62 mg 3,1 mg
(NH4)6Mo7 3 nM
H3BO3 400 nM
CoCl2 30 nM
CuSO4 10 nM
MnCl2 80 nM
ZnSO4 10 nM

Preparation of Esvelt medium for turbidostat study on the influence of casamino acids on the bacterial growth

Esvelt medium was prepared according to the instructions found in (Esvelt et al., 2011). Different amounts of casamino acids were added, evaluating the effect on the growth by OD measurement at different timepoints (one measurement per hour).
Several transformations were performed in order to gain the Dickinson plasmids or rather identify the problematic aspects of the triple transformation implemented in week 27. Therefore, MP1 was transformed in different amounts (1 ng and 100 ng) in electrocompetent S1030, as well as in chemocompetent Top10 cells. As a control, pSB1C3 was also transformed in S1030 for investigation of the general competence of the cells. 100 mM Glucose was added to each of the recovery media

Tuesday 07/11/2017

Restriction digest of MP1, MP4 and MP6

For digestion of MP1, as well as MP4 and MP6, 0,5 µl of the appropriate enzymes (MP1 – EcoRI; MP4 – BsmBI; MP6 – BsaI and SpeI) were added into 1 µl of the appropriate buffer (EcoRI – Cutsmart; BsmBI – 3.1; BsaI and SpeI – Cutsmart). 1 µl DNA was added to the mixture. All samples were filled up to 10 µl with ddH20 and digested for 30 min at 37 °C, or rather 55 °C for BsmBI.
Growth and preparation of MP1, MP4 and MP6
Agar stocks previously ordered from addgene were scratched out on LB-agar plates with 25 mM glucose and chloramphenicol resistance. Cells were grown on these plates overnight and picked for single colony growth in 2xYT-medium, containing 25 mM glucose and chloramphenicol resistance as well. Miniprep of plasmid DNA was performed using Qiagen Plasmid Mini Kit following the instructions of the attached protocol. All samples were eluted in ddH20 and concentration of the product was measured using the Implen spectrometer. For further information about the growth, please visit the awesome_sheet/growth.

PCR of Phage_backbone_b for PAM-PACE, fragment1 and fragment2 of destructase PACE and gIII for enzyme-PACE

PCR was performed in a total volume of 50 µl containing 25 µl 2x Phusion Flash Master mix, 2.5 µl of each 10 µM primer, as well as 1 µl of the appropriate template. 2,5 µl DMSO and 10 µl 5M Betaine were added to the mixture. PCR was performed under the following conditions: Initial denaturation at 98 °C for 10 secs, followed by 35 cycles of denaturation at 98 °C for 1 s, annealing at the respective temperature for 5 s, and primer extension at 72 °C for the associated time. The last extension step was carried out at 72 °C for 1 min. PCR product was evaluated by agarose gel electrophoresis. The annealing temperature for each primer pair was estimated using the Phusion Flash Tm calculator. For further information visit the awesome_sheet/PCRs ID35 - 41.

Wednesday 07/12/2017

PCR of Phage_backbone_b for PAM-PACE, fragment1 and fragment2 of destructase PACE and gIII for enzyme-PACE

PCR was performed in a total volume of 50 µl containing 25 µl 2x Phusion Flash Master mix, 2.5 µl of each 10 µM primer, as well as 1 µl of the appropriate template. 2,5 µl DMSO and 10 µl 5M Betaine were added to the mixture. PCR was performed under the following conditions: Initial denaturation at 98 °C for 10 secs, followed by 35 cycles of denaturation at 98 °C for 1 s, annealing at the respective temperature for 5 s, and primer extension at 72 °C for the associated time. The last extension step was carried out at 72 °C for 1 min. PCR product was evaluated by agarose gel electrophoresis. The annealing temperature for each primer pair was estimated using the Phusion Flash Tm calculator.
PCR Products 29 (65.9°C), 34 (65.4°C), 35 (66.2°C) and 36 (71.6°C) – awesome_sheet/PCRprod
Agarose gel electrophoresis
DNA fragments were separated using agarose gel electrophoresis. The gels were prepared by dissolving 1% agarose in 1 % TAE buffer. For visualization of the DNA, 5 µl Roti Safe gel stain was added. Gel electrophoresis was performed for 20 min at 140 V. For image acquisition, BioRad Gel imaging and documentation software was used. Size of the DNA was determined by comparison with the NEB 2-log ladder.
Gel extraction
Gel extraction has been conducted by using Monarch DNA Gel Extraction Kit executing the attached protocol. All samples were eluted in ddH20 and concentration of the product was measured using a nanodrop.

Thursday 07/13/2017

PCR of Phage_backbone_b for PAM-PACE, fragment1 and fragment2 of destructase PACE and gIII for enzyme-PACE + different gene III variants

PCR was performed according to the instructions of the previous day (Wednesday 07/12/2017). Several adjustments concerning the annealing time and the extension time were made, because of the non-functional PCR of the last day. In addition, three other PCRs were performed using the same reaction layout as the ones of the previous day.
For further information, please visit the awesome_sheet/PCRs, ID 88 - 108
Gibson Assembly of AP_Cas9_pSB1A3
For Gibson assembly reaction, 10 µl Gibson Assembly MasterMix is mixed with 0.2 – 1 pmol DNA of the different fragments. The reaction mix is filled up to 20 µl with ddH2O. The substance amount is calculated using NEBioCalculator with the length of the fragment and the concentration per µl. For assembly of the AP_Cas9_pSB1A3, five different fragments were mixed in equimolar concentrations. According to this, 15 µl of Gibson assembly MM, 13 µl of equimolar DNA fragments and 2 µl of ddH2O was mixed and incubated at 50 °C for 60 min.

Friday 07/14/2017

Preparation of media supplements for esvelt medium used in first PACE experiment

Substance Stock concentration [mg/ml] Volume [ml] Mass [g]
Ampicillin 50 250 12,5
Chloramphenicol 34 250 8,5
Kanamycine 10 250 2,5
All antibiotics were steril filtered using 0.45 µl steril filter and a syringe. The stock concentrations were all diluted 1:1000. Due to that, the respective volumes of the 250 ml stock concentrations were added to the media.

Supplements for esvelt medium:

Substance 1 liter 60 liters 80 liters 140 liters
Glucose 4.5 g 270 g 360 g 630 g
Sodium Citrat 0.5 g 30 g 40 g 70 g
L-Leucine 0.025 g 1.5 g 2 g 3.5 g
MgSO4 0.025 g 1.5 g 2 g 3.5 g
Casamino acids 5 g 300 g 400 g 700 g
The volume indicated in the heading marks the amount of ddH2O used for solvation of the substances.
One-liter arabinose solution was prepared in a concentration of 5 w/v (50 g/1 l).

Saturday 07/15/2017

Dotblot for identification of Phage contamination

Dotblot for investigation of the PACE run
2 µl of samples are spotted on the nitrocellulose membrane. A grid is drawn to indicate different samples. The membrane is drown and all non-specific sites are blocked by using 5 % BSA in TBS-Tween for 30 to 60 min. Afterwards, the membrane is incubated with the 1:200 dilution of M13 coat protein antibody which is directly bound to Alexa 488 for 30 min. The dilution was prepared using 2 ml BSA/TBS-Tween and 10 µl antibody. Following this, the membrane was washed three times with TBS-Tween for 5 min and imaged using a BioRad chemidoc with Image Lab 5.2.1.
PCRs for mNeonGreen amplification, Ampicillin resistance from pSB1A3 and
In order to optimize the mentioned PCRs listed above, the reaction mix was
Sunday 07/16/2017

Reference:

Esvelt, K.M., Carlson, J.C., and Liu, D.R. (2011). A system for the continuous directed evolution of biomolecules. Nature 472, 499-503.
Calendar Week 30

Modeling

The theoretical model for PALE is built, but not solved yet. Probably finding a solution symbolically is not possible. For that reason the next step is to write a python script that solves the system numerically.
A basic model is running in browser, written in javascript with numjys and plotly.js. Further models will be written the same way, old ones translated.

Software

DeeProtein

  • We implemented 1d convolutions to better capture the 1-dimensional character of protein sequences.
  • We generated a new test case to harness the accuracy issue for DeeProtein (the CNN was not able to grasp differences in the classes for the toplevel EC dataset). The new test case looks like this:

EC-number Function n(samples)
1.1.1.- With NAD(+) or NADP(+) as acceptor 2558
2.1.1.- Methyltransferases 3739
2.7.1.- Phosphotransferases with an alcohol group as acceptor 2137
4.2.1.- Hydro-lyases 3242
6.1.1.- Ligases forming aminoacyl-tRNA and related compounds reduced to 2913 (og: 11652)

Embedding

A first version of the rewritten generate batch function for a more memory efficient version of the embedding is ready. Testing and further improvements have to follow. This makes training an embedding on the full uniprot dataset possible.

Split_Polymerase

  • Monday: "Nachschicht" The psp-genIII AP of S1059 (pJC175e) strain was purified and transformed into electro-competent S2060 cells, in order to gain a better readout using plaque assays as phage detection method. S2060 has a beta-Galactosidase integrated in its genome, due to which blue plaques are resulting from a phage infection on soft agar with IPTG.
  • Tuesday: General concept of an addition to the awesome_sheet was prepared and adjusted. Further tests and progress report of Catharina are needed for further development.
    • In order to quantify a possible phage contamination, PCR using two different primer pairs were implemented. Master mix was prepared using OneTaq MM as well as primers JM_074, JM_082, JM_076 and JM_083. For further details about the detection PCR, please visit the PCR sheet of the awesome_sheet.
  • Wednesday: "Nachschicht" For quantification of phage contamination of the turbidostat, qPCR was performed following the standard protocol. Primers JM_074 and JM_082 were used for PCR reaction as well as ABI Sybr Green master mix.

1-3 4-6 7-9 10-12
1 H20 T34 T35 T36
2 T39 T41 L1.34 L1.35
3 L1.36 L1.39 L1.50 L2.34
4 L2.35 L2.36 L2.39 L2.41
5 L2.50 Pos 10-4 Pos 10-8 Pos 10-12
6 Neg unv. 2xYT LB Glu
The qPCR did not yield an result. In the following, qPCR was not referred to as a option for phage detection during the PACE runs.
As an alternative quantification method, DotBlots were tested. Due to too much noise during detection, they were not considered as detection method as well. For further information, please have a look at Catharina´s weekly report.
  • Thursday: PACE device was cleaned and autoclaved in order to change strains of split-polymerase PACE.
  • Friday: PACE device was build again and installed for usage (have a look at the provided video). Stock was thawn and a culture with S1030 + MP1 + Jin177 + Jin172 was grown. Inoculation of the turbidostat was performed and culture was grown until OD 0,6 - 0,8. Thereby, the last phase of our Dickinson-PACE experiment was started.
  • Saturday: PREDCEL - phage related dis-continous evolution concept was designed together with Tobias Stadelmann. Tobias tested the new experimental design by using a 50 ml S1030 + MP1 + Jin69 culture. Culture was prepared using LB-Medium, Cmp, Amp and 25 mM Glucose. The culture was grown until OD 0,6 - 0,8 and separated into nine 5 ml cultures (in 14 ml PP tubes). Afterwards, they were infected with different amount of phages (10 µl, 100 µl and 1 ml) to evaluate which of the amounts has to be used for following experiments. Beyond that, the infected cultures were centrifuged after different amounts of time (20, 30 or 60 min). The supernatant containing the phages was then transfered into a new 5 ml culture.

General concept of PREDCEL:

Idea: Using PACE in an experimental design, suitable to be carried out by every iGEM Team.
PREDCEL describes a manual implementation of PACE.
The incubation of E.coli with phages in a reaction tube (max. 14 ml) using a normal incubator should be performed.
First we wanted to test PREDCEL under as equivalent reaction conditions as possible (Dickinson-PREDCEL). PREDCEL allows us to test more “PACE-runs” manually and in parallel. It also enables to test different reaction conditions (different media, inducer, etc.) while having a
simple handling, suitable for every iGEM Team.

Pre-experiments:

  1. Test of conditions - amount of phage supernatant to be transferred and time until transfer between each "run"
    • Results: 1 ml phage supernatant should be transferred
      • should prevent washout of phages in experimental set-up
      • time independent - BUT to small reproduction size for a validated statement about the amount of time needed for propagation (20, 30, 60 min)
        • according to Julius phages need 10 min to infect and reproduce themselves for the first time
    • analysis was performed using plaque assays
  2. Test of medium which should be used for experiments
    • LB, 2YT and Esvelt medium
    • infectivity of the phages should be observed in dependence of the medium
      • analysis was performed using plaque assays
      • the three different times of experiment number 1. were adopted (20, 30, 60 min)
    • six rounds were implemented for 2YT and LB
    • three rounds were implemented for Esvelt medium

Experiments to go:

  1. Dickinson-PACE reproduction using the real phage (unevolved)
    • using new strains with Jin69_xxxx_with_LuxAB and MP1, Jin177_LuxAB + Jin172_LuxAB + MP1, Jin177_LuxAB + Jin173_LuxAB + MP1
    • Aim: measure LuxAB activity and sequencing of Dickinson phage for proof of concept
      • MP4 or MP6 activity should be tested as well
  2. First experiment of PAM-PACE functionality
    • does the PAM-PACE circuit work like it is supposed to

Prerequesite:

  • knowledge about time and medium
  • Sunday: Qiagen Midi Prep protocol for M13 phages was implemented for purification of our selection phage. Following this further, oxford nanopore standard protocol for testing the MinION sequencing was tested for feasibility. Several kits, supplied by NEB were ordered for preparation of the test library. A cooperation with AG Conrad of the Eils Labs should be started. They work with single cell sequencing and should therefore know a lot about library preparation and NGS.
Calendar Week 31

Modeling

We have a first set of equations for the PREDCEL Model, but did not solve it yet. Probably finding a solution symbolically is not possible. For that reason the next step is to write a python script that solves the system numerically.
Modelling the specialised variants did not start this week, as well as the overview.
However we have a basic model running in the browser, written in javascript with numjys and plotly.js. The next set of models should be written the same way, also the old ones have to be translated.

PACE

The current PACE setup was documented. Future PACE setups were discussed. Especially a setup that allows for host strains to be switched would be interesting for PamPACE. The idea is to use two turbidostats alternating. The host supply for the lagoons has a Y-connector that potentially connects both Turbidostats to the lagoons. Only one of the two tubings to the lagoons is allowed to be open at a time. When strains are switched, the block is switched. The Turbidostat that is not in use anymore, can be disconnected, the resulting open Luer-connector has to be sprayed with incidin, H2O2 or whatever kills phage and bacteria. Then the unused Turbidostat can be cleaned, autoclaved and reinstalled. It should be possible to switch strains every 12h theoretically. This setup needs either a second bottle of medium or changed lids, that hold two tubings that collect medium. New tubing to the turbidostat waste is needed.
This setup makes it possible to switch back and forth between two strains, that can live on different media.
The overall principle of a chemostat is that the population is kept constant by adjusting the amount of an essential resource. This helps to save medium but decreases the metabolic activity of the E. coli which may be problematic for PACE. But there are papers (Dickinson et al.) that used chemostats. Probably their protocol can be used with a few adjustments.

Phage_Propagation_and_Quantification

Phage Quick-Test

To find an easy and quick detection and quantification method for phages, the activity of beta-galactosidase, which is under control of a psp-promoter can be measured.
Therefore, 100 µl of the E. coli strain S2208, stock ID: 47 (OD600 ~ 0.5) was pipetted to 8 wells of a 96-well plate and 1 µl X-gal stock solution (40 mg/ml) was added to all wells except of the first. Different volumes of phage suspension were added according to the following table:

Number Phage Volume [µl] X-gal
1 Dickinson_Phage_001 10 no
2 No Phage - yes
3 Dickinson_Phage_001 1 yes
4 Dickinson_Phage_001 10 yes
5 Positive Ctrl 1 yes
6 Positive Ctrl 10 yes
The 96-well plate was placed in the plate reader and the phage detection protocol was started measuring the absorption at 600 nm every 10 min for 2 h.
There was no significant difference in absorption between the infected wells and the non-infected control as well as the well without X-gal. The results demonstrated that the conversion of X-gal is not a suitable method for phage quantification by spectroscopy as the product absorbs light at 650 nm, which interferes with the absorbance of cells. This explains the continuous increase of absorption during 2 h of cultivation. Further tests are planned using the beta-galactosidase substrate ONPG as the product absorbs at 420 nm.

Plaque Assay Test

To simplify the method, plaque assays were performed using only top agar. Top agar (0.7 %) supplemented with carbenicillin and X-gal stock solution was mixed with E. coli (ID:47) and phage suspension, according to the protocol: Blue Plaque Assay, and was poured into an empty well of a six well plate.
Plate layout:

Number Phage Dilution Comments
1 Positive Ctrl 10-11 -
2 Positive Ctrl 10-12 sample was frozen
3 - - Cells Only
4 - - Top Agar Only
The plate was stored at 37 °C overnight.
On the next day top agar with E. coli was slightly blue, but only the positive control 10-11
showed 12 clear plaques. This proves the assumption that frozen phage samples are not longer infectious.

Software

DeeProtein

  • We implemented 1d convolutions to better capture the 1-dimensional character of protein sequences.
  • We generated a new test case to harness the accuracy issue for DeeProtein (the CNN was not able to grasp differences in the classes for the toplevel EC dataset). The new test case looks like this:

EC-number Function n(samples)
1.1.1.- With NAD(+) or NADP(+) as acceptor 2558
2.1.1.- Methyltransferases 3739
2.7.1.- Phosphotransferases with an alcohol group as acceptor 2137
4.2.1.- Hydro-lyases 3242
6.1.1.- Ligases forming aminoacyl-tRNA and related compounds reduced to 2913 (og: 11652)

Split_Polymerase

Dickinson-PACE:

Several plaque PCRs were performed using plaque assays of Dickinson-PACE samples. The sequence of interest was amplified using PCR reaction and sended to GATC for sequencing analysis. Sequencing revealed that we got an evolved phage from the Dickinson lab. This phage obviously not showed any mutations or differences after the seven days of evolution in our PACE device. As results we can retain the facts that we could not determine the danger of washout. On top of that, we have learned how to handle the turbidostats, as well as the lagoons. Beyond that, we are now familiar with the workflow the apparature. We know how to propagate phages and moreover avoid turbidostat contamination with phages.

PREDCEL:

First PREDCEL runs and the associated plaque assays showed that we are not running into wash out for all of the three different amounts of transfered supernatant (10 µl, 100 µl and 1000 µl). On top of that, the phages are propagating well, as we can have similar phage titers for the different timepoints (20, 30 and 60 min). Our second experiment should evaluate, whether 2YT, LB or Esvelt-Medium should be used for further experiments. Since Kevin Esvelt reported in his PhD thesis, that medium influences the amount of F-pili an E.coli expresses. This should be directly correlated to the number of phages which can infect one of the E.colis, since M13 phages are using f-pili for their transport in and out of the cell. The results of these tests showed a high number of single plaques for 10-2 and 10-4 dilutions of the supernatant from LB medium. In contrast, 2YT and Esvelt medium showed uncountable numbers for these dilutions, with a slightly higher propagation using Esvelt medium. Exceptions have been made, since we had poblemes with the negative control in all cases. Since we want to observe the general feasibility of each medium to propagate phages, this should not be important in this case.

Reporter systems:

Four different reporter systems should be integrated in our assesory plasmids. mNeonGreen is currently available in all our APs, while the Nano Luciferase, provided by Promega, will be cloned in one of Jan´s APs. This reporter should supply a bright signal while having a moderate size. LuxAB was former used in the Bt toxine paper from Badran et al.. All our important strains contain the LuxCDE genes in their genome (S1059, S1030, S2060, S2208). We should consider it as a reporter system. At last, LacZ is considered as another possibility for a simple readout before infecting the laggons. A possible reaction could include the addition of phages to a culture which contains AP with LacZ. After an amount of time, for example ten minutes, phage infections should be observed. For the next week, cloning of the different reporter systems is planed.

MinION sequencing:

  • spoke to Jan-Philipp Mallm who uses the MinION in AG Rippe
  • he is sitting in the fourth floor (R425 or something)
  • told me that the device is, like we already know, leaky in terms of error probability (around every tenth base is called incorrect)
    • first he was pessimistic to detect small mutations
    • BUT one could calculate the number of random mutations by comparing the occured mutations over time
      • one should find the directed mutations by the number of occurences and the occurence over time
  • due to that we should not need barcording - is not useful
    • maybe barcoding for different sample pools from different timepoints - simultanous sequencing of several samples
  • 99.7 % accuracy is possible using a 50x coverage (Oxford nanopore)
  1. DNA Lobind tubes from (n irgendwas) are essential and will be needed
  2. first part of each oxford nanopore protocol is for the generation of cDNA
    • we should start with generating phage DNA (at least 1 µg) - in a first experiment Max purificated 3 µg of phage ssDNA
    • Jan told us to use random priming for generation of dsDNA for sequencing
  3. NEB Ligase Kit is needed as well as the TA-tailing Kit, which we have samples for (RNA library prep kit for Ilumina by NEB) in -20 °C
  4. we should simply follow the protocol which is straight forward in this concern

To do:

  • order all the needed stuff (with Roland´s permission)
    • DNA low binding tubes
    • NEBNext Ultra II End-repair/dA-tailing module
    • NEB Blunt/TA Ligase Master Mix
    • ThermoFisher Scientific Dynabeads
    • Klenow Fragment
    • Random hexamers
  • think about a possibility to use

Tetracyclin_Destructase

Destructase PACE Cloning I:

Weekly Report 23.07.17 - 30.07.17 PP

Destructase PACE consists, as all the other PACE experiments, of three plasmids, namely the accessory plasmid (AP), the selection phage (SP), and the mutagenesis plasmid (MP).
Design of the AP and SP is described in the following lines.

SP

The gene of interest in this PACE run is an anhydrotetracycline destructase named Tet(x) from Bacterioides fragilis.
We aim at destructing Tet(x) and thereby increasing anhydrotetracycline (aTC) concentration and gIII expression. The sequence was codon optimized for E. coli and a strong, self-created RBS was calculated (Salis Lab). The selection phage is re-assembled because stop-codon creating mutations within tet(x) should be prevented.
Therefore gII which is essential for phage replication by introducing a nick into the circular phage genome lin1972role was linked by an GGSGGS directly behind tet(x).

AP

For phage and anhydrotetracycline dependent geneIII (gIII) expression a fusion promoter (Carlson et. al) consisting of a phage shock promoter (psp) and anhydrotetracycline induced domain was used (Ppsp-tet).
Gene III is followed by mNeonGreen, a reporter fluorophore to measure gIII expression level. The backbone consists of a E. coli origin of replication (ColE1 ori) and an ampicilline resistance. A non-coding sequence has been added between backbone and Ppsp-tet because all fragments are assembled by Gibson assembly following the cloning standard.

Destructase PACE Cloning II:

1. AP

  • The psp-tet promoter sequence is taken from Carlson et. al (Negative Selection) and was ordered as a gblock. psp-tet promoter is flanked by gibson overlap 1 on the 5´end and gibson overlap 2 on the 3´end.
  • GeneIII is extracted from a Dickinson plasmid (Jin_177:_p15A_origin_pG). On the 5´end geneIII is concatenated to gibson overlap 2 and on the 3´end to gibson overlap 3.
  • mNeonGreen or any other reporter system is taken from the PAM-PACE accessory plasmid. mNeonGreen gblock includes gibson overlaps 3 and 4 at 5´and 3´end.
  • The backbone includes a ColE1 ori and AmpR and is extracted from pSB1A3. By using primers with extensions gibson overlap 4 is added to the 5´end and gibson overlap 5 to the 5´end.
  • The non-coding fragment between Gibson overlap 5 and 1 is extracted from the iGEM Distribution Kit and contains 256 bp from BBa_K808003, a small subunit B1 of the tripartite tricarboxylate transporter family.
    With primers PP_005_fwd and PP006_rev a fragment from coding sequence is extracted (bp 78-334). The extracted fragment consists of 265 bp and with the named primers gibson overlap 5 is added to the 5´end and gibson overlap 1 is added to the 3´end.
    (PCR Protocol PP_001)
Primers 10 µM Template 10ng/µl

Name Volume [µl]
Template BBa_K808003 0.5
Primer forward PP_007_fwd 2.5
Priner reverse PP_008_rev 2.5
Master Mix Phusion Flash 2x 25
ddH2O ad 50 µl 19.5
The PCR product was purified with die Quiagen PCR Purification Kit. To assemble the five fragments using the gibson overlaps a CPEC (circular polymerase extension cloning)
was performed.
  • low ramping temperature 0.5 °C/s
  • Q5

Name amount [ng]
backbone pSB1A3 75
insert 1 BBa_K808003 2x equimolar
insert 2 psp-tet promoter 2x equimolar
insert 3 geneIII equimolar
insert 4 mNeonGreen equimolar

step temperature [°C] time [s]
initial denaturation 98 30
denaturation 98 20
annealing 55 20
elongation 72 150
final extension 72 300
To test whether the CPEC worked out the CPEC product was transformed by chemical transformation (iGEM protocol) into TOP10 cells and plated out on ampicilline agar plates. The plates incubated at 37 °C overnight but after 20 hours of incubation no colonies grew up.
Because the CPEC did not work after several attempts the whole fragment between gibson overlap 2 and 5 was amplified from PAM-PACE accessory plasmid (AP_Cas9_pSB1A3_NGG_SD4) following this PCR protocol

Name Volume [µl]
Template AP_Cas9_pSB1A3_NGG_SD4 0.5
Primer forward JM_053 1.25
Priner reverse JM_073 1.25
Master Mix Q5 Master Mix 2x 12.5
ddH2O ad 25 µl 9.5

step temperature [°C] time [s]
initial denaturation 98 30
denaturation 98 10
annealing 66 20
elongation 72 150
final extension 72 120
The PCR did not work out and therefore, the same PCR was performed with different concentrations of Betain (0,5 M and 1 M) and DMSO (3 %, 5 % and 10 %)
This did not work as well and therefore, new primers are designed to anneal at 72 °C (PP_013 and PP_014). The PCR worked out and a gel extraction was performed and the PCR product AP_destructase_g2-g5_PP013-PP014 is used for further CPEC.

weekly_report_organosilicones

Organosilicones

weeks before

After we decided in our group to pursue the challenge to build Si-C bonds in vivo building on the intriguing paper >„Directed evolution of cytochrome c for carbon–silicon bond formation Bringing silicon to life“ by Francis Arnold (2016) and use PACE to lastly modify the catalytic enzyme, we began forming a strategy. Combining various literature methods, our plan more and more assumed a form. We wanted to use a riboswitch, detecting the Si-C bond containing small molecule Rma cytochrome c forms, to logically link the specificity and activity of our evolving protein to the Phage production.
To start this, we wanted to use the mutated version of cyt c to enantio- and chemoselectively produce our target molecule.
Next, the so generated molecule should enable the direct evolution of our RNA riboswitch with the method SELEX (Systematic Evolution of Ligands by EXponential Enrichment). Lastly, the respective implementation of the riboswitch and cyt c DNA in PACE should give rise to new mutants with an improved catalysing mechanisms for our desired molecule. The educt for the reaction requires a molecule with a diazo functional group plus one with a Si-H bond. The ethyl-diazo propionate forms, under nitrogen split-off, a carben intermediate with the iron ion of the heme cofactor of the cytochrome c. This in turn insertes itself specifically in the Si-H bond, giving rise to a product with a chiral carbon which is now covalently bound to the Si atom of the other molecule. After we consulted Reiner Tacke for a deeper insight of the relevance of Si-C bonds in industry, we decided that especially enantio- and chemselective biocatalyst are of high importance, narrowing the pool of possible substrates for our reaction. Me-EDA, are used in the paper, was commercially available, whereas 4-(dimethylsilyl)-anilin is not. The expert opinion of various professors of the chemical institute of our university implied that the synthesis is to dangerous for our team internal biochemist, who has a bit of organic lab experience. So we lastly found a group willing to cooperate and synthesising our molecule. Later on they also supported us regarding general chemical matters of this subproject. For the SELEX part, we got help of Dr. Murat Sunbül from the Jäschke Lab, who is experienced in aptamer development.

week 10-16.07.2107

The plan was to produce the enantioselective molecule we want to SELEX against later on. To start with, we needed the protein. As described in the Supplementary Materials of >"Directed evolution of cytochrome c for carbon–silicon bond formation Bringing silicon to life“, further on referred to as S.M. (p.3 (B)), we expressed and purified the Rma cytochrome c TDE (triple mutant. Thereby IPTG initiated the expression of the T7 Polymerase in the bacterial strain BL21 (DE3), in turn triggering the production of cyt c from the pet22(b) plasmid. Purification was facilitated though the apparent colour of the heme cofactor of the protein, which coloured the sample brown-red (Bild1) and therefore enabled an easy identification of the wanted fractions.

week 17-23.07.1207

The purity of the previously extracted protein was determined via a NuPAGE () and the concentration through the ferrous assay (S.M). The reactions of Me-EDA and 4-(dimethylsilyl)-anilin were performed with purified cyt c as well as whole cell solution. Both samples were treated as described in the >“Preparative-Scale Whole-Cell Biocatalytic Reaction“, whereby anaerobic conditions were mimicked through Nitrogen Gas washing of the reaction vial. Purification of the product through silica column chromatography delivered 25 fractions. Simultaneously, for the SELEX part, the RNA-library given to us by the Jäschke Lab was amplified. Therefore, a large scale PCR was conducted with primer ordered by IDT, in whole taking a volume of 10 mL (Bild). Next, the DNA was purified through Ethanol precipitation.

week 24-30-07.2017

The from silica column chromatography obtained fractions were analysed first with DC and building on that, the most prosperous fraction was sent to the cooperating Grebe group for GC-MS. The spectra was analysed with the help of a chemist from the group, Fabian Ebner, revealing that the reaction did not successfully take place as only educt was detected by the GC-MS. After further informing ourselves about the exact reaction mechanisms and the conditions as far as published in the technical literature, we derived at two suspicions for the failure of there reaction. One was that the incativation of the catalytic activity of the protein through carben transfer to a protein side chain as described in >“Identification of mechanism-based inactivation in P450-catalyzed
cyclopropanation facilitates engineering of improved enzymes“ (H.Renata, 2016). Here the author states, that >“inactivation of whole-cell catalyst was previously found to be slower than purified enzyme“.
The other approach suggested that our protein was in its oxidised state due to lack of sufficient oxygen exclusion and no addition of reduction reagent to the samples, implying that Carben Intermediat formation is only possible with the reduced Heme group. For this reason, a set of test reactions were performed with the commercially available Dimethylphenylsilan as a replacement of the 4-(dimethylsilyl)-anilin, always requiring custom synthesis. Hereby, we referred to the procedure of >“Lysate-small scale biocatalytic reactions“, mainly differing from the previous experiment by the addition of Dithionite solution as a reducing agent and GC Septum vials as reaction vials. Besides four samples with purified cyt c under various conditions, two samples still contained the whole cell solution.
Calendar Week 32

Modeling

Interactive modeling works, the information on hover does not work yet. Interactive modeling is expanded to plotting user data in the models.
Numeric modeling of PredCel is implemented but only generates oscillations so far.

Phage_Propagation_and_Quantification

Detection of Destructase and Cas9 Phages

As we assumed, phage production is very inefficient Blue Plaque Assays were performed using more
phage suspension. 100 µl E. coli suspension (stock ID:47) and 200 µl phage suspension was mixed with 1500 µl top agar supplemented with carbenicillin and X-gal stock solution and plated on LB/Amp-plates.

Plate Phage
1 SP_Destructase (GG1)
2 SP_Destructase (GG2)
3 SP_dCas9 (GG1)
4 Positive Ctrl
5 Negative Ctrl
Positive control exhibited plaques and negative control was clear, but only one plaque was visible for SP_Destructase (GG2). Blue Plaque Assay was repeated using a 15 cm dish. 500 µl E. coli suspension and 1 ml phage suspension was mixed with 5 ml top agar supplemented with carbenicillin and X-gal stock solution and plated on 15 cm LB/Amp-plates.

Plate Phage
1 SP_Destructase (GG1)
2 SP_Destructase (GG2)
3 SP_dCas9 (GG1)
4 SP_dCas9 (GG2)
5 Positive Ctrl
6 Negative Ctrl
On the next day, positive control exhibited plaques and negative control had no plaques. SP_dCas9 (GG2) showed around 50 plaques, but no plaques were visible on the other plates.

Software

Generation of datasets for GO-term classification

  • Filtered the Swissprot file and rewrote it to a .csv containing: name, seq, EC_nr, GO-terms
  • Then GO-term specific .csv files were generated. The GO-DAG (directed acyclic graph) was completed for each GO annotation in each protein, and the respective protein was then added to all GO-term specific files, that it has GO-annotations for.
  • Extremely generic GO-labels like: GO:0003674 molecular function, or GO:0005488 binding were excluded.

Testing of the padded input-pipeline

  • A lot more noise through the different padding sizes
  • The system is a lot more complicated as we're feeding directly from the queues and therefore have two different graphs, that share parameters

    Name Volume [µl]
    Template SP_CBXTAL supernatant 2
    Primer forward PP_001/PP_003 1.0
    Priner reverse PP_002/PP_004 1.0
    Master Mix Q5 Master Mix 2x 10
    ddH2O ad 25 µl 6

    step temperature [°C] time [s]
    initial denaturation 98 300
    10x
    denaturation 98 20
    annealing 64 20
    elongation 72 120
    25x
    denaturation 98 300
    annealing 72 20
    elongation 72 120
    final extension 72 300
    Both fragments were successfully amplified as to see gel electrophoresis picture XXXX. To assemble these two fragments and the Tet(x) gblock all primers included extensions with BsaI restriction sites and Golden Gate overhangs.
    The Golden Gate was performed following the Golden Gate protocol:

    Name amount [ng]
    backbone SP_2 75
    insert 1 SP_1 equimolar
    insert 2 Tet(x) promoter equimolar

    temperature [°C] time [min]
    15x
    37 2
    16 1
    1x
    37 30
    55 5
    The Golden Gate product is purified using the ZymoResearch Kit and then 1 µl is
    electroporated into 25 µl of S1059 E. coli. After 4 hours incubation in 3 ml SOC medium the cell culture is used for an plaque assay. A plaque PCR showed that the Tet(X)-insert is not included in the Golden Gate product and
    therefore new Golden Gate overhangs are designed.
Calendar Week 33

Cytochrome_Engineering

1. Cloning of SP Cytochrome PACE:

All required parts for the selection phage were amplified in sufficient concentrations by PCR to be assembled by Golden Gate cloning.

1 2 3
Q5 MM (2x) 25 µl 25 µl 12 µl
CG_011 (10 µM) 2.5 µl 2.5 µl 2.5 µl
SP_R (10 µM) 2.5 µl 2.5 µl 2.5 µl
SP_ATMTAL 1 µl 1 µl 1 µl
DMSO 2.5 µl - 2.5 µl
Betain - 5 µl 5 µl
The cycler conditions were as follows:

Phase Temperature [°C] Time [min] Cycles
Initial Denaturation 98 0:30 1
Denaturation 98 0:10 35
Annealing 72 0:25 35
Extension 72 3:00 35
Final Extension 72 2:00 1
Hold 10 1
The PCR fragments were analyzed by gel electrophoresis. There were no bands on the gel. Maybe the primer concentration, which was 1 µl instead of 10 µl was too small. PCR was repeated with two different polymerases, Q5 and Phusion Flash and two different conditions in a 25 µl mix.

1 2
Polymerase MM (2x) 25 µl 25 µl
CG_011 (10 µM) 2.5 µl 2.5 µl
SP_R (10 µM) 2.5 µl 2.5 µl
SP_ATMTAL 1 µl 1 µl
Betain - 5 µl
The cycler conditions for Q5 were the same as above, the conditions for Phusion Flash were as follows:

Phase Temperature [°C] Time [min] Cycles
Initial Denaturation 98 3:00 1
Denaturation 98 0:01 35
Annealing 72 0:05 35
Extension 72 1:30 35
Final Extension 72 1:00 1
Hold 10 1
The reaction mix with Q5 did not exhibit any bands on the gel, but with Phusion Flash there were bands of the expected length (6.198 kb) visible on the gel. The reaction mix without betain exhibited a stronger band. The PCR product was purified with the QIAquick Gel Extraction Kit and eluted in 30 µl ddH20 with 9.6 ng/µl DNA.

APs Opto

PCR amplification of AP fragments:

For the generation of the two Opto APs, AP_light and AP_dark three fragments are necessary:

AP_light AP_dark
1 gIII_luxAB gIII_luxAB PCR product
2 AP_light_BB AP_dark BB PCR product
3 pBLind pBLrep oligo
gIII_luxAB: In order to amplify gIII_luxAB, primers gIII_luxAB_F and gIII_luxAB_R were used together with AP_ATML (plasmid ID:5) as template in a 50 µl Q5 PCR reaction. The cycler conditions were set as follows to amplify the 3.784 kb gIII_luxAB fragment:

Phase Temperature [°C] Time [min] Cycles
Initial Denaturation 98 0:30 1
Denaturation 98 0:10 35
Annealing 72 0:25 35
Extension 72 2:00 35
Final Extension 72 2:00 1
Hold 10 1
Analysis by gel electrophoresis demonstrated a band of the expected length. The PCR product was purified with the QIAquick Gel Extraction Kit and eluated in 30 µl ddH20 with 52.5 ng/µl DNA. AP_light_BB and AP_dark_BB: To amplify AP_light_BB and AP_dark_BB a 50 µl Q5 reaction mix was prepared using AP-ATML as template and the respective primers for the APs:

AP Primer forward Primer reverse
AP_light_BB BB_F CG_005
AP_dark_BB BB_F CG_008
The cycler conditions were set as follows to amplify the 2,74 kb AP_light_BB and AP_dark_BB fragments:

Phase Temperature [°C] Time [min] Cycles
Initial Denaturation 98 0:30 1
Denaturation 98 0:10 35
Annealing 72 0:25 35
Extension 72 2:00 35
Final Extension 72 2:00 1
Hold 10 1
After the cycler run, 5 µl of the reaction mix was load on an agarose gel and detected by UV light. AP_light_BB and AP_dark_BB exhibited a strong and clear band at 2,7 kb, therefore PCR product was purified with the QIAquick PCR Purification Kit.

AP concentration [ng/µl]
AP_light_BB 166.0
AP_dark_BB 146.8

Phage_Propagation_and_Quantification

Phage propagation of the unevolved Dickinson phage

Phage supernatant of the unevolved Dickinson phage: target_133_N-term_T7-C was received from Dickinson group. To propagate the phages, 4 ml E. coli culture (Stock ID: 47) was cultivated to an OD600 of 0.6 (in LB media + 25 mM glucose + amp) and infected with 4 µl of the phage supernatant. Culture was shaked at 37 °C overnight. On the next morning, culture was centrifuged at 6000 xg for 5 min and supernatant, which contains the phages, was stored at 4 °C. A Blue Plaque Assay was performed to determine the phage titer of the supernatant. 143 plaques were counted at the 10-10 dilution, which leads to a phage titer of 1.43*1015 PFU/ml. A plaque of this plate was picked to infect a 4 ml E. coli culture (Stock ID: 47) with an OD600 of 0.4. This culture was cultivated for two hours shaking at 37 °C and was subsequently transferred to 100 ml fresh 2xYT medium. After 1 hour carbenicillin (1000x) was add. On the next day, culture was centrifuged at 3640 xg for 20 min. Supernatant was stored at 4 °C. A Blue Plaque Assay was performed to determine the phage titer of the supernatant. The monoclonal Dickinson phage target_133_N-term_T7-C exhibited a phage titer of 1.85*109 PFU/ml.

Software

KW33

Generation of datasets for 427-classes and 299 classes

In order to further simplify the problem of having a high amount of weakly represented classes we generated two additional datasets:
427 EC classes with a minimum sample size of 100, and 299 EC classes with a minimum samplesize of 200.

Implementation of first generative Model

To reach the lowest hanging fruit in generative modeling, we switched DeeProtein to the generative mode. Generative modes take an input (or noise input) and return a modified version of the inout maximizing the class probability for the target class. In the context of protein sequences, generative modeling could for instance generate new funcitonalities by "shifting" a sequence in the multidimensional featurespace aling the direction of the target class.
The easiest way to achieve functionality is by just "training" the input of the network. By freezing all the weights and allowing only the input to be modyfied iteratively, a optimized sequence is obtained. Results have been disentchanting.

Split_Polymerase

Dickinson-PACE:

For implementation of the LuxAB reporter system in the Dickinson strains, golden gate cloning was implemented. The resulting strains had LuxAB integrated directly behind gene III. This should allow us to investigate the phage infection directly correlated to the LuxAB expression. These slightly different versions were transformed into S2060 cells for use in the PACE apparature. Another change of the experimental design was implemented using MP6 instead of MP1. MP6 was transformed into the S2060 used for the Dickinson transformations. This change was made in order to generate more mutations than before. On top of that, MP6 allows to implement less genetic drift, as well as no stepping stones for the evolution. The general design of the PACE experiment itself was held constant, including strain swapping after three and four and a half days. Currently we have no wash out and no phage contamination in the turbidostat according to the phage detection PCR.

Reporter system:

The first testing of the S2060 + MP6 + alternative Dickinson-APs showed only "special" results. Generally, we can detect a luminescence signal. However, the signal of the phage supernatant is higher than the signal of the infected cells. Beyond that, the turbidostat displayed an higher luminescence signal in both lagoons. Based on previous results from Jan, we suppose that a phage titer of 10hoch10 is required for a sufficient "real-time" readout. Due to that, "real-time" PACE will not be possible. Based on other results from different LuxAB testing, where LuxAB was tested under the psp-AR2 promotor, one can suppose that leaky expression of the T7 promotor could be the reason. In contrast to the T7 promotor, the psp-AR2 promotor did not show any background signal. This underlines the previous hypothesis.

Toolbox impressions:

ThermoFisher:

ThermoFisher supplies a table containing different parts of which you can choose.

part number/name public name (gene name) part family/class of the part part category
Abbrevation (BBa_K747000) LacZ gene CDS reporter system

BlueHeron:

BlueHeron offers a broad amount of different vectors for several purposes (protein purification detection - epitope tags, protein visualization - fluorescence tags, resistances, expression vectors). Each of the higher categories is divided into the purposes. Every subclass is then characterized by the following table:

Entry vector SKU FP Feature N-Tag C-Tag E.coli Selection Cell selection Expression
name of the plasmid (pLenti-N-tRFP) abbrevation (BBa_K747000) RFP tRFP Myc-DDK resistance resistance Bacteria

IDT:

IDT only offers six different vectors, where your gene of interest will be inserted. Your sequence can be pasted and will be inserted in the vector after complexity scan and codon optimization.

vector name Size selection marker application sequence
name of the plasmid (pUCIDT (Amp)) in bp ampicillin cloning, in vitro transcription, etc. pUCIDT.txt/genbank etc.

Agilent:

SureVector Tool is the best example for our potential toolbox:
<code>* several fragments can be choosen
1. Origin (Kanamycine)
2. Promotor/Prom-tag (p15A)
3. Tag (yARS)
4. Bacterial selection (URA3)
5. XP1 Slot (PgalHis)
6. XP2 Slot (Gene of interest)
7. Upload sequence of interest (as file or by copy paste)
</code>

Several steps to create the vector map:

  1. Select a bacterial origin of replication
  2. Select a promotor or a promotor-tag fusion
  3. Select a tag (skip this step if you selected a promoter.tag fusion)
  4. Select a bacterial selection marker
  5. Select additional options for functionality in yeast or mammalian cells.
  6. Supply your 5´-3´ gene-of-interest in the sense strand direction by uploading a fastA file or pasting the sequence in the box
  7. Click create, review your configuration and order the corresponding SureVector kits.
Making your gene-of-interest Sure Vector compatible:
If amplifying from an existing template:
Forward primer and reverse primer:
  • two primers are provided for insertion into the Vector

To Do:

Protocol every different part we have:

1. Aktivator + Promotor

Each user can insert the required input sequence but can use several ready-to-use promotor systems like the psp-promotor.

2. RBS + geneIII

Different RBS sites in combination with gene III (different stringency)

3. RBS + reporter system

LuxAB, mNeonGreen, Lacz, as well as one or two other reporter system will be provided (estimation of leaky expression of the used circuit - NOT RT-PACE)

4. RBS + additional CDS with terminator and RBS

Another input sequence for the user depending on the protein of interest.

5. Resistance + origin of replication

Three different ORIs in combination with all our available resistances except from Chloramphenicol (always on the MPs)
  • for every PACE circuit document the different possibilities
  • Activator + Promotor:

  • PAM-PACE gRNA (binding sequence) with Promotor (muss noch kloniert werden)
  • pBlind(v2) (promotor Blue-light induced) (muss noch kloniert werden)
  • pBLrep(v2) (promotor blue-light repressed) (muss noch kloniert werden)
  • Upstream activating sequence II (UASII) + UASI + integrating host factor (IHF) + psp-tet (alles gesamt psp-tet) (muss noch kloniert werden)
  • T7 Promotor (muss noch kloniert werden)
Amplification with primers containing gibson overlaps 1 and 2 (g1 and g2) Fwd_primer_g1: TGAAATTCTGGGCCCTGCATTACGAGAACTGA
Rev_primer_g2: TACTGCACCTGTAGGTCGACGATGATCTGATA

2. RBS + geneIII:

  • SD8 (always with geneIII)
  • sd8
  • sd6
  • SD4
  • sd2
  • sd5 (muss noch kloniert werden)
  • Theophylline Riboswitch 12.1 (muss noch kloniert werden)
Amplification with primers containing gibson overlaps 2 and 3 (g2 and g3) Fwd_primer_g2: TATCAGATCATCGTCGACCTACAGGTGCAGTA
Rev_primer_g3: TTAGTACTACTCGAGTTGATAGGCACCGACCA

3. RBS + reporter system:

  • mNeonGreen-RBS + mNeonGreen + BBa B1002 and BBa B1006 Terminator
  • LuxAB + BBA B1002 and BBa B1006 Terminator
  • NanoLuciferase
  • Beta-Galactosidase
Amplification with primers containing gibson overlaps 3 and 4 (g3 and g4) Fwd_primer_g3: TGGTCGGTGCCTATCAACTCGAGTAGTACTAA
Rev_primer_g4: TGGTACTTCAAATGCGGCTTGGCTCCAGACAA

4. RBS + additional CDS with terminator and RBS:

  • RBS + sgRNA_expressioncassette_Cas9_RFP_TrrnB_terminator
  • BBa_J23105 + HDJ_Chaperone + BBa_B0010
  • Lac UV5 promotor SD8 RBS cI - Ser/Gly Linker-FKBP + SD8 RBS DHFR - Ser/Gly Linker - RPA
  • BBa_K808003 - (non sense)
  • Kat Promotor - linker - C-terminal T7-RNAP - Terminator
Amplification with primers containing gibson overlaps 5 and 1 (g5 and g1) Fwd_primer_g5: TCAGTTCTCGTAATGCAGGGCCCAGAATTTCA
Rev_primer_g1: TTATCCACTGGCGAGCTCTGTAACGAAACGTA The additional CDS dependent on the actual circuit which will be used. Due to that, the user has to insert a second sequence, where we can put an RBS in front of and a terminator behind.

5. Resistance + origin of replication:

  • pSC101 + Ampicillin
  • pBR322/ColE1 (=pSB1A3) + Ampicillin
  • p15A + Ampicillin
  • pSB1A3 + Tetracyclin (muss noch kloniert werden)
  • pSC101 + Tetracyclin (muss noch kloniert werden)
  • p15A + Tetracyclin (muss noch kloniert werden)
  • p15A + streptomycin (muss noch kloniert werden)
  • pBR322 + streptomycin (muss noch kloniert werden)
  • pSC101 + streptomycin (muss noch kloniert werden)
Amplification with primers containing gibson overlaps 4 and 5 Fwd_primer_g4: TTGTCTGGAGCCAAGCCGCATTTGAAGTACCA
Rev_primer_g5: TGAAATTCTGGGCCCTGCATTACGAGAACTGA
  • check for typeIIs restriction sites

Recommendation:

Use Ampcillin for the AP, Kanamycine for the SP and Chloramphenicol for MP
  • Chloramphenicol cannot be used because of the mutagenesis plasmids which contain a Cmp Resistance
  • tetracyclin cannot be used in combination with the psp-tet promotor

Selection Phage:

  • linearized selection phage without inserted gene of interest - fixed overhangs for golden gate cloning
  • Amplification of your gene of interest with primers A and B for cloning into the selection phage (either paulines primer - destruction or jans primer - for evolution of an gene-of-interest)
  • check for typeIIs restriction sites
    • provide only resistances without typeIIs
      • will be tested (muss noch gemacht werden)
      • cloning of SPs with different resistances (muss noch gemacht werden)
  • Alternative for selection phage if you want to analyze the essential amino acids:
    • use paulines primer with backbone
      • Primers contain BsaI restriction sites
  • The sequences in brackets are the overhangs for the new primers, the rest of the sequences are specific primers for paulines and Jans PACE (can be ignored for the toolbox)

Paulines Primer overhangs (for cloning of SP):

  • gene of interest:
    Fwd: [TTGGTCTCATATC]
Rev: [TTGGTCTCACAAT]TGAACCTCCAGAGC
  • backbone part_1:
    Fwd [GCTTGGTCTCAATTG]ACATGCTAGTTTTACGATTACC
Rev: [GCTCGGTCTCATAGA]TTAGAAAAACTCATCGAGCATCAAATG
  • backbone part_2:
    Fwd: [TTGCGGTCTCATCTA]GAAGGAGATTTTCAACATGCT
Rev: [TTGCGGTCTCAGATA]ACCCCGGTTGATAATAAG

Jans Primer overhangs (for cloning of SP):

  • gene of interest:
    Fwd:[CACCACAGGTCTCGAAAA]ATGGATAAGAAATACTCAATAGGCTTAGC
Rev: [CACCACAGGTCTCGGTGC]GTCACCTCCTAGCTGACTC
  • backbone part_1:
    Fwd: [CACAGGTCTCATCAG]GTCAGAAGGGTTCTATCTC
Rev: [CACCACAGGTCTCGTTTT]TTTTCCTCCTAAGCTGTTAGAAAAACTC
  • backbone part_2:
    Fwd: [CACCACAGGTCTCGGCAC]GCGTAACTGTTCAG
Rev: [CACAGGTCTCTCTGA]AAGCGTAAGAATACGTGGC
  • backbone should be supplied with kanamycin, ampicillin, streptomycin and tetracyclin
  • the antibotics cannot be the same in AP and SP

Examples of tested genes of interest:

  • Tet(X)
  • rpoZ-Cas9
  • CYP1A2
  • EL 222
  • Split-Cas9
  • T7-RNAP 1-179 (N-terminal) - Dickinson
  • RNAP Cytochrome C
For the integration of the toolbox in our project, in silico cloning using Geneious was performed. Primers were ordered to clone the different parts in our cloning standard as well as building up parts for the iGEM registry.

Tetracyclin_Destructase

AP

CPEC is performed with PCR product AP_destructase_g2-g5_PP013-PP014 as backbone and psp-tet promoter and and BBa_K808003 as inserts.
The CPEC product (AP_destructase) is electroporated S2060 electrocompetent cells an plated on an ampicilline plate. Incubation over night was successful and a colony PCR is performed.

step temperature [°C] time [s]
initial denaturation 94 30
denaturation 94 25
annealing 51 40
elongation 68 50
final extension 68 5

Name Volume [µl]
Template AP_destructase pipette tip
Primer forward JM_072 0.5
Priner reverse JM_065 0.5
Master Mix OneTaq 2x Master Mix 12.5
ddH2O ad 25 µl 11.5
PCR proofed that the insertes are present in AP_destructase and a glycerol stock is created.
Calendar Week 34

Cytochrome_Engineering

1. Cloning of AP Cytochrome PACE:

As the assembly of the AP worked, the next aim was to produce AP variants with different ORIs. To achieve this, PCRs of the backbone (all plasmid parts except the ORI and AmpR) were performed. Unfortunately problems occured, which is the reason why series of following PCRs have been performed under different conditions (different polymerases, different annealing temperatures, adapted cycler runs). After all those actions had no promising effect we assumed, that there must have been a mistake in the previous procedure.
In order to repeat the CPEC, one part (LuxAB) was still needed. So in the following PCRs were performed to amplify the missing LuxAB. This again beared problems, so different conditions were tested again (different polymerases, different primer concentrations, DMSO, Betain, different cycler runs, touch down PCR, slow ramping after denaturation)
After none of this worked, we decided to chose mNeonGreen as a reporter instead of LuxAB, as this one has already been amplified.
Now all necessary parts for the CPEC were abundant. The CPEC product was then transformed into chemically competent cells. On the next day colony PCRs have been conducted with 6 clones and simultaneously the corresponding growths were purified.

2. Cloning of SP Cytochrome PACE:

In parallel to the Plaque Assys, plaque PCRs have been performed to check for the correct SP. These PCRs unfortunately did not show the correct outcome. Also new Primers have been designed and ordered, that amplify a smaller fragment of the SP making the test PCRs more efficient.
When the new primers arrived test PCRs of the phage supernatant were performed, as an inefficient PCR could not be clearly excluded as source of error at this point. As this PCR also did not show the correct bands we decided to repeat the SP assembly with new backbone template. Those templates needed to be generated first. So PCRs from SP CBX TAL were conducted. These PCRs did not lead to the desired result and were repeated with different conditions (different polymerase, DMSO, Betain, cycler run). The changed conditions helped optimize the PCRs so that the products could be isolated from the agarose gel.
Afterwards the Golden Gate could be repeated and the isolated GG product was tranformed into 2060 electro competent cells. When recovering after the transformation, after 4 hours of incubation at 37�C a centrifugation step followed to gather the supernatant containing SP. Using the supernatant, a Plaque Assay was performed.

Modeling

A first numeric model of PredCel works without oscillations. Graphs look reasonable so far. Modelling was performed on three levels: On the lowest level one Step of Predcel monitoring Phage concentration as well as uninfected, infected and phage producing E. coli concentration Graph of level 1 Predcel model. One level above all the concentrations were tracked over 100 iterations of Predcel Graph of level 2 Predcel model. And on the third level different sets of values for starting fitness, starting phage concentration and starting E. coli concentration were tested. In this case we only monitored how long the phage titer at the end of each iteration of Predcel stayed above 1 pfu/mL and below 1e8 pfu/mL.Graph of level 3 Predcel model
At least the two higher levels probably only work in python, but maybe an interactive version of what happens during one iteration is possible. A final more comfortable version of the script was started.

Optogenetics

SP Opto

PCR amplification of SP fragments:

The blue-light mediated transcription activator EL222 was amplified from a gblock (EL222_expression_cassette) using CG_012_fwd and CG_013_rev as primers in a 25 µl Phusion Flash reaction mix.
The cycler conditions for the 678 bp fragment were as follows:

Phase Temperature [°C] Time [min] Cycles
Initial Denaturation 98 0:10 1
Denaturation 98 0:01 35
Annealing 72 0:05 35
Extension 72 0:11 35
Final Extension 72 1:00 1
Hold 10 1

Oligo annealing and phosphorylation of the AP promoters:

pBLind: CG_010_fwd, CG_009_rev
pBLrep: CG_007_fwd, CG_006_rev 5 µl of each primer (CG_010_fwd, CG_009_rev for pBLind and CG_007_fwd, CG_006_rev for pBLrep) (100 µM) was added to 1.1 µl T4 ligase buffer in a PCR tube. Reaction mix was heated up to 98 °C and slowly cooled down (0.1 °C/s). The annealed oligos were phosphorylated using the T4 PNK according to the neb protocol "Non-radioactive Phosphorylation with T4 PNK or T4 PNK (3´ phosphatase minus)".

Golden Gate assembly of SP and AP:

GG of Opto SP (20 µl Reaction):

Reagents Volume [µl]
SP_BB (Purification ID: 515) 15
EL_222 (Purification ID: 520) 2 (of a 1:10 dilution)
BsaI 0.5
T4 Ligase 0.5
t4 Ligase Buffer 2
GG of Opto APs (20 µl Reaction): AP_light

reagents Volume [µl]
gIII_luxAB (Purification ID: 514) 2
AP_light_BB (Purification ID: 516) 0.5
pBLind 1 (of a 1:100 dilution)
BsaI 0.5
T4 Ligase 0.5
t4 Ligase Buffer 2
AP_dark

reagents Volume [µl]
gIII_luxAB (Purification ID: 514) 2
AP_dark_BB (Purification ID: 517) 0.5
pBLrep 1 (of a 1:100 dilution)
BsaI 0.5
T4 Ligase 0.5
t4 Ligase Buffer 2
Cycler conditions were as follows:

Temperature [°C] Time [min] Cycles
37 3:00 15
16 4:00 15
37 30:00 1
65 0:05 1

Trafos of APs and SP:

APs were transformed into Top 10 cells (chemically competent) and SP was transformed into S1059 (electrocompetent). AP trafos were plated on LB-plates containing Amp. SP trafos were cultivated in SOC for 4 h and a Blue Plaque Assay was subsequently performed with the phage supernatant. The phage supernatant was also used to perform a PCR of the phage insert. Therefore, a 25 µl Q5 PCR reaction was prepared. To amplify the 1.9 kb fragment containing EL222, the primers JM_068_rev and JM_069_fwd were used together with 2 µl of the phage supernatant as template. The cycler conditions were as follows:

Phase Temperature [°C] Time [min] Cycles
Initial Denaturation 98 3:00 1
Denaturation 98 0:10 35
Annealing 64 0:25 35
Extension 72 0:50 35
Final Extension 72 2:00 1
Hold 10 1
Analysis by Gel electrophoresis revealed two bands at 1 kb and 1.9 kb.

Software

Word2Vec Embeddings on Proteinsequences

We rewrote a word2vec implementation from tensorflows tutorials that implements Efficient Estimation of Word Representations in Vector Space, ICLR 2013 (Mikolov, et. al.). The model is a skipgram model with negative sample that uses custom ops written in C. The code was adapted to our needs, mainly by changing datatypes in the C kernels and writing a different evaluation function based on predicting the nearest words to the most frequent words instead of using analogies. Two new datasets were generated based on both swissprot and uniprot. Training of 4mer embeddings in 50, 100 and 200 dimensions were started but have not been calculated yet.
Visualisation of the first checkpoints is possible via tensorboard Visualisation of an example embedding via tensorboard.

IMPLEMENTATION OF SQUEEZENET Architecture

With implamentation of a new architecture based on Sequeeze-net (Forrest N. Iandola, 2017), relying on 1x1 convolutions we were able to grasp the 299 as well as the 637 classes dataset. The new model architecture looks the following:
  • InputLayer model_valid/input_layer_valid: (64, 20, 1000, 1)
  • PadLayer model_valid/block1/pad_layer_valid: paddings:[[0, 0], [0, 0], [3, 3], [0, 0]] mode:CONSTANT
  • Conv2dLayer model_valid/block1/cnn_layer_valid: shape:[20, 7, 1, 128] strides:[1, 5, 1, 1] pad:VALID act:prelu
  • Conv1dLayer model_valid/block2/cnn_layer_valid: shape:[6, 128, 128] stride:1 pad:SAME act:prelu
  • Conv1dLayer model_valid/1x1_I/1x1_valid: shape:[1, 128, 64] stride:1 pad:SAME act:prelu
  • BatchNormLayer model_valid/1x1_I/batchnorm_layer_valid: decay:0.900000 epsilon:0.000010 act:identity is_train:False
  • Conv1dLayer model_valid/block3/cnn_layer_valid: shape:[5, 64, 256] stride:1 pad:SAME act:prelu
  • PoolLayer model_valid/block3/pool_layer_valid: ksize:[2] strides:[2] padding:VALID pool:pool
  • BatchNormLayer model_valid/block3/batchnorm_layer_valid: decay:0.900000 epsilon:0.000010 act:identity is_train:False
  • Conv1dLayer model_valid/block4/cnn_layer_valid: shape:[5, 256, 256] stride:1 pad:SAME act:prelu
  • PoolLayer model_valid/block4/pool_layer_valid: ksize:[2] strides:[2] padding:VALID pool:pool
  • BatchNormLayer model_valid/block4/batchnorm_layer_valid: decay:0.900000 epsilon:0.000010 act:identity is_train:False
  • Conv1dLayer model_valid/1x1_II/1x1_valid: shape:[1, 256, 128] stride:1 pad:SAME act:prelu
  • BatchNormLayer model_valid/1x1_II/batchnorm_layer_valid: decay:0.900000 epsilon:0.000010 act:identity is_train:False
  • Conv1dLayer model_valid/block5/cnn_layer_valid: shape:[5, 128, 256] stride:1 pad:SAME act:prelu
  • PoolLayer model_valid/block5/pool_layer_valid: ksize:[2] strides:[2] padding:VALID pool:pool
  • BatchNormLayer model_valid/block5/batchnorm_layer_valid: decay:0.900000 epsilon:0.000010 act:identity is_train:False
  • Conv1dLayer model_valid/block6/cnn_layer_valid: shape:[5, 256, 512] stride:1 pad:SAME act:prelu
  • PoolLayer model_valid/block6/pool_layer_valid: ksize:[2] strides:[2] padding:VALID pool:pool
  • BatchNormLayer model_valid/block6/batchnorm_layer_valid: decay:0.900000 epsilon:0.000010 act:identity is_train:False
  • Conv1dLayer model_valid/1x1_III/1x1_valid: shape:[1, 512, 256] stride:1 pad:SAME act:prelu
  • BatchNormLayer model_valid/1x1_III/batchnorm_layer_valid: decay:0.900000 epsilon:0.000010 act:identity is_train:False
  • Conv1dLayer model_valid/block7/cnn_layer_valid: shape:[5, 256, 516] stride:1 pad:SAME act:prelu
  • PoolLayer model_valid/block7/pool_layer_valid: ksize:[2] strides:[2] padding:VALID pool:pool
  • BatchNormLayer model_valid/block7/batchnorm_layer_valid: decay:0.900000 epsilon:0.000010 act:identity is_train:False
  • Conv1dLayer model_valid/block8/cnn_layer_valid: shape:[5, 516, 1024] stride:1 pad:SAME act:prelu
  • PoolLayer model_valid/block8/pool_layer_valid: ksize:[2] strides:[2] padding:VALID pool:pool
  • BatchNormLayer model_valid/block8/batchnorm_layer_valid: decay:0.900000 epsilon:0.000010 act:identity is_train:False
  • Conv1dLayer model_valid/1x1_IV/cnn_layer_valid: shape:[1, 1024, 512] stride:1 pad:SAME act:prelu
  • BatchNormLayer model_valid/1x1_IV/batchnorm_layer_valid: decay:0.900000 epsilon:0.000010 act:identity is_train:False
  • Conv1dLayer model_valid/block9/cnn_layer_valid: shape:[5, 512, 1024] stride:1 pad:SAME act:prelu
  • PoolLayer model_valid/block9/pool_layer_valid: ksize:[2] strides:[2] padding:VALID pool:pool
  • BatchNormLayer model_valid/block9/batchnorm_layer_valid: decay:0.900000 epsilon:0.000010 act:identity is_train:False
  • Conv1dLayer model_valid/outlayer/cnn_layer_valid: shape:[1, 1024, 637] stride:1 pad:SAME act:prelu
  • BatchNormLayer model_valid/outlayer/batchnorm_layer_valid: decay:0.900000 epsilon:0.000010 act:identity is_train:False
  • MeanPool1d global_avg_pool: filter_size:[7] strides:1 padding:valid
The architecture is fully convolutional, ending in an average pooling layer as outlayer, with the channels dimension corresponding to the number of classes. All inputs were 1-hot encoded and zero padded to a boxsize of 1000 positions.

Model lr classes Comment restored maxstep boxsize ACC
0.01 299 NO 220000 1000 0.8 (valid)
0.01 637 NO 180000 1000 0.55 (valid)
0.01 637 YES 35000 1000 0.75 (valid)

References:

  1. Iandola, F. N., Han, S., Moskewicz, M. W., Ashraf, K., Dally, W. J., & Keutzer, K. (2016). SqueezeNet: AlexNet-level accuracy with 50x fewer parameters and< 0.5 MB model size. arXiv preprint arXiv:1602.07360.
  2. Mikolov, T., Sutskever, I., Chen, K., Corrado, G. S., & Dean, J. (2013). Distributed representations of words and phrases and their compositionality. In Advances in neural information processing systems (pp. 3111-3119).

Tetracyclin_Destructase

AP

To test under which conditions the selection phage propagates best, 12 variations of the destrucatase accessory plasmid are assembled by gibson assembly.
The variations include different RBS and origins of replication. Because time was running out we decided to end this project at that point and concentrated on other PACE and PREDCEL experiments.
Calendar Week 35

Cytochrome_Engineering

1. Cloning of AP Cytochrome PACE:

After the purification showed extremely low DNA concentrations, the transformation of the AP was repeated using the previously purified AP product. Afterwards the performed colony PCRs did not show the expected results. So a test digest was conducted which also did not show the desired results. The next day the colony PCR was performed again using a different polymerase. After this did not lead to the expected result, the CPEC was conducted one more time. Loading the CPEC product on an agarose gel sadly showed that this repeated CPEC also did not work.
Starting right from the beginning, new PCRs have been performed in order to generate all needed parts to repeat the CPEC once again. After those PCRs finally worked, the products were extracted from the agarose gel to avoid the impact of primer dimers in the following CPEC.
The CPEC was conducted under a different strategy this time as it did not work several times before. Instead of using all fragments/parts in one CPEC reaction, the procedure was split in 3 CPECs as a whole. In the first 2 CPECs half of the parts were assembled together respectively. Then the AP parts were assembled toghether in the third CPEC to form the full length AP. This product was then again transformed into chemically competent cells.
After the transformation did not work, the 3 separate CPECs have been repeated once again. This time the third CPEC from the previous attempt and from the new attempt have been transformed into two different strains of chemically competent cell respectively to make sure that the cells used for the latest transformation are ok.
This time the transformations with the second strain of chemically competent cells worked in both cases, the new and the previous CPEC. Test PCRs followed of several clones of both transformations. Some of those clones showed bands at the correct height. Subsequently, the corresponding growhts containing the AP plasmids were preparated in order to isolate the APs and one of them was sent in to GATC to be sequenced.

2. Cloning of SP Cytochrome PACE:

On the next day, one Plaque was picked to be tested by PCR using the new Primers. The result was not what we were hoping for. So 10 more Plaques were picked and tested by PCR. Three of these Plaques showed the correct bands and a PCR purification was performed with those three PCRs. Two of the purifications had suitable concentrations and were sent to GATC for sequencing. The Plaque Assay was repeated twice with the picked plaques that were sent in for sequencing using different dilution steps.
The results of these Plaque Assays were promising but did not show single plaques. So in order to determine the phage titer the Plaque Assay was repeated one more time. On the next day plaque PCRs have been conducted with the new Primers again with 5 chosen plaques.

3. Test constructs for Riboswitch and CYP1A2:

Test constructs were designed in order to seperately test the Riboswitch of the AP and the CYP1A2 on the SP.
1) Riboswitch test construct:
Taking the backbone of the pSB4K5 Plasmid which contains RFP and cloning the Riboswitch right infront of the RFP replacing the previous RBS. This construct enables an easy readout when adding theophylline to see whether theophylline induces the Riboswitch.
2) CYP1A2 test construct:
Taking pSB1C3 and replacing the RFP by the CYP1A2 gene. When this construct is transformed into competent cells already containing the Theophylline PACE AP (contains the chaperone) and the substrate (caffein) is added, the activity of the enzyme can be checked.
So far the test construct have only been implemented in silico. Also all needed primers were designed and ordered.

Modeling

Calculations of medium consumption for iGEM goes green were performed and made interactive. A heatmap visualizes medium consumption, calculation of ideal turbidostat and lagoon sizes are possible. The user can annotate own experiments and compare them to others.
Heatmap with default values
Experimentation with elements of a Predcel model based on distributions instead of scalars were started. The idea is that a population of phages does not have one fitness between 0 and 1 but rather has individuals that have different fitness values. In this more complex model the concentrations have to be calculated for each phage fitness value, depending on the amount of phages that have that fitness value. The fitness distribution is changed by mutation and by selection. The first naive approach for mutation was programmed. It simply substracts a given percentage of the difference between the amount of a fitness value and the mean amount from the amount of a fitness value. Obviously this is oversimplified and will therefore be replaced by a model based on the idea that every sequence that mutates gets better or worse with normally distributed changes.

Optogenetics

SP Opto

The plaque assay of the SP trafo revealed plaques, but the negative control also. Nevertheless, 5 plaques were picked and a SP PCR with Q5 was performed in the same way as the day before. Analysis by gel electrophoresis exhibited bands of incorrect length. Probably a contamination of the TALEN phage occurred.

AP Opto

As the AP trafos exhibited no colonies on the Amp-plates, the Golden Gate was repeated using a slightly different cycler protocol:

Temperature [°C] Time [min] Cycles
37 3:00 15
16 4:00 15
16 5:00 1
65 2:00 1
Trafos in chemically competent Top 10 cells were performed and cells were plated on amp-plates.

Protein_Interaction

Week 35 and previous:

The idea of building a new genetic gene circuit for evolution of protein interactions in context of phage assisted continuous evolution (PACE) came up. Lots of different proteins were proofed whether they match the requirements of PACE evolution in iGEM context or not. Major focus at the beginning was the evolution of an anti-body-like protein. Several ideas and the respective cloning concepts were discussed. Even thus the imagined gene circuit should have worked several concerns let to further research.

Software

Performance of the Squeezenet Architecture - singlelabel 599

The model was run successfully on the old 599 classes dataset.
Parameters: lr = E-2, batchsize=64, epsilon=0.1
ROC
Precision

Performance of the Squeezenet Architecture - singlelabel 679

The model was run successfully on the 679 classes dataset.
Parameters: lr = E-5, batchsize=64, epsilon=0.1
ROC
Precision

Performance of the Squeezenet Architecture - Multilabel 1084

The model was run successfully on the 1084 GO-classes dataset.
Parameters: lr = E-3, batchsize=64, epsilon=0.1
ROC
Precision

Corrected datasets for missing classes, reworte eval() to enclude the whole validation set

  • Dataset 637, was missing 138 classes due to the min. length requirement in the DatasetGenerator class. The requirement was lowered to 175AA. Further the DatasetGenerator class was rewritten, to ensure to contain 5 samples from every class in the validation set.
  • the eval() function of DeeProtein was rewritten to perform the validaion on the whole validation set at given steps.
Performance on 679 classes with minlength 175:
lr=0.01, e=0.1, batchsize=64
ROC
Precision lr=0.001, e=0.1, batchsize=64
ROC
Precision Reinitialization with pretrained parameters and lower learning rate allowed finetuning of the classifier. Especially as the validation set is uniformally distributed (in contrast to the training set) the classifier can be considered as trained.

ROC/ACC/AUC-metrics

ROC and AUC was added to be calculated on the fly (after validation on the whole validation set.).

Training models on the embedded sequences

We generated batches from the word embeddings (dim=100, kmer-length=3) for the 679(EC) and the 1084 mulilabel network. However training proceeds much more slowly as the parametersize is 5 times the size of the one-hot network.

Multilabel-classification

In order to be able to perform multilabel classification, we rewrite the input pipeline (DatasetGenerator, BatchGenerator, TFrecordsgenerator) and generated two datasets with 339 and 1084 classes respectively. The considered labels were chosen solely based on their polulation. As the GO-term hierarchy follows a directed acyclic graph (DAG) we looked up all parent nodes for each leaf nodes and included the total set of annotations for each sequence. First models were run after extending the network for 2 convolutional and 2 1x1 layers on the 1084 classes dataset. Results were disenchanting.

Comparison of datasets

Total seqs after filtering (EC): 220488
Total seqs after filtering (GO): 235767

Datatype Dataset Samples % of filtered sequences considered % of total sequences considered
EC 679 165658 75.13 63.18
GO 1084 233386 98.99 89.00

Word2Vec

The new word2vec adaption we implemented last week was optimized by a few minor changes. Metadata of the embeddings can now be used to analyse the embeddings with tensorboard. Search for a single k-mer, finding the nearest neighbours, annotating frequencies or other properties to the points and search for groups of k-mers defined by regular expressions work. Embeddings of 3-mers and 4-mers in 50, 100 and 200 dimensions were calculated on whole swissprot and whole uniprot. Principle Component Analysis (PCA) was performed on a 100-dimensional embedding of 3-mers from swissprot and showed interesting properties. However, reduction of 100 dimensions to two or three may account for at least some of these.
Frequencies of 3-mers
Here the 3-mers marked darker are more frequent, those seem to cluster together.
All 3-mers containing a specific amino acid
For each amino acid all 3-mers containing it are marked in read. Some of those selection cluster together, others do not. This may also be due to the dimension reduction.
All 3-mers containing cysteine
All 3-mers containing cysteine are marked in red.
All 3-mers containing lysine
All 3-mers containing lysine are marked in red.
All 3-mers containing proline
All 3-mers containing proline are marked in red.
All 3-mers containing valine
All 3-mers containing valine are marked in red.
Especially cysteine, proline and lysine containing 3-mers cluster together, while the ones containing valine appear to be more equally distributed.
Clustering of k-mers based on amino acid content is a first hint, embeddings can be useful for input of neural networks. However their performance can probably only be measured indirectly by comparing the performance of a single architecture on different inputs.

Split_Polymerase

Dickinson-PACE 2.0:

Due to problems with the documentation in the awesome_sheet, we were using different strains for our PACE run. The results of two PCR reactions showed the abundance of a positive and a negative AP, which should not be a problem based on the findings of the properties of the MPs one week before. Follwing this further, we performed several days of PACE using this strain.
After seven days, plaque assays were performed and a plaque PCR was used for amplification of the N-terminal domain of the T7 RNAP. Even the sequencing of over 10 plaques did not yield any mutations in the essential domain. Due to this, we should go back to the old experimental design of Dickinson-PACE.

Dickinson-PACE:

We returned to the old experimental design, using MP1 instead of MP6 as well as the normal Dickinson strains without LuxAB behind gene III. For the organisation of our runs, we now use a better version of our google docs sheet. We are now documenting any problems and any results of our plaque assays and our phage detection PCRs. Plaque assays are implemented every 24 hours, while the Phage detection PCR is carried out every 12 hours. We could perform this PACE run for seven days, like planed before. After seven days, plaque assays were performed and single plaques were picked. The N-terminal domain of the T7 RNAP was amplified using PCR reaction. PCR product was send to GATC for sequencing analysis. This run did also not yield any mutations. On the one hand we therefore used an PCR on the phage supernatant, while on the other hand using single plaques as implemented before. Both methods showed no mutations, which is why we evaluated most of the parts. Especially the mutagenesis plasmid could be or should be the most probable problem. Calculations of the amount of arabinose used for induction in our PACE device and the original Dickinson-PACE leaded to the finding that we are using only 4 - 6,5 mM Arabinose instead of 30-40 mM which is used in Dickinson-PACE. This underlines the hypothesis that we are using not enough arabinose to induce the MP in a sufficient way.

PREDCEL:

First experiment:
Test the general induction of MP with different arabinose & glucose concentrations
Inoculate 20 different cultures with a bacteria from Stock
0, 5, 12.5 and 25 mM Glucose
0, 6, 12.5, 20 and 31 mM Arabinose
Read-out: time until the cultures reaches OD 0.6 - 0.8
Results:
only the cultures without arabinose reached the OD in sufficient time
all other cultures needed more than 30 hours for growing up Second experiment:
Let a culture with Glucose grow until it reaches OD 0.6
Let it cool on ice and split into 4 ml cultures with different concentrations of glucose (0 mM, 12,5 mM, 25 mM) and arabinose (0 mM, 6 mM, 25 mM)
Read-out: Incubate again and measure OD at different time points
Results:
differences in OD are small (probably due to high glucose concentrations)
Higher arabinose concentrations result in lower ODs Use PREDCEL for induction of mutagenesis and mutations
transformation of different Dickinson strains in bacteria containing MP6
even high selection pressure should be surmountable
the probability for directed mutation should be high
Afterwards - propagation of the phages having the essential mutations by strains containing no MP Random mutagenesis
Alternative:
Transformation of MP1, MP4 and MP6 in S2208 strain (= S2060 - PACE strain - and pJC175e - Plasmid which is used for phage propagation)
several implementations of PREDCEL
propagation of phages should lead to several mutations - “phage-library”
Afterwards - propagation of the phages having the essential mutations by strains containing no MP

Toolbox:

Most of the PCRs for amplification of the parts in our gibson-standard are done. Further steps will include the golden gate cloning for some parts and the integration into pSB1C3 for the iGEM standard. The amplifications which could not be implemented, will be carried out again or alternative primers will be ordered.

MinION Sequencing:

We implemented our first two sequencing runs with the MinION. According to the ONT Community, shearing can be performed using 26 gauge needles, instead of using the Covaris g-tubes. The library preparation was perfomred following the lambda control experiment protocol. Following this, the MinION was loaded according to Oxford Nanopores instructions. Each run was implemented for 10 h according to the supplied ONT script. Live basecalling was performed, yielding 2000 reads and 11 million events. The accuracy of both runs was approximately 65 %, which is not good enough for sequencing of our amplicons. Possible reasons for this fact, could be the age of the flow cells as well as deviations from the standard protocol, normally using AMPure beads instead of our DNA Dynabeads from ThermoFisher. For this reason, another experiment using a newer flow cell will be implemented. In this experiment, we will use PCR products from supernatant from different timepoints of Dickinson-PACE, in order to investigate the behavior of the MinION sequencing in this case. Further knowledge about the sequencing should be acquired by using the ONT community. In analogy to the Bt toxin paper, it should be possible to implement a sequencing workflow, enabling a reliable basecalling accuracy. A scheme of this workflow is provided as a picture.
Calendar Week 36

Cytochrome_Engineering

1. Cloning of AP Cytochrome PACE:

As the colony PCRs of this clone looked promising further expermients were carried out with it in the mean time. The aim of those experiments was to create AP variants with different ORIs. In order to do so, PCRs were performed to generate PCR products from Gibson overhang 5 to Gibson overhang 4 (3.6 kb fragment lacking the pBR ORI and AmpR). As those PCRs did not work, different conditions have been tested according to Phusion Flash polymerase trouble shooting protocol in order to optimize PCR conditions (different template concentrations, different cycler runs, touch down PCRs, Betain, DMSO, decreased amount of total cycles, different primer concetrations). After none of this worked either, further test PCRs with the template AP were conducted over the encoded chaperone protein, to see, whether the corresponding Gibson overhang was correct. The results indicated, that some of the APs carry the correct chaperone and some of the APs do not.
So the AP was send back in to GATC for sequencing with a different primer this time. Also the two other AP purifications were send in for sequencing to check, whether they fit better.
Sequencing results showed, that the AP clones from which the PCRs for the exchange of ORIs have been tried, is contaminated. The other two AP purifications showed no contamination as one of them showed clear results over the whole sequencing length, this one was used in all further experiments.
Since the aim of generating different AP variants could not be achieved yet, the next experiments should enable this aim. To increase the chance of success of the following PCRs, the 3.6 kb fragment was split in 2 smaller fragments (2.7 kb and 900 bp). Those fragments were meant to be amplified in separate PCR reactions. Only one of these PCR reactions succeeded and could be isolated from the agarose gel. For the other PCR the trouble shooting protocol was consulted again to optimize PCR conditions. Several different conditions (annealing temperature changed, annealing temperatur + primer concentration, cycler run have been tried also in combination but did not lead to the desired outcome.
Taking one step back, the PCR of the 900 bp fragment was repeated with a different AP template and finally worked. So the second fragment could be isolated from the agarose gel subsequently. The following Gibson assembly was conducted with both backbone fragments (Gibson overhang 5-3, Gibson overhang 3-4) and the two different ORI variants (p15A and pSC101) respectively. Afterwards the particular Gibon assembly products containing the different ORIs were transformed into chemically competent cells and plated on LB-Carbenicillin plates. To test the grown colonies, colony PCRs were performed which showed promising results so the corresponding growths were purified and sent in to GATC for sequencing. The sequencing results confirmed the desired outcome. So Glycerol stocks have been prepared from all different AP variants with pBR ORI, p15A ORI and pSC101 ORI.

2. Preparation of AP pBR ORI for PREDCEL:

In order to test the circuits in PREDCEL, the AP with pBR ORI was transformed into 1030 and 2060 electro competent cells and plated. The colonies of these trafos were picked and corresponding growths were inoculated. In parallel colony PCRs were performed to check the Trafo outcome. All of the picked clones showed excellent results. After the colony PCRs were repeated from the growths and did not show clear results, one of the other sequenced APs was transformed into 1030 and 2060 electro competent cells and into chemically competent cells as well. These Trafos have also been checked by colony PCR the next day. As this showed the desired results, glycerol stocks were prepared from the corresponding growths and the growths were plated on LB-Tet plates to select for the F-Plasmid, which is crucial for phage infection.
After the sequence of the other two AP variants (ORI p15A and pSC101) were confirmed, they were transformed into 1030 electro competent cells as well. The three different ORI strenghts (high copy-low copy) allow gradual adaption in the following PREDCEL experiments.

3. Test constructs for Riboswitch and CYP1A2:

Started to transform the pSB4K5 (iGEM Distribution Kit 2017, Plate 4, well 6H) and pSB1C3 (iGEM Distribution Kit 2017, Plate 4, well 4B). But still waiting for all primers necessary to build test constructs.

Modeling

Bugfixes and Discovery of new bugs. A first set of experiments was designed to adjust the model to the real behaviour of the experiment. Focus lies in the development of the phage titer under different conditions. The numeric PREDCEL model contained a mistake in the binding constant between phage and E. coli, k and was therefore rerun. Compared to the last results the corrected version has better results for longer time in the lagoons and higher starting concentrations of phage.
A first version of the iGEM goes green page with calculations for the medium consumption went online. This also was the first real application of in line LaTex as well as plotly.js on the wiki.
The numeric model with a fitness distribution instead of a single fitness value for a lagoon was expaneded by adding a mutation model based on normally distributed changes in fitness that is symmetric to the previous value for each value. The next improvement will probabbly be a parameter that allows for assymetric changes, as most mutations decrease fitness.

Optogenetics

SP Opto

The SP Golden Gate was repeated using a slightly different protocol: GG of Opto SP (30 µl Reaction):

reagents Volume [µl]
SP_BB (Purification ID: 669) 19.3
EL_222 (Purification ID: 520) 0.65 (of a 1:10 dilution)
BsaI 2
T4 Ligase 2
t4 Ligase Buffer 3
Cycler conditions were as follows:

Temperature [°C] Time [min] Cycles
37 3:00 25
16 4:00 25
50 5:00 1
80 5:00 1
About 5 µl of the GG product was transferred into a fresh 1.5-ml reaction tube, the remaining 25 µl were purified using the QIAquick PRC Purification Kit and eluted in 30 µl ddH20.
The GG products (unpurified and purified) were transformed into electrocompetent S1030 cells and recovered with SOC. Half of the culture was cultivated for 6 h, the other half was cultivated more than 14 h at 37 °C in a shaker. All four samples were used to perform a Blue Plaque Assay.
On the next day each sample revealed plaques. Plaques were picked and a SP detection PCR was performed.
The primers JM_068_rev and JM_069_fwd were used and a pipette tip containing phages of a single plaque was plunged into the 25 µl OneTaq reaction mix. The cycler conditions were as follows:

Phase Temperature [°C] Time [min] Cycles
Initial Denaturation 94 5:00 1
Denaturation 94 0:25 30
Annealing 51 0:40 30
Extension 68 2:00 30
Final Extension 72 5:00 1
Hold 10 1
Analysis by gel electrophoresis revealed the expected band of 1.9 kb (purified GG product, cultivated for more than 14 h). Four plaques of this plate were used to infect four 4 ml (2xYT + carbenicillin) E. coli culture (stock ID: 47) with an OD600 of 0.4. Cultures were cultivated for around 14 h and centrifuged at 8000 xg for 5 min to receive the phage containing supernatants.

AP Opto

The trafos of the APs revealed 2-5 colonies on each plate. A growth of two colonies of each AP was prepared in 4 ml 2xYT medium with carbenicillin. On the next day, a 2 ml glycerol stock (25 %) of each clone was prepared and the remaining cell suspension (3 ml) was purified using the QIAprep spin miniprep kit.
The AP plasmid purification was analyzed by sanger sequencing (GATC). The sequence was not as expected. It turns out, that for amplification of the AP_BB_light and AP_BB_dark the primers CG_005_rev and CG_008_rev were mixed-up. The purification of AP_BB_light is actually AP_BB_dark and the other way around. As the overhangs in the GG did not fit, it is hardly surprising that the transformation efficiency was that low.

Golden Gate assembly of APs:

The Golden Gate of the APs was repeated using the right fragments: GG of Opto APs (15 µl Reaction): AP_light

Reagents Volume [µl]
gIII_luxAB (Purification ID: 514) 1.9
AP_light_BB (Purification ID: 517) 0.5
pBLind 1 (of a 1:100 dilution)
BsaI 1
Bsa 1.5
T4 Ligase 1
t4 Ligase Buffer 1.5
AP_dark

Reagents Volume [µl]
gIII_luxAB (Purification ID: 514) 1.9
AP_dark_BB (Purification ID: 516) 0.5
pBLrep 1 (of a 1:100 dilution)
BsaI 1
Bsa 1.5
T4 Ligase 1
t4 Ligase Buffer 1.5
Cycler conditions were as follows:

Temperature [°C] Time [min] Cycles
37 3:00 25
16 4:00 25
50 5:00 1
80 5:00 1
The GG products were transformed into chemically competent Top 10 cells and plated on amp-plates.
On the next day each plate revealed many single colonies. Two colonies of each plate were picked and cultivated in 4 ml 2xYT with carbenicillin overnight. A glycerol stock (25 %) was prepared and the remaining 3 ml of the overnight culture was purified with the QIAprep spin miniprep kit.

Protein_Interaction

Week 36:

After some weeks of doing research to find a fitting protein to evolve in protein-interaction(PI)-context, one decided to evolve split variants of T7RNAP and dCas9 to reassemble more efficiently. Both of them are commonly used proteins in biotechnological context and thereby characterized well. Protein structures of the proteins are accessible and basic literature is available. Nevertheless, split variants have not been evolved to reassemble better, yet. Four matching split sizes for T7 RNAP and one for dCas9 were chosen and in consequence first primers have been designed to produce the respective APs and SPs. Primers were tested for primer-dimerization and secondary structures. Primer products will be fused using GoldenGate cloning based on BsaI and BsmBI restriction sites in primer extensions.

Software

Embedding

The embedding script did not really support restoring, that was fixed. The first embedding on uniprot that reached the second epoch was trained. Finding the best hyperparameters, especially the embedding dimension, and the length of the k-mers the embedding is based on, still remains. Checkpoints are saved as picklefiles now, making them accessible for applications without tensorflow.
Calendar Week 37

Cytochrome_Engineering

1. Preparation of APs for PREDCEL:

After all AP variants have been generated and their sequences were confirmed, electro competent cell strains were produced (according to our protocol for electro competent cells), each containing one ORI variant.
To start the first PREDCEL attempt the AP 1030 pBR ORI competent cells were chosen, as this is the strongest ORI and serves as a control to see whether something does happen or not. At the same time the AP 1030 pBR competent cells were transformed with MP6 to also generate a strain that produces mutations for later PREDCEL experiments.

2. First PREDCEL run:

The first PREDCEL run was conducted with AP 1030 competent cells (lacking MP6). The aim of this experiment was to see whether the addition of theophylline or caffeine to the AP strain (containing gene 3 and the chaperone protein for the CYP1A2) lead to an increasing phage titer, as gene 3 is crucial for phage propagation. Based on literature, theophylline and caffeine concentrations have been defined that were used in this experiment. So we chose to add 100 �M theophylline, 500 �M caffeine and 5mM caffeine. Also two different sample time points have been defined, after one hour and after 16 hours. The outcome of this experiment was checked by Plaque Assyas. Not to forget to mention, that besides the Theophylline SP also Dickinsons SP was used as a control.

Phase Temperature [°C] Time [min] Cycles
Initial Denaturation 94 5:00 1
Denaturation 94 0:25 30
Annealing 51 0:40 30
Extension 68 2:00 30
Final Extension 72 5:00 1
Hold 10 1
Both samples exhibited the expected 1.9 kb band on the agarose gel. The remaining PCR product was purified with the QIAquick PCR Purification Kit and sent to GATC for Sanger sequencing (using primer JM_068_rev).
The sequencing results revealed the right sequence! We have a infectious phage containing EL222. To determine the phage titer of the Opto_SP_Klon-8 supernatant, a Blue Plaque Assay was performed.
According to this plaque assay the titer of the Opto_SP_Klon-8 supernatant is 2.6*1015 PFU/ml.

AP Opto

The purified APs (2x AP_light, 2x AP-dark) were sent to GATC for Sanger sequencing (using primer JM_065_rev). The sequencing results revealed the right sequence! We have the right AP containing the promoters pBLind and pBLrep in front of gene III. Next, growths of AP_light-1 and AP_dark-1 were prepared in 4 ml 2xYT medium with carbenicillin. Plasmids were purified using QIAprep Spin Miniprep Kit. Trafos were performed transforming AP_light and AP_dark in S1030 and S1030 containing MP6. After one hour recovery in SOC, cells were plated on the respective plate.
Only the plate with AP_dark in S1030 cells contained cultures, which were very close and had to be plated on an other plate. Trafos were repeated with different strains, recovery media, plates, and purifications without success. A restriction double digest of the APs with HindIII-HF and SpeI-HF was performed and analyzed by gel electrophoresis. The expected bands were slightly visible on the gel, but also smear.

Phage_Propagation_and_Quantification

Simplification of Blue Plaque Assays

A 50 ml E. coli culture (ID:47) was cultivated in 2xYT with carbenicillin until an OD600 of 0.8. 3 ml of this culture were used to perform a Blue Plaque Assay. The remaining culture was put on ice for 10 min and was afterwards stored at 4 °C in the fridge. On the following five days, a Blue Plaque Assay was performed in the evening. Every Plaque Assay contained the 10-11 dilution of the Dickinson phage target 133 (1x1015 PFU/ml) in triplicates and a negative control. The plaque assays were performed using the E. coli culture from the fridge. Plaque assays worked for four days with a similar number of plaques. Therefore, it is confirmed that E. coli cultures for plaque assays can be stored in the fridge for at least four days.

Protein_Interaction

Week 37:

After primers had been delivered, DNA was resolved and diluted. First, all planned PCRs (see Awesomesheet for detailed conditions and protocol for general procedure for all named methods) were performed to get the primer products (PCR product IDs 357-375), that are needed for Golden Gate Cloning (starts with PCR ID 1321) using Q5 polymerase. Afterwards, PCR products were analyzed on gel and purified via gel extraction or PCR purification. PCRs that did not work in the beginning were optimized by altering annealing temperatures, temperature ramp rates and DMSO and Betain concentrations. Switching from Q5 to Fusion polymerase worked too in some cases. This week PCR products with IDs 358; 359; 373; 374; 375 could be produced with matching band length on gels.

Software

Embedding

Training of the embedding on Uniprot continues. First approaches to optimize the current multilabel classifier by adding parameters to it, either in depth or in width started.
The concept of a generator that models natural evolution was designed. There, a starting sequence is exchanged by randomly chosen point mutations, which are then scored by the classifier with a sigmoid output for each label. Sequences that score higher are more likely to be chosen to continue working with. Those will then be changed by point mutations again. In this scenario the selection is modeled by a function that uses the classifier. However it remains unclear, if the classifier supports gradual improvements in a direction that improves the function in reality and in the classifier.
Calendar Week 38

Cytochrome_Engineering

1. Big PREDCEL Run:

After the first PREDCEL run did not show any plaques in the Plaque Assay it was repeated. Both Plaque Assays then showed a positive result with blue plaques. The results have been confirmed by a conducted plaque PCR. A phage titer of 1,2 * 10^9 could be determined.
After the first PREDCEL Run was conducted with only one ORI variant (AP 1030 pBR ORI competent cells lacking MP6), the experiment was repeated using all three ORI variants under the same conditions (adding 100 �M theophylline, 500 �M caffeine and 5mM caffeine). Also the sample drawing was performed at the same time, after one hour and after 16 hours, for better comparison with the previous experiment. The outcome of this experiment was checked by Plaque Assyas. Again the Dickinsons SP was used as a control. In all three ORI experiments the PREDCEL run seemed successful, but the most promising results could be achieved with the pSC101 ORI. So for further experiments this ORI variant was chosen.
As in this big experiment the dilutions fort he control (adding neither theophylline nor caffeine) have been missed in the plaque Assay it was repeated fort he p15A ORI and the pSC101 ORI for both sample time points.
Another PREDCEL run was conducted adding a number of different caffeine concentrations (0 �M, 500 �M, 200 �M, 100 �M, 50 �M, 30 �M, 10 �M, 5 �M) to find the threshold where we can detect how much caffeine is needed to be converted into theophylline to have an effect on the phage propagation (evidence for activity of CyP1A2).

Riboswitch and CYP1A2 testing:

Next, we planned to use the HPLC to see how precisely we can detect theophylline and caffeine. So all necessary buffers and samples have been prepared.

Modeling

Equations for modeling the glucose concentration in PACE lagoon and in PALE lagoon were developed. The model that remains to be build can help finding the ideal glucose concentration in the medium. That means glucose concentration that are high enough to suppress the MPs before their induction, and low enough to allow for their induction.

Optogenetics

AP_dark_Ppsp-EL222-BR

To ensure that gIII expression is only activated after phage infection, a new Promoter: Ppsp-EL222-BR was developed. This is necessary to prevent cells from becoming resistance
towards phage infection as a consequence of the continuous gIII expression. The psp-Promoter and pBLrep were combined based on the paper "Negative selection and stringency modulation in phage-assisted continuous evolution".
The tet-binding region of the psp-tet promoter was exchanged by the EL222 binding region, which should lead to a gIII expression only in the dark and after phage infection. The psp-promoter was amplyfied using AP_destructase_SD2-luxAb (Purification ID: 536) as template and the primers CG_014 and CG_015 in a 50 µl Q5 reaction mix.
The cycler conditions were as follows:

Phase Temperature [°C] Time [min] Cycles
Initial Denaturation 98 0:30 1
Denaturation 98 0:10 30
Annealing 67 0:25 30
Extension 72 0:04 30
Final Extension 72 2:00 1
Hold 10 1
PCR product was analyzed by gel electrophoresis and the 170 bp fragment was purified with the QIAquick gel extraction kit (Purification ID: 852). Golden Gate of AP_dark_Ppsp_EL222-BR To generate the AP with gIII under control of Ppsp-EL222-BR, Golden Gate cloning was performed. AP_dark_Ppsp_EL222-BR

reagents Volume [µl]
gIII_luxAB (Purification ID: 514) 1.9
AP_dark_BB (Purification ID: 516) 0.5
Ppsp-EL222 (Purification ID: 857) 1.2
BsaI 1
BSA 1.5
T4 Ligase 1
t4 Ligase Buffer 1.5
Cycler conditions were as follows:

Temperature [°C] Time [min] Cycles
37 3:00 15
16 4:00 15
37 30:00 1
65 0:05 1
Transformation
Golden Gate product was transformed into Top 10 (chemically transformed) and cells were plated on a LB-agar-Amp plate.

AP_light and AP_dark

Glycerol stocks of AP_light and AP_dark were prepared from overnight cultures.

Phage_Propagation_and_Quantification

Purification

Theophylline SP and Dickinson SP were propagated and purified according to the NEB protocol "M13 Amplification". Phage titer was determined by Blue Plaque Assays. It was possible to reach phage titer of (1x1015 PFU/ml) and increase the number of phages per ml by a factor of 100
by step 4 to 6.

Quantification

A quantification method by spectrophotometry was tested according to "m13 bacteriophage production for large-scale applications", but as the phage titer differed more than 103 PFU/ml from plaque assay results the method was considered as unsuitable.

Protein_Interaction

Week 38:

PCRs for primer products that could not be produced in the previous week (PCR product IDs 357; 360 - 372) were repeated under altered conditions. Thereby, all the PCRs except PCR products IDs 358, 369 and 373 were performed successfully. As PCR products ID 369 (AP of Split T7RNAP backbone with GenIII) is not absolutely necessary because it can be replaced equivalently by PCR product ID 370 (AP of Split T7RNAP backbone with GenIII and YFP) all needed PCR products for T7 RNAP cloning were present to start GoldenGate assembly.
In contrast PCR products with IDs 358 and 373 are crucial for the dCas9 experiment and could not be produced in any way. Therefore focus moved on cloning T7 RNAP constructs alone rather than trying to manage both, the T7 RNAP and the dCas9 cloning, to work.

Software

PixelCNN

Adapting the pixelCNN from a 2-dimensional problem with three channels to a 1-dimensional one with 20 channels began. The network itself works but the logistic sampling function that generates the sequence in the end and the custom loss function are hard to adapt.

GAIA

The concept of a genetical artificial intelligent algorithm (GAIA) that is based on the classifier was developed. Starting from a given sequence it maximizes given labels and minimizes other given labels. At the same time the garbage score of the sequence is minimized. The model starts with one sequence but always works with B sequences. Therefore the starting sequence is in silico mutated B times to yield B different sequences. Those are then classified by DeeProtein and the outputs are used to calculate a sequence score. A given number of best sequences is selected and again mutated in silico to yield B sequences. This cycle of mutation and selection repeats and optimizes the starting sequence.
Additionally it is possible to also use recombination after the mutation to combine good subsequences
Calendar Week 39

Cytochrome_Engineering

1. PREDCEL:

The plaque assay with different caffeine concentrations did not work as the control (Dickinson phage) could not be detected. The results for the SP theophylline could not be clearly regarded as reliable then. To check for the SP Dickinson and to validate the results for SP theophylline the plaque assay was once more repeated. But unfortunately the negative control in this plaque assay showed a contamination. In order to get rid of the contamination, we performed several following plaque assays using completely new prepared materials.
Also new SP theophylline needed to be reproduced (from phage supernatant + stock 47 in 2YT + Carb)

2. Test constructs for Riboswitch and CYP1A2:

The next step, was a double transformation of the pSB4K5 (riboswitch) and pSB1C3 (CYP1A2) in Top 10 chemically competent cells. The colonies were picked to grow growths that are used for further testing. The growth was diluted 1:100 and caffeine or theophylline in different concentrations were added (5 mM, 500 �M, 50 �M, 5 �M, 0,5 �M) respectively. To ensure induction of the Lac Promotor IPTG was added (1:1000) as well as the required antibiotics (Kanamycin and Chloramphenicol).

3. HPLC Test:

Several samples have been run over the C18 column of the HPLC (theophylline 1 mg/ml, caffeine 1 mg/ml and a mix of theophylline and caffeine 1:1). The mix showed isolated peaks at about 5.92 min for theophylline and 9.23 min for caffeine. The test runs with the pure samples could not be conducted by implication (probably operating mistakes due to little experience with the HPLC).

Optogenetics

AP Opto

A colony of Top 10 cells, which was transformed with AP_dark_Ppsp_EL222-BR was picked and grown overnight in 4 ml 2xYT + carbenicillin. Out of this culture and an overnight culture of AP_light in S1030 cells the plasmids were purified according to the QIAGEN mini prep protocol and sent to GATC for sequencing. The DNA sequence of the AP_dark_Ppsp_EL222-BR out of Top 10 cells was correct, whereas the cells, which were thought to be AP_light (in S1030) exhibited the sequence of AP_dark. There must be a contamination or mix-up of AP_light and AP_dark. The purified plasmid AP_dark_Ppsp-EL222-BR was transformed in S1030 by electroporation and plated on Amp/Tet plates. Tet is necessary to make sure that the E. coli still carry the F-plasmid for phage infection.
On the next day, colonies were picked an grown in 4 ml 2xYT overnight. A glycerol stock of the overnight culture was prepared, and the remaining culture was used for plasmid purification. Plasmid purifications of AP_light, AP_dark and AP_dark_Ppsp_EL222-BR out of S1030 were send to GATC for sequencing. The sequence of AP_dark_Ppsp_EL222-BR was correct, but AP_dark as well as AP_light exhibited the sequence of AP_dark.

Protein_Interaction

Week 39:

GoldenGate cloning for all T7 RNAP constructs were performed and transformed into S1030 cells via electroporation. Transformation for phages worked properly, but AP transformation first did not show any cultures on plate to pick for glycerol stocks (plasmid IDs 305 + 306). Therefore, AP-transformation was repeated with 2060 cells, which worked.

Software

Embedding

Embedding on Uniprot is calculated over 15 epochs, vectors that are close to each other in the high dimensional space seem to be similar on sequence level on first sight.
Genetic artificial intelligent algorithm is implemented. Features are mutation of a given maximum of residues, mutation of specific residues. Those are inserted in the model in small letters while the rest of the sequence is written in capital letters. Furthermore weighting arbitrary numbers of goal- and avoid- GO-terms are possible. The variance of the given classes, a measure for the authenticity of a protein, is substracted from the overall score. Recombination of sequences was not implemented as it corrupts maintaining less then the maximum allowed mutations in a sequence.
Calendar Week 40

Cytochrome_Engineering

1. PRDECEL:

After the plaque assay with different caffeine concentrations did not work and also the repeated experiment, it needed to be conducted from the beginnig. Therfore, new phages had to be reproduced to have a sufficient amount of phages to inoculate the new PREDCEL culture. The phage reproduction worked well, a titer of 1*10^10 could be achieved. This could then be used to inoculate the new PREDCEL culture in order to repeat the whole process. The outcome was checked by plaque assay afterwards. The plaque assay worked well this time, but no significant difference between caffeine (in different concentrations) and no caffeine could be detected. So we came to the conclusion, that in order to meassure the CYP1A2 activity another assay has to be planned.
Furhtermore, the first PREDCEL run with passaging was conducted (in triplicates, and 100 �M and 1000 �M theophylline were used). As most of the reproduced phages were used up for this experiment, new phages had to be reproduced again. The titer was then checked by plaque assay as well as the outcome of the PRDCEL run. Checking all triplicates from all sample timepoints would have been to laborious for a frist check, so only the first, the second and the last time point were checked by plaque assay. There one could see, that the last timepoint did not show plaques, which indicates that the five passages could not be conducted successfully. Nevertheless, some of the passages might be successful and the samples inbetween the second and the last timepoint were checked by plaque assay afterwards.

2. Transformmation:

Previously electrocompetent cells containing the AP and all three ORI variants have been generated. Those were now transformed with the expression plasmid containing the CYP1A2 gene (which was used in further testings before). Right after recovering the transformed cells, they were plated on LB-chloramphenicol- ampicillin-plates.
All three variants could be successfully transformed and colonies were picked respectively in order to produce glycerol stocks. The next day an experiment was started to check for the CYP1A2's correct tertiary structure. To do so, growths have been inoculated from the glycerol stocks and the cells containing the CYP1A2 were then lysed by ultrasonic waves and centrifuged. Then 50 �l Disodiumdithionide (10 mg/ml) were added to 700 �l of the supernatant. The Disodiumdithionide reduces the heme group, which is essential for the catalytic activity of the CYP1A2. The heme group can only be bound, when the proper three dimensional structure of the CYP1A2 is formed. So indirectly, this experiment can indicate, whether our CYP1A2 is in its active conformation or not.

3. HPLC:

As the HPLC test with a mixture of caffeine and theophylline (1:1) worked but the pure solutions of theophylline and caffeine did not show any peaks, the experiment was repeated several times. We assumed that operating mistakes or wrong settings are the reason for the negative test results.
But as the experiments have been run under the same conditions this week again, there is probably another source of mistakes. Only in a few cases peaks could be detected, but then distinct peaks were detected.
Outlook: probably an assay to meassure CYP1A2 activity.

Modeling

A simple model for calculating the glucose concentration in a turbidostat from the glucose concentration in the medium and vice versa was made. A tool with a short summary of the math behind it was uploaded to the tools page of the wiki.
The same model was adapted for a flask instead of a turbidostat, changing it from a steady state model to one that integrates over time. Bugfixes remain to be done.

Optogenetics

AP Opto

As the S1030 cells, containing AP_light, were contaminated or mixed-up with AP_dark a new AP_light strain
had to be prepared. The purifications of AP_light and AP_dark, which were already checked
by sequencing, were used to transform new S1030 cells by electroporation and plated on Amp/Tet plates. On the next day a
colony was picked and grown in 4 ml 2xYT + carbenicillin overnight. A glycerol stock was prepared, and the remaining culture was used to isolate and purify the plasmid.
The plasmid purifications were sent to GATC for sequencing. The sequencing results proved that the cells carry the right promoter in front of gIII.

AP Testing

To test if the SP_Opto phage can propagate in S1030 containing AP_light and AP_dark,
cultures with the respective cells were grown in 20 ml 2xYT + carbenicillin to an OD600 of 0.6 and infected with either SP_Opto of the Dickinson phage target_133
and grown in the dark and under blue light After three and six hours samples were taken. A plaque assay was
performed to check if the Opto phages can propagate, but due to a phage contamination of the S2208 cells the plaque assay failed. Three additional plaque assays failed, as the glycerol stock and pipettes were contaminated.

Protein_Interaction

Week 40:

This week was mainly about checking the phages to have the right SP inserts, which came along with lots of troubleshooting. Test PCRs did not or not only show the expected band lengths. As a band appeared, which length matches phage backbone size without insert and as well PCR product of Dickinson phage, a phage construct used for validation of our PACE device, contamination of samples could not be excluded. Therefore at the beginning PCR was optimized and transformation was repeated. As problem could not be solved, transformation was repeated a second time and new primers were designed showing significant differences in product length. New PCR needed optimization as well, but finally worked well and robustly, so it became standard PCR for contamination detection of PI-constructs.

Software

GAIA

GAIA was expanded by allowing multiple mutations before selection takes place, a linear decay of the number of mutations before a round of selection is applied. Results look slightly better compared to one mutation before selection. With this mode epistatic mutations can be found, that means mutations that are beneficial together but decrease the score of a protein when occuring alone.
Calendar Week 41

Cytochrome_Engineering

1. PREDCEL

Another plaque assay was conducted with the samples of the firste passaging experiment to find out whether we are able to propagate phages over several passaging steps. The outcome showed that 3 passages have been conducted successfully. Subsequently, all triplicates of the last timepoint were screened in order to see if a significant difference between Theo100 and Theo1000 and no Theo can be detected (-> AP works, more phages when theophylline is added, but storing samples in the refigerator reults in lower phage titers).
As only three of five passages have been successfull, a next PREDCEL run was started where more pahges were transferred (1ml instead of 0.5 ml) to see whether this can lead to more successful passaging steps. This time Theo100 was tested along with different caffeine concetrations. The samples have again been checked by plaque assay. Unfortunately, there was a phage contamination in all applied dilutions, so no conclusions could be drawn.
So another PREDCEL was started using the same conditions plus adding 5 mmol of 5-Aminolevunilic acid, a precursor of the heme molecule which displays an important cofactor for the CYP1A2 (this was added, as in the mean time CYP1A2 testing showed visible colour change when adding the substance, indicating higher activity of the CYP1A2). The plaque assay of this run confirmed that adding theophylline leads to faster phage propagation and showed that caffeine concentrations between 5-500 �M can also impact phage propagation.

2. CYP1A2 Test

To find out whether our CYP1A2 is functional or not, we repeated the dithionide experiment as no sufficient peak could be detected the first time. Small changes took place: a change of buffer was performed before ultrasonifying the cells. Also a dithionide solution was directly prapared befor meassuring the optical density at 550 nm. As this did not lead to a colour reaction, it was repeated once again adding 5-Aminolevunilic acid, which is a precursor of the heme molecule which in return displays an important factor for CYP1A2 activity. We added 5-Aminolevunilic acid, as we assumed, that mybe ther would not be sufficient amounts in e.coli available. This time a colour change could already be detected by eye, as the colour of the growths was red/brown. The optical density meassurement at 550 nm then showed a peak, indicating, that our CYP1A2 is functional.

Modeling

Logistic growth was implemented for glucose modeling to be able to model larger time scales with E. coli leaving exponential phase. Plotting was improved.

Optogenetics

Weekly summary 09.10-15.10.2017 CG

Opto PACE

AP Opto

To prevent phage contamination all instruments and tools were cleaned, all solutions were prepared freshly and
filter tips were used more consequently. It turned out that there was a phage contamination in the glycerol stock of AP_dark_Ppsp-EL222 containing S1030 cells. As these cell
have a psp-promoter, phages which don't contain EL222 can propagate with and without blue light.
A new colony was picked from the plate, but this colony was also already contaminated with phages. Therefore, AP_dark_Ppsp_EL222-BR had to be transformed again into S1030 cells.
The same AP_dark_Ppsp_EL222-BR purification was used to transform S1030 cells by electroporation and cells were plated on Amp/Tet-plates. On the next day, a colony was picked and grown in 4 ml 2xYT + carbenicillin overnight. A glycerol stock was prepared, and the remaining culture was used to isolate and purify the plasmid.
The plasmid purifications were sent to GATC for sequencing. The sequencing results proved that the cells carry the right promoter in front of gIII. AP_dark_Ppsp_EL222-BR is further referred as AP_dark.

AP Testing

Cultures of S1030 cells containing AP_light and AP_dark were prepared and grown to an OD600 of 0.6 in 2xYT with carbenicillin. Then, this culture was split into 20 ml cultures and infected with 10 µl of either SP_Opto or Dickinson phage target_133 (titer: (~1015 PFU/ml) according to the following table:

Strain Phage Condition
AP_light SP_Opto blue light
AP_light SP_target_133 blue light
AP_dark SP_Opto blue light
AP_dark SP_target_133 blue light
AP_light SP_Opto dark
AP_light SP_target_133 dark
AP_dark SP_Opto dark
AP_dark SP_target_133 dark
After infection the culture was mixed and samples were taken. Cultures were incubated for one hour at 37 °C in the shaker, half of them in the dark, the other half was illuminated with blue light pulse 25 % (15 sec. on, 45 sec. off). Afterwards samples were taken. All sampels were centrifuged at 8000 xg for 3 min in a table centrifuge. The supernatant was transferred into a fresh reaction tube and stored at 4 °C in the fridge. These samples were used to perform a plaque assay using 10 µl phage supernatant per quarter of the plate. Plaque assays was performed according to the Blue Plaque Assay protocol. On the next day, plaques were counted to determine the phage titer:

Strain Phage Time [h] Condition Phage titer [PFU/ml]
AP_light SP_Opto 0 - 0.7x107 PFU/ml
AP_light SP_target_133 0 - 0.6x107 PFU/ml
AP_light SP_Opto 1 blue light 2.2x107 PFU/ml
AP_light SP_target_133 1 blue light 1.4x107 PFU/ml
AP_light SP_Opto 1 dark 2,3x107 PFU/ml
AP_light SP_target_133 1 dark 1.9x107 PFU/ml
AP_dark SP_Opto 0 - 1.1x107 PFU/ml
AP_dark SP_target_133 0 - 2.1x107 PFU/ml
AP_dark SP_Opto 1 blue light 2.3x107 PFU/ml
AP_dark SP_target_133 1 blue light 3.6x109 PFU/ml
AP_dark SP_Opto 1 dark 6.2x107 PFU/ml
AP_dark SP_target_133 1 dark 5.7x109 PFU/ml
There was the right tendency for AP_light and AP_dark after one hour. For further experiments the incubation time should be increased. The Dickinson phage target_133 exhibited higher phage titer than SP_Opto. For AP_dark this is not surprising as this phage don't contain EL222 and therefore activates the psp-promoter. But it should not activate the promoter of AP_light, pBLind. As it was suspected that the Dickinson phage is superior in its propagation cycle due to e.g. packing ability, a phage containing a non-binding EL222 should be created.

Protein_Interaction

Week 41:

Parts of phage genome with insert sequence were amplified, sequenced and thereby approved after two to three sequencings. While SP sequencing worked AP sequencing did not show clear results in three sequencings, which is why GoldenGate on AP related PCR products was performed again. Then, APs were transformed and purified, again. As own APs were not approved by sequencings Jin69 plasmid (plasmid ID 17) was used for phage propagation testing over 24 hours.

Software

GAIA

GAIA was used to generate mutations in a glucoronidase sequence, beneficial for galactosidase activity. A subset of the generated sequences was ordered as oligonucleotides for kinetic assays. In order to show that the score from the used classifier correlates with the real-world activity, a set of less active beta-lactamase sequences was predicted by gaia. Those were also bought as oligonucleotides and will be tested for their activity.
A randomization function was implemented to examine how the different scores develop when the number of mutations increases. The results look sigmoid, showing almost no change in the score for the first mutations, then rapidly dropping below the threshold for positive classification, after which they only decrease slowly.
A combination function to examine all possible combinations of a set of mutations was implemented in order to aid the decisions on which sequences to test in the lab.
Calendar Week 42

Modeling

Modeling the phage titer based on a distribution of fitness was improved. Mutation modeled by a skewed normal distribution was enabled.
The model was compared to phage titer data from the lab and performed well. However, improvements in the scaling are needed. All available interactive tools were unified and put into the wiki layout.

Optogenetics

Weekly summary 16.10-22.10.2017 CG

Opto PACE

SP_Opto_EL222_non-binding

A primer was ordered that binds within the helix-turn-helix motive of EL222 and should therefore
amplify a truncated non-binding variant of EL222.
EL222_non-binding was amplified using the gblock EL222 expression cassette (PCR Product ID: 242) as template and the primers CG_012 and CG_016 in a 25 µl Phusion Flash reaction mix.
The cycler conditions were as follows:

Phase Temperature [°C] Time [min] Cycles
Initial Denaturation 98 0:10 1
Denaturation 98 0:01 25
Annealing 72 0:05 25
Extension 72 0:09 25
Final Extension 72 1:00 1
Hold 10 1
PCR product was analyzed by gel electrophoresis and the 548 bp fragment was purified with the QIAquick gel extraction kit (Purification ID: 1006). Golden Gate of SP_Opto_EL222_non-binding To generate an SP with a non-binding EL222 variant Golden Gate cloning was performed. SP_Opto_EL222_non-binding

reagents Volume [µl]
SP_Opto_EL222_non-binding (Purification ID: 1006) 3.56
SP_BB (Purification ID: 1007) 3.33
BsaI 1
BSA 1.5
T4 Ligase 1
T4 Ligase Buffer 1.5
Cycler conditions were as follows:

Temperature [°C] Time [min] Cycles
37 3:00 15
16 4:00 15
37 30:00 1
65 0:05 1
Golden Gate Product was transformed into S2208 cells by electroporation and were grown for 16 h in SOC after recovery.
Culture was centrifuged at 8000 xg for 3 min and phage supernatant was transferred into a fresh tube. This supernatant was then used to perform a Blue Plaque Assay.
On the next day a plaque was picked and used for an insert PCR (JM_068, JM_069) as well as for infection of a new culture of S2208 cells (OD600 of 0.6). The PCR and subsequent gel electrophoresis exhibited a fragment of the expected length. The infected overnight culture was centrifuged at 8000 xg for 3 min and phage supernatant was transferred into a fresh tube stored at 4 °C.

Protein_Interaction

Week 42:

By plaque assay performance phage propagation on Jin69 was proved to work, so gene circuit can be used in PACE and/or PREDCEL. After several attempts of positive AP sequencing AP were finally proven right. Related stocks were proven not to be contaminated by contamination detection PCR. S1030 Cells containing approved APs were made electrocompetent and MP1 and MP4 were transformed into these cells separately by electroporation. For proper and complete sequencing of T7 RNAP to detect mutations new sequencing primers were designed and ordered. Furthermore to exchange SD8 RBS sequence with weaker sd8 RBS sequence GoldenGate primers were bought and used to perform GoldenGate cloning.

Software

GAIA

For further insight in properties of a sequence, a function was implemented, that scores every possible single mutant and plots changes in the score. Potentially active site and other crucial sequence regions can be predicted by this.
Calendar Week 43

Optogenetics

SP_Opto_EL222_non-binding

A SP_insert PCR (JM_068, JM_069) was performed and sent to GATC for sequencing.
The sequencing results proved that the phage contains a truncated version of EL222.
This phage is further used as control for AP_testing.

AP_testing

Cultures of AP_light and AP_dark were prepared and grown to an OD600 of 0.6. A sample was taken to
prove that the cultures were not already contaminated with phages. Then, the two cultures were split into 10 ml cultures and infected with 10 µl of either SP_Opto or SP_Opto_EL222_non-binding (titer: (~1011 PFU/ml) to reach a MOI of 1.
All strain, phage and cultivation conditions were prepared in duplicates. Cultures were cultivated for 3 h either in the dark or under blue light (25 %: 15 sec. on, 45 sec. off). After 3 h 2 ml of each culture was centrifuged at 6000 xg for 3 min and supernatant was transferred into a fresh reaction tube. 1 ml of this supernatant was used to infect a fresh E. coli culture (OD600 ~ 0.6)
which was again cultivated for 3 h under the same conditions as before. This transfer step was repeated twice and final samples of T=4 were analyzed by a Blue Plaque assay (T=4: after 4 x 3 h of cultivation, 3 transfer cycles). On the next day plaques were counted and phage titers were determined:

Strain Phage Transfer cycle Condition Duplicate Phage titer [PFU/ml]
AP_light SP_Opto 0 - 6x106 PFU/ml
AP_light SP_Opto_EL222_non-binding 0 - I 10x106 PFU/ml
AP_light SP_Opto 3 blue light I 7x104 PFU/ml
AP_light SP_Opto 3 blue light II 2x104 PFU/ml
AP_light SP_Opto_EL222_non-binding 3 blue light I 1x104 PFU/ml
AP_light SP_Opto_EL222_non-binding 3 blue light II 1x104 PFU/ml
AP_light SP_Opto 3 dark I 4x104 PFU/ml
AP_light SP_Opto 3 dark II 1x104 PFU/ml
AP_light SP_Opto_EL222_non-binding 3 dark I 4x104 PFU/ml
AP_light SP_Opto_EL222_non-binding 3 dark II 4x104 PFU/ml

Strain Phage Transfer cycle Condition I Phage titer [PFU/ml]
AP_dark SP_Opto 0 - I 3x106 PFU/ml
AP_dark SP_Opto_EL222_non-binding 0 - I 8x106 PFU/ml
AP_dark SP_Opto 3 blue light I more than 1014 PFU/ml
AP_dark SP_Opto 3 blue light II 5.3x1013 PFU/ml
AP_dark SP_Opto_EL222_non-binding 3 blue light I more than 1014 PFU/ml
AP_dark SP_Opto_EL222_non-binding 3 blue light II more than 1014 PFU/ml
AP_dark SP_Opto 3 dark I 9x1012 PFU/ml
AP_dark SP_Opto 3 dark II more than 1014 PFU/ml
AP_dark SP_Opto_EL222_non-binding 3 dark I more than 1014 PFU/ml
AP_dark SP_Opto_EL222_non-binding 3 dark II more than 1014 PFU/ml

Protein_Interaction

Week 43:

In Week 43 strain with AP and MPs were used to perform 2 PACE runs and several PREDCEL rounds according to the respective protocols. Besides that RBS exchange was finished and plasmid transformation was done, again, by electroporation.

Software

GAIA

The usability of the script was improved by adding supports for passing parameters via flags. Plotting heatmaps of scores with the sequence on x-axis and the amino acids to which exchanges were scored on the y-axis was improved.