Team:Waseda/Project

Project - iGEM 2020 Team:Waseda

Project

zombie model

We based the preparation of E. coli extracts for cell-free reactions on a paper by Jewett [Kwon et.al 2015]. The basic protocol is as follows.

In this COVID-19 crisis, we, iGEM Waseda, tried to attract high school students and general public to synthetic biology through three tactics. First, since cell-free system is relatively easy to operate, we wanted to provide an opportunity high school students and general public to experience synthetic biology hands-on. We thus proceeded with our iGEM project adopting the cell-free system. This was an advantage for us to continue with our wet lab experiment in this COVID-19 pandemic situation. Second, we have developed a comprehensible tool and given an interesting story to modeling. Modeling plays an essential role in synthetic biology, but it is difficult to explain accurately to the general public. These stories and tools will help them to understand and enjoy synthetic biology. Last but not least, the theme we are dealing with must not only be interesting as a story, but also biologically. This year, we thus decided to study decision making by the toggle switch gene circuit, which is an important basis in developmental biology and regenerative medicine; in those research fields, cell-fate decision making is quite an important notion.

Fig.2-1-1 Benefit of cell-free system
Fig.2-1-1 Benefit of cell-free system
Fig.2-1-2 Image ou our project
Fig.2-1-2 Image ou our projec
Fig.2-1-3 Image of landscape of cell fate
Fig.2-1-3 Image of landscape of cell fate

The story we developed was the war between zombies and Samurai. Both Zombie and Samurai are essentially human beings, but their states switch by some factor (Fig.2-1-4 top). Zombies are able to make Samurai join their group by releasing communication molecules (Fig.2-1-4 middle). On the other hand, Samurai also try to save the companion. At the end of the battle, the winner side is defined by the concentration of communication molecules in the field (Fig.2-1-4 bottom).

Fig.2-1-4 State switching with factors
Fig.2-1-4 State switching with factors
Fig.2-1-5 our genetic circuit
Fig.2-1-5 our genetic circuit

A genetic circuit for the Zombies and Samurai war (or cell-fate decision system) consists of two sub-systems: toggle switch and cell-cell communication module. The two states, Samurai and Zombie, are achieved by two stable states of a toggle switch circuit composed of two repressors each of which inhibiting gene expression of the other (Fig.2-1-6 left). A Zombie-state cell produces not only Zombie repressor but also Zombie communication molecule produced by an enzyme encoded downstream of the repressor. The increase of the Zombie cells increase the concentration of Zombie communication molecule designed to convert a Samurai cell. For that conversion mechanism, all of the cells, even a cell in a Samurai state has Zombie repressor gene regulated by Zombie communication molecule. For the symmetry of the scenario, all of the cells have genes for Samurai repressor and enzyme producing Samurai communication molecule.

To demonstrate this model, we solved the ordinary differential equations containing Hill repression and activation terms. The ideal result for the war between Samurai and Zombie is in the right panel of Fig.2-1-6 If there are more samurai cells in the initial state, almost all cells will be samurai at the end point, and vice versa. This has previously been achieved with ordinary toggle switch genetic circuits without cell-cell communication.

Before modeling of in vitro systems with cell-cell communication, we established differential-equation based modeling for the same circuit which works in living E. Coli cells (Fig.2-1-7) because of accumulated cases for such in vivo models. In this modeling, we assumed the behavior of E. coli in vivo, and used parameters adjusted somewhat to behave as designed. From the appearance of the nullcline (Fig.2-1-7), the equilibrium points are arranged in the same way as a normal toggle switch. Thus, once our parameter search found the adjusted parameters for cell-free systems, we were able to realize bistable system.

Fig.2-1-7 the equation of our modeling and nullcline
Fig.2-1-7 the equation of our modeling and nullcline

Fig.2-1-8 illustrates the behavior of our gene circuit. When the zombie state cell is induced with a samurai signal, the cell fate will change depending on the size of the signal. If the zombie receives a large number of samurai signals, the zombie will not be able to withstand the signal and will turn into a samurai. However, if samurai signal received by the zombie state cell is not enough, the zombie can revert back to being a zombie again, without completely transforming into a samurai. As a result, we were able to create two final states simply by changing the initial concentration of communication molecules.

Fig.2-1-8 The concept behavior of our gene circuit
Fig.2-1-8 The concept behavior of our gene circuit

In contrast to cellular system, a cell-free system lacks cell growth and has much smaller degradation- dilution term, which lead to longer operating time. By using parameters from reference (Maurizi MR,1992), we thus estimated operating time for state-switch between zombie and samurai in our cell-free system.

Fig.2-1-9 Degradation ratio and time of transition between equilibrium points
Fig.2-1-9 Degradation ratio and time of transition between equilibrium points

Fig.2-1-9 shows the the relationship between this degradation rate and operation time for state switching. The blue point is about normal protein degradation, and the green point shows the case of the LVA degradation tag are introduced on the protein. Considering a system based only on normal protein degradation(1%/hour =1.6*10^-4/min )of a cell-free system without growth of volume, it takes too much time (approximatery 2x104 min = 300 hours) to operate a state-switch between Samurai and Zombie. To solve this problem, we first investigated the size of the required degradation term. Operation time for the state switch was estimated from the time that high concentration of Samurai repressors at one equilibrium point are reduced by degradation to the repressor concentration at another equilibrium point. Although we did not consider the term of protein synthesis in this calculation, there is no problem because this assumption without protein synthesis shows shorter operation time than a case with protein synthesis.

Although cell-free system is selected as our project platform for implementation of the zombie vs. samurai scenario using a toggle switch circuit, cell-growth rate is known to be important for in vivo toggle switch. Its behavior can be predicted by a mathematical modeling which include a dilution-dependent decrease of concentration of repressors. For cell-free system, our modeling described above indeed have shown that cell-free toggle switch will not work if the system has low degradation rate.

In order to increase the degradation rate for specific proteins in our cell-free system, we decided to add LVA degradation tag to the repressor proteins. LVA degradation tag is a 9~11 amino acids sequence attached to the end of a protein. The tagged proteins are well degraded by a certain combination of proteins: clpX, clpP and sspB.

Fig.2-1-6 The ideal results of our gene circuit
Fig.2-1-6 The ideal results of our gene circuit

To demonstrate this model, we solved the ordinary differential equations containing Hill repression and activation terms. The ideal result for the war between Samurai and Zombie is in the right panel of Fig.2-1-6 If there are more samurai cells in the initial state, almost all cells will be samurai at the end point, and vice versa. This has previously been achieved with ordinary toggle switch genetic circuits without cell-cell communication.

Before modeling of in vitro systems with cell-cell communication, we established differential-equation based modeling for the same circuit which works in living E. Coli cells (Fig.2-1-7) because of accumulated cases for such in vivo models. In this modeling, we assumed the behavior of E. coli in vivo, and used parameters adjusted somewhat to behave as designed. From the appearance of the nullcline (Fig.2-1-7), the equilibrium points are arranged in the same way as a normal toggle switch. Thus, once our parameter search found the adjusted parameters for cell-free systems, we were able to realize bistable system.

Fig.2-1-7 the equation of our modeling and nullcline
Fig.2-1-7 the equation of our modeling and nullcline

Fig.2-1-8 illustrates the behavior of our gene circuit. When the zombie state cell is induced with a samurai signal, the cell fate will change depending on the size of the signal. If the zombie receives a large number of samurai signals, the zombie will not be able to withstand the signal and will turn into a samurai. However, if samurai signal received by the zombie state cell is not enough, the zombie can revert back to being a zombie again, without completely transforming into a samurai. As a result, we were able to create two final states simply by changing the initial concentration of communication molecules.

Fig.2-1-8 The concept behavior of our gene circuit
Fig.2-1-8 The concept behavior of our gene circuit

In contrast to cellular system, a cell-free system lacks cell growth and has much smaller degradation- dilution term, which lead to longer operating time. By using parameters from reference (Maurizi MR,1992), we thus estimated operating time for state-switch between zombie and samurai in our cell-free system.

Fig.2-1-9 Degradation ratio and time of transition between equilibrium points
Fig.2-1-9 Degradation ratio and time of transition between equilibrium points

Fig.2-1-9 shows the the relationship between this degradation rate and operation time for state switching. The blue point is about normal protein degradation, and the green point shows the case of the LVA degradation tag are introduced on the protein.

Considering a system based only on normal protein degradation(1%/hour =1.6*10^-4/min )of a cell-free system without growth of volume, it takes too much time (approximatery 2x104 min = 300 hours) to operate a state-switch between Samurai and Zombie. To solve this problem, we first investigated the size of the required degradation term.

Operation time for the state switch was estimated from the time that high concentration of Samurai repressors at one equilibrium point are reduced by degradation to the repressor concentration at another equilibrium point. Although we did not consider the term of protein synthesis in this calculation, there is no problem because this assumption without protein synthesis shows shorter operation time than a case with protein synthesis.

Although cell-free system is selected as our project platform for implementation of the zombie vs. samurai scenario using a toggle switch circuit, cell-growth rate is known to be important for in vivo toggle switch. Its behavior can be predicted by a mathematical modeling which include a dilution-dependent decrease of concentration of repressors. For cell-free system, our modeling described above indeed have shown that cell-free toggle switch will not work if the system has low degradation rate.

In order to increase the degradation rate for specific proteins in our cell-free system, we decided to add LVA degradation tag to the repressor proteins. LVA degradation tag is a 9~11 amino acids sequence attached to the end of a protein. The tagged proteins are well degraded by a certain combination of proteins: clpX, clpP and sspB.

To see the improvement in degradation effect by LVA degradation tagging, we constructed an improved part (BBa_K3580003) by a modification of an existing part: Plux/tet-GFP( BBa_K934025). Reduction of the concentration of proteins by both degradation/ dilution and increase of the concentration by production play an important role in the behavior of gene circuit. Although protein concentration decrease by dilution caused by the growth of the cell, cell-free systems can’t because its volume is fixed. Therefore, an improvement is necessary to incorporate degradation into the cell-free system.

In order to compare in vivo activity of these parts (BBa_K3580003 and BBa_K934025)), we first measured the fluorescence of GFP 240 minutes after the start of induction. The fluorescence of tagged GFP (BBa_K3580003 ) was lower than that of normal GFP (BBa_K934025)) at 240 min point (fig 2-1-10). Although GFP is a stable protein with a β-barrel and difficult to be degraded, this result shows that tagged GFP was successfully degraded as we planned.

Fig 2-1-10 ssrA degradation tag assay in vivo
Fig 2-1-10 ssrA degradation tag assay in vivo

Then, we compared the fluorescence of GFP in a cell-free system which consisted of cruedextract of E.coli containing luxR protein (Fig 2-1-11). Because of the programmed degradation, the fluorescence of tagged GFP(BBa_K3580003) showed slight signal nearly equal to a negative control where neither the template GFP DNA nor the inducer AHL existed. The results show that LVA degradation tagged protein can be degraded exceptionally both in vivo and in vitro. We will be able to enjoy state switching in a test tube. So we further promoted our modeling.

Fig 2-1-11 ssrA degradationtag assay in vitro
Fig 2-1-11 ssrA degradationtag assay in vitro


Fig.2-1-12 Basic war in cell-free system
Fig.2-1-12 Basic war in cell-free system

Our simulation showed fate of the battle between zombies and samurais in the test tube Fig.2-1-12 Depending on initial concentration of components, we can prepare a zombie test tube and a samurai test tube which are stable state. Then, by mixing zombie and samurai cell-free solutions in various ratios, the battle started. When zombie state cell-free solution and samurai state cell-free solution were mixed at a ratio of 2:8, the mixed solution was stable in the samurai state. This means that the samurais have beaten the zombies. However, a slight change of the ratio completely turned their fate. When the zombie and samurai solutions were mixed at a ratio of 3:7, the mixed solution was stable in the zombie state and the zombies willdestroyed the Samurai.

Cross talk assay in vitro

When we construct and model the “zombie vs samurai” system, a potential problem was crosstalk, by which a promoter in a quorum sensing (QS) systems is activated unexpectedly by another type of QS components. Consequently, we some experiments to reveal which combination of QS systems have less crosstalk ability in a cell-free system.

Fig 2-1-13  Crosstalk in the Zombie vs Samurai system
Fig 2-1-13 Crosstalk in the Zombie vs Samurai system
Fig 2-1-14  the possible crosstalks
Fig 2-1-14 the possible crosstalks

Among quorum sensing systems, we put focus on rhl system because of absence for in vitro crosstalk regarding Rhl system. Our candidate for the counterpart for Rhl system is lux system because Rhl system in Pseudomonas aeruginosa is known to be acticvated by Las system. Although combination between Las and Lux is also attractive candidate, Andrew.D et., al have shown crosstalk between them. (2018, Andrew.D). Therefore, we investigated whether there is a crosstalk between the lux and rhl QS systems in vitro. As RhlR-responsible promoters, we used a promoter set which showed improvement in S/N ratio in vivo (Tokyo tech iGEM 2014).

2-1-15  kinds of Prhl
2-1-15 kinds of Prhl

For the crosstalk experiments, we constructed two kinds of cell-free translation systems containing a substantial amount of R protein (luxR,rhlR). In those systems, we put signaling molecule (3OHSL-C4 or 3OHSL-C6) and target genes (plux/tet-GFP or prhl-RR-GFP), respectively. We measured the fluorescence of GFP expressed from plux/tet-GFP and prhl(RR)-GFP on a RT-qPCR machine. By measuring the fluorescence of FITC atthe same time, we were able to quantitatively compare the fluorescence values of GFP in different runs.

2-1-16 Result of crosstalk assay of the RhlR in cell-free system
2-1-16 Result of crosstalk assay of the RhlR in cell-free system

Results

Firstly, we did not fined crosstalk between RhlR activator and Plux promoter Fig 2-1-16. In contrast to strong activation in the combination among 3OC4HSL , RhlR, and Prhl(RR) promoter, the other combinations showed little fluorescence.

2-1-17 Result of crosstalk assay of the luxR in cell-free system
2-1-17 Result of crosstalk assay of the luxR in cell-free system

Unfortunately, however, we found Severe crosstalk from LuxR activator. As shown in two red bars, LuxR protein with OC6HSL activates not only Lux promoter but Rhl promoter. Althogh we tested another LuxR-responsible promoter, PRhlLR, it showd another cross talk from RhlR protein (三階層リンク). This project thus rejected this promoter variant.

 Fig 2-1-18 summarize our crosstalk analysis using cell-free system.
Fig 2-1-18 summarize our crosstalk analysis using cell-free system.

Even with the identified severe crosstalks between Rhl promoter and LuxR-3OC6HSL complex, reduction of maximum activity of Rhl promoter allowed us to prepare test tubes for Samurai and Zombie. Because of the same activation by samurai signaling molecule for the both repressors, Samurai state seems to be not stable. As describe above, nullclines for our system without crosstalk showed three intersections and thus implied bistability (Fig.2-1-19 left). By additional effect from the messy cross talk, shape of nullclines changes dramatically. Now, two null clines have only one intersection. This means loss of bistability and we cannot discriminate zombie and samurai state (Fig.2-1-19 center).

However, we could reduce CI repressor production from Rhl promoter by mutations at its DNA sequence. Even though such reduction had a risk of unitability of Zombie state, adequate modulation of maximum activity of Rhl promoter kept the stability. When we reduce the maximum expression by Rhl promoter to 1/10 compared with Lux promoter activity, phase space analysis show both stabilities of Zombie and Samurai state; the reduction did not change the number of intersections of nullclines from the case without crosstalk (Fig.2-1-19 right). Moreover, for this model we drew Fig.2-1-20 as a time course of mixed Zombie and Samurai. We found that even if we consider the crosstalk, the two powers can compete and be biased to either side depending on their ratios. In other words,

 Fig.2-1-20 War in the cell-free system containing crosstalk between R proteins and Promotes
Fig.2-1-20 War in the cell-free system containing crosstalk between R proteins and Promotes

[1] Maurizi MR. Proteases and protein degradation in Escherichia coli. Experientia. 1992 Feb 15 48(2):178-201
[2] Kigawa, T., Yabuki, T., Matsuda, N. et al. Preparation of Escherichia coli cell extract for highly productive cell-free protein expression. J Struct Func Genom 5, 63–68 (2004).

Entrepreneurship

SAMURAI System

Fig2-3-1 Overview of SAMURAI system
Fig2-3-1 Overview of SAMURAI system

The smells we have synthesized is not only for Zombie and Samurai in the story but for us in real-world. From initial stage of our project, our team has thought that synthetic biology can attract general public if we establish an interesting system.

Combined with our synthesis of limonene and sabinene, the incoherent feedforward loop (IFFL) system allowed us to design “Switching And Modulated Utilization of Refreshing Aroma Integration system”, SAMURAI system, where pulse expression of each of the two monoterpene synthases switches time-dependently.

IFFL

in order to change the scents with a time lag, we designed fragrance production regulated by IFFL(Incoherent Feed Forword loop), which can generate pulse activation of a gene. A typical IFFL consists of an activator protein (X), and a repressor protein (Y) ,and a taget gene (Z).

A target gene Z is directly activated by protein X. Aditionally, gene Z is indirectly repressed by protein X due to accumulation of repressor Y produced by X. Because this direct inhibition requres longer time than direct activtion, the time course of Z production shows a pulse

Fix 2-3-2 Explanation of incoherent feed forward loop
Fig 2-3-2 Explanation of incoherent feed forward loop
Fig2-3-3 Explanation of some point of the concentration pulse
Fig2-3-3 Explanation of some point of the concentration pulse

For time-dependent change of fragrances, we designed an expanded IFFL system containing two target genes (Z1,Z2) each of which has different peak time to the other. By adjusting the parameters for each activation and repression, we can provide a time difference in the generation of each pulse.

Fig2-3-4 Multi IFFL gene circuit
Fig2-3-4 Multi IFFL gene circuit
Fig 2-3-5 ideal behavier of multi IFFL
Fig 2-3-5 ideal behavier of multi IFFL

By using that system, we can change fragrance automatically as time passed

In order to check whether the designed multi-target IFFL can really generate two pulse with a time delay, we extensively searched parameter space for the system

Fix 2-3-6 Independent param search of Kx(left) and Ky(right)
Fig 2-3-6 Independent param search of Kx(left) and Ky(right)

We first investigated how the concentration of Z changes with time when the parameters Kx and Ky are changed independently in Formula xx.1.

Kx, Ky, are the binding coefficients of X and Y to the promoter of gene z. We varied those parameters because actual coefficient values can be controlled by changing the sequence of promoter.

As a result, we found that when the parameter Kx is high, the peak of the pulse become small, and when Ky is high, the sharpness of the peak is lost.

Logo
Fig 2-3-7 time heatmap

Next, we constructed a heatmap to check how much time (t) of the pulse peak is delayed as we change the parameters Kx and Ky simultaneously. Fig. 2-3-7 is the results, which show that the lighter the heatmap’s color is, the more the peak time delays. However, when the values of parameters become too large, pulse lost it’s sharpness as the the difference between the maximum value and the steady-state concentration decreases.

Logo
Fig 2-3-8 ratio heatmap
Logo
Fig 2-3-9 cross heatmap

Furthermore, in order to evaluate the balance of the pulse shape, we calculated the ratio of the peak concentration and the following steady-state concentration for each parameter and presented it in a heat map. This is called the ratio heatmap. (Fig 2-3-8). The darker the color, the smaller the ratio is and the pulse will be in an indistinct form as shown in Fig 2-3-7, Kx=1000.

Fig 2-3-8 shows that the ratio becomes smaller when the Kx value is too high. We evaluated each parameter based on the peak duration in Fig 2-3-7 and the balance of the pulse shape in Fig 2-3-8, respectively. Based on these two indicators, we multiplied the cross heatmap by the value of the time heatmap and the ratio heatmap for each parameter in order to find out which parameter has both time difference and balance, comprehensively (Fig 2-3-9). To make this cross heatmap, we subtracted a certain constant number of time values from the time heatmap because we have to make the time delay above a certain level. If the time heatmap is negative, it is treated as 0.

This is called the time processed heatmap. Then, when multiplying the values of the time processed heatmap and the ratio heatmap, the value of the time processed heatmap was raised to the power of 1.2 and the value of the ratio heatmap was raised to the power of 10 to adjust for the contribution of the respective heatmap values. We selected the optimal parameter for the multi IFFL in the resulting match heatmap as the one with the highest score. Two pulses were generated using the optimal parameter in practice(Fig 2-3-10). Yellow indicates the temporal variation of the concentration of Z1 and green indicates that of Z2. As a result, We have succeeded in creating pulses with a time difference. This demonstrates the feasibility of our SAMURAI system.

Although we sure that you have enjoyed our SAMURAI, however, we also know that real-world implementation requires another effort: business model creation. This is why we've been working on the business model from the early time of the project.

Logo
Fig 2-3-10 results


Entrepreneurship

Logo
2-3-11 Overview of business model innovation

As well as planning for scientific research, business model creation should learn from previous success. From Mr. Watanabe’s point of view, a business has a higher probability of success when two of a four points are changed from the existing business success. The kind of change is called business model innovation.

Logo
Fig 2-3-12 Interview of business model

With with the knowledge we acquired from lectures on business model planning, we created some business models based on our cell-free based iGEM project. Although we created detailed plan which can take advantage of our success in monoterpene syntheses , another interview to experts during the above planning gave us more attractive and valuable business model where a cell-free system is used for substance detection.

Cell-free monoterpene synthesis

Background

How can we observers distinguish between zombies and samurai? In our project we propose one unique identification method, though there is also a method using fluorescent protein. We propose a method of distinguishing the scents of zombies or samurai. For this project, we selected monoterpenes, especially limonene and sabinene, as the scented substances produced by zombies and samurai.

Experiment

Parts_and_metabolic_pathways_in_this_experiment_and_a_schematic_diagram_of_the_experiment
Fig. 2-2-1 Parts and metabolic pathways in this experiment and a schematic diagram of the experimen

In this cell-free monoterpene synthesis, we mixed two E. coli extracts each of which has either first 7 or last 2 enzymes of a pathway from Ac-CoA, which is a major intermediate of cell central metabolism. Through mevalonate pathway, the former extract one (derived from E. coli into which pBbA5c-MevT-MBI has been introduced) can provide IPP and DMAPP, which can also be used as intermediates for other important biosynthesis.

Here we indeed supplemented only glucose and acetate as carbon sources. We obtained expression system for those seven genes from addgene and have converted this into Biobrick RFC 1000 format by synonymous replacement ( BBa_K3580103). In order to take advantage of an engineering principle of synthetic biology we provided two biobrick parts ( BBa_K3580101, BBa_K3580102 ) for the source for the latter extract. BBa_K3580101 has GPP synthase (GPPS) and limonene synthase. Although GPP synthase is shared with BBa_K3580101 , BBa_K3580102 has sabinene synthase, which has one point mutation in limonene synthase (Srividya Narayanan et al 2015) and a new coding sequence for Parts registry of iGEM (See here for more details on this experiments).

Results

Similar to a very recent study from Jwett lab (Dudly et al 2019) who mixed 7 extracts, we synthesized limonene using only two extracts. Although fine-tuning can be possible when a larger number of extracts is prepared, we are sure that entry projects in iGEM should be simple but has engineering principle in order to expand iGEM sucess to an educational tool. This is why we selected the division into two extracts.

Fig. 2-2-2 GC/MS analysis results of limonene synthesis system
Fig. 2-2-2 GC/MS analysis results of limonene synthesis system

We confirmed limonene synthesis GC/MS analysis with SIM. In this SIM analysis, ions with four m/z values characteristic in limonene (68 and 93) and sabinene (77 and 91, 93) were analyzed. By using authentic limonene standard, we confirmed a retention time for GC and the characteristic limonene SIM signal at the specific m/z values.

At the same retention time with the standard, limonene-specific m/z value (68, 93) ions were detected in the selected ions (The upper right figure of Figure. 2-2-2). We also draw a GC chart by summation of the signals from the selected ions (The lower left figure of Figure. 2-2-2). By comparison with a negative control experiment which we omitted the extract containing GPP synthase and limonene synthase, we found clear peak from our limonene synthesis.

Fig. 2-2-3 GC/MS analysis results of sabinene synthesis system
Fig. 2-2-3 GC/MS analysis results of sabinene synthesis system

By similar GC/MS analysis with SIM, we confirmed the world's first report of sabinene synthesis using cell extracts. Although elaborative reconstituted system by mixing of purified enzymes has been reported (Srividya Narayanan et al 2015), (Korman Tyler P et al 2017). Extract-base system with not only simple but engineering principle must provide much contribution to iGEM and DIY biology which as large number of players. Furthermore, extract systems do not require expensive coenzymes (Dudly et al 2019). As well as limonene detection in Fig 2-2-2, we firstly confirmed sabinene retention time and SIM signal at the m/z values (77, 91, and 93).

From our cell-free production, we then detected sabinene SIM signal and chromatogram peak which cannot find from the negative control. Also, from the system containing sabinene synthase and GPP synthase, a peak with the same retention time as the standard limonene product and ions with m/z values characteristic in limonene at that retention time were detected. This was consistent with the fact that the paper that referred to the sabinene synthase used in this experiment also reported the synthesis of limonene as a by-product (Srividya Narayanan. et al 2015).

The part BBa_K3580102 used in our sabinene production is a new part for Parts Registry. One amino acid mutation was introduced into limonene synthase of Mentha spicata to create sabinene synthase, by using information reported in a 2015 paper by Srividya and Narayanan et al.

We regard our BBa_K3580101 as an improved version from pre-existing BBa_K3580201, which has only limonene synthase. Our improved part can use the set of IPP DMAPP both of which are provided not only mevalonate pathway used here but the other famous non-mevalonate pathway. Thus, the set is produced in so many engineered pathways, including the addgene plasmid used here, as an important intermediate set for bioproduction. From this view point, we are sure that the simple combination of two enzymes in this part satisfied easy use modularity in synthetic biology than pre-existing BBa_K3052001.

Fig. 2-2-4 Yields of monoterpene per reaction solution in this experiment
Fig. 2-2-4 Yields of monoterpene per reaction solution in this experiment

Finally, each of monoterpenes were quantified based on peaks of substance having 93 m/z and each monoterpene standard curves. Taking advantage of the modularity of the combination of extracts, we confirmed whether the yield of monoterpenes could be changed by changing the mixing ratio of the extract containing the enzyme of the mevalonate pathway and the extract containing GPP synthase and monoterpene synthase. As a result, changes in the yield of monoterpenes due to the mixing ratio of the extracts were observed (Figure. 2-2-4). The best yield of limonene per reaction solution with limonene synthase contained system was 0.73 µM, and the best yield of sabinene per reaction solution with sabinene synthase contained system was 4.6 µM.

Reference

(1) Alonso-Gutierrez, J. et al. (2013). Metabolic engineering of Escherichia coli for limonene and perillyl alcohol production. Metabolic engineering, 19, 33-41.
(2) Dudley, Q. M. et al. (2019). Cell-free biosynthesis of limonene using enzyme-enriched Escherichia coli lysates. Synthetic Biology, 4(1), ysz003.
(3) Srividya, N. et al. (2015). Functional analysis of (4S)-limonene synthase mutants reveals determinants of catalytic outcome in a model monoterpene synthase. Proceedings of the National Academy of Sciences, 112(11), 3332-3337.
(4) Korman, T. P. et al. (2017). A synthetic biochemistry platform for cell free production of monoterpenes from glucose. Nature communications, 8(1), 1-8.
(5) Mass Spectrometry Data Center, William E. Wallace, "Mass Spectra" in NIST Chemistry WebBook, NIST Standard Reference Database Number 69, Eds. P.J. Linstrom and W.G. Mallard, National Institute of Standards and Technology, Gaithersburg MD, 20899, https://doi.org/10.18434/T4D303, (retrieved October 22, 2020).

Racemase ~War In The Mirror~

Introduction

In the Zombie vs. Samurai storyline, we envisioned a scenario in which one side deprives their opponent of the food they need. A battle between two enantiomers, that is to say, one side has D amino acid and the other has L amino acid derived food as their food source, and both sides converts their opponent’s amino acid and deprives the of their food (Fig.2-4-1). Although there is not yet a translation and transcription system that efficiently incorporates D amino acids, it is known that the originally natural enzyme chemically synthesized with D amino acids would be active.

Fig.2-4-1 Images of War In The Mirror in the Zombie vs Samurai storyline
Fig.2-4-1 Images of War In The Mirror in the Zombie vs Samurai storyline

Although there is not yet a translation and transcription system that efficiently encorporates D amino acids, it is known that the natural enzyme chemically synthesized with D amino acids would be active.A protein made up of only d-amino acids act on substrates with opposite chirality, as demonstrated using the HIV protease [Milton et.al 1992].

Besides, practical usage of D-amino acids is advancing in the field of Biotechnology today. Peptides composed of D-amino acids are not susceptible to degradation by protease in L-body organisms, thus are highly stable [Tugyi et.al 2005][Garton et.al 2018]. Consequently, researchers are seeking a way to make the target of D-amino acid Peptide Aptamers(A type of antibody) L-body organisms’ protein [Oberthür et.al 2015][ Majier et.al 2016][ Schumacher et.al 1996][ Funke et.al 2009]. Currently, target proteins are synthesized with D-bodies, and L-body peptide aptamers that bonds to the protein are artificially evolutionized. When protein is expressed with the obtained amino acid sequence with D-amino acids, the binding site would consist of L-bodies [Oberthür et.al 2015][ Majier et.al 2016][ Schumacher et.al 1996][ Funke et.al 2009].

Such storyline can be applied to our real experimental system by adopting racemase, an enzyme which catalyzes interconversion of D amino acid and L amino acid. Racemase catalyzes interconversion between L-amino acid and D-amino acid. There are many kinds of racemases and each corresponds to respective amino acid. This year, we focused especially on alanine racemase. Alanine racemase (AR) is a fold type III racemase enzyme catalyze the conversion of L-alanine to D-alanine[Walish 1989]. This enzyme provides D-alanine consuming L-alanine and using pyridoxal 5- phosphate (PLP) as a cofactor. AR is unique to prokaryotes and plays an important role in a biosynthesis[Walish 1989] [Azam et.al 2016]. Since D-alanine is used as a necessary component of the peptidoglycan layer of bacterial cell walls, lack of alanine racemase can result in the inhibition of growth of prokaryotes [Azam et.al 2016]. As it is generally absent in higher eukaryotes such as human but is ubiquitous throughout bacteria. Therefore, alanine racemase is recognized as an attractive target for antibacterial drug development [Azam et.al 2016].

Summary & Experiment

In this experiment, alanine racemase (AR) gene carried on our new part (BBa_K3580200 ) on a T5 promoter-controlled pCA24N vector was transformed into the BL21(DE3) star strain. After expression induction by IPTG, the BL21 cells were harvested by centrifugation, then sonicated and His-Tag purification was performed (Fig.2-4-2).

Fig.2-4-2 Purification of alanine racemase expressed in E. Coli
Fig.2-4-2 Purification of alanine racemase expressed in E. Coli

To assess the activity of purified racemase, fluorescence values of GFPS1 synthesized by a combination of GFPS1 plasmid, L-alanine, D-alanine, and purified AR were measured in the Pure System, a reconstituted CFPS that contains no enzymes other than the translation system [Shimizu et.al 2001] (Fig.2-4-3).

Fig.2-4-3 Cell-free protein synthesis in the PURE system
Fig.2-4-3 Cell-free protein synthesis in the PURE system

The natural translation system is difficult to take up D-amino acids.[Kuncha 2019].19 amino acids other than alanine are L-bodies, so if only alanine is D-body, even if the sequence is correct, it has no activity as an enzyme (Fig.2-4-4).

Fig.2-4-4 Cell-free protein synthesis under each alanine conditions
Fig.2-4-4 Cell-free protein synthesis under each alanine conditions

Since GFP is composed of L-amino acids, when GFPS1 takes up L-alanine, the system would give off fluorescent effect, and when D-alanine and AR are added, the racemized L-alanine can be used for translation and the fluorescence value was expected to be restored(Fig.2-4-5). The evaluation of the effect of amino acids on the synthesis of a reporter protein is the most appropriate methodology for this scenario as a "food source". As PURE frex does not contain any extra metabolic enzymes other than those involved in transcription and translation, we would be able to perform the assays without having to control the metabolic system, which could plague the assays with extract based CFPS.

Fig.2-4-5 Cell-free protein synthesis using racemization of alanine by racemase
Fig.2-4-5 Cell-free protein synthesis using racemization of alanine by racemase

We decided to assay AR by adopting PURE frex donated by GeneFrontier, Inc. After identifying L-alanine contamination problems of commercial D-alanine, the assay with AR was confirmed by the combination of each alanine and racemase as well as the assay with extracts. The results of this experiment are as follows (Fig.2-4-6).

Fig.2-4-6 Effect of L-alanine substrate repletion for translation by racemase
Fig.2-4-6 Effect of L-alanine substrate repletion for translation by racemase

The results demonstrated that GFP could not be synthesized by D-alanine alone, but L-alanine produced by racemase-mediated racemization of D-alanine and L-alanine could be used for translation, creating a situation in which fluorescence was restored by the synthesized GFP. In other words, we succeeded in creating a situation in which the opponent’s food was converted into food for themselves in the Zombie vs. Samurai scenario.

We also used a cell-free transcription-translation system using cell extracts prior to our experiments with the pure system, but we were unable to obtain clear results due to the production of alanine by the metabolic enzymes in the extracts.

Reference

(1) [Milton et.al 1992] Milton, R., Milton, S., & Kent, S. (1992). Total chemical synthesis of a D-enzyme: The enantiomers of HIV-1 protease show reciprocal chiral substrate specificity [corrected]. Science, 256(5062), 1445-1448. doi: 10.1126/science.1604320
(2) [Tugyi et.al 2005] Tugyi, R., Uray, K., Iván, D., Fellinger, E., Perkins, A., & Hudecz, F. (2005). Partial d-amino acid substitution: Improved enzymatic stability and preserved ab recognition of a MUC2 epitope peptide. Proceedings of the National Academy of Sciences, 102(2), 413-418. doi: 10.1073/pnas.0407677102
(3) [Garton et.al 2018] Garton, M., Nim, S., Stone, T. A., Wang, K. E., Deber, C. M., & Kim, P. M. (2018). Method to generate highly stable D-amino acid analogs of bioactive helical peptides using a mirror image of the entire PDB. Proceedings of the National Academy of Sciences, 115(7), 1505-1510. doi: 10.1073/pnas.1711837115
(4) [Oberthür et.al 2015] Oberthür, D., Achenbach, J., Gabdulkhakov, A., Buchner, K., Maasch, C., Falke, S., Rehders, D., Klussmann, S., Betzel, C. (2015). Crystal structure of a mirror-image L-RNA aptamer (spiegelmer) in complex with the natural L-protein target CCL2. Nature Communications, 6(1), 6923. doi: 10.1038/ncomms7923
(5) [Majier et.al 2016] Maier, K., & Levy, M. (2016). From selection hits to clinical leads: Progress in aptamer discovery. Molecular Therapy.Methods & Clinical Development, 5, 16014. doi: 10.1038/mtm.2016.14
(6) [Schumacher et.al 1996] Schumacher, T. N. M., Mayr, L. M., Minor, D. L., Milhollen, M. A., Burgess, M. W., & Kim, P. S. (1996). Identification of d-peptide ligands through mirror-image phage display. Science, 271(5257), 1854-1857. doi: 10.1126/science.271.5257.1854
(7) [Funke et.al 2009] Funke, S. A., & Willbold, D. (2009). Mirror image phage display—a method to generate d-peptideligands for use in diagnostic or therapeutical applications. Mol.BioSyst., 5(8), 783-786. doi: 10.1039/B904138A"
(8) [Walish 1989] Walsh, C. T. (1989). Enzymes in the D-alanine branch of bacterial cell wall peptidoglycan assembly. The Journal of Biological Chemistry, 264(5), 2393-2396.
(8) [Azam et.al 2016] Azam, M. A., & Jayaram, U. (2016). Inhibitors of alanine racemase enzyme: A review. Journal of Enzyme Inhibition and Medicinal Chemistry, 31(4), 517-526. doi: 10.3109/14756366.2015.1050010
(9) [Shimizu et.al 2001] Shimizu, Y., Inoue, A., Tomari, Y., Suzuki, T., Yokogawa, T., Nishikawa, K., & Ueda, T. (2001). Cell-free translation reconstituted with purified components. Nature Biotechnology, 19(8), 751-755. doi: 10.1038/90802
(10) [Kuncha et.al 2019] Kuncha, S. K., Kruparani, S. P., & Sankaranarayanan, R. (2019). Chiral checkpoints during protein biosynthesis. The Journal of Biological Chemistry, 294(45), 16535-16548. doi: 10.1074/jbc.REV119.008166

Education

Overview

At initial stage of our iGEM activity, we realized that not only those with little scientific knowledge but also students who aspired to be scientists didn’t recognize synthetic biology. As iGEMers, who mainly work with synthetic biology, we felt the need to spread the knowledge of its potential to produce world-changing results in both engineering and physiological sciences. This year, we’re attracting attention with an interesting story whose subject is gene circuit. We tried to spread toggle switch, which was core of the circuit and starting point of synthetic biology.

Logo
Fig 2-6-1 EducationApplication student

We decided to develop a smartphone application (Fig. 2-6-1) for better understanding of mathematical modeling of synthetic biology for the synthetic biology beginners. Mathematical modeling is fundamental in synthetic biology, for engineering principle requires such modeling. We fortunately had a chance to lecture Waseda University High School students. The lecture was held through a widely used video communication service, zoom.

fig 2.7 toggle switch simulator
fig 2-6-2 toggle switch simulator

The application shows how genetic circuit works, especially about toggle switch. Rates of change in repressor concentrations move according to parameter values(⑤). Resulting in visual information such as graph or color transitions(③,④,⑥). Each student can move the parameter values(⑤), watching the feedback, and learn how each parameter affects the result. By operating the app hands-on, we received positive comments from the students. Some said they felt fresh about the way they recognized organisms by using mathematical model, and others said they found utility in this method. Afterwards, this smartphone application was incorporated in a university class because of its usefulness.

2020 — iGEM Waseda