Team:UT-Tokyo/Project

INTRODUCTION

Project Overview

Mechanisms for development of living things have been investigated for a long time. In 1952, Alan Turing made a key contribution. He showed that the interaction of two substances with different diffusion rates can generate spatial concentration patterns, which are called Turing Patterns, without any preformed pattern[1]. This idea has been applied to explain hair follicle distribution of mice[2]<, stripe formation on zebrafish[3], and the determination of pigment pattern within avian feather follicles[4].

Here, We tried to reconstruct Turing Pattern by E. coli in a way of synthetic biology to understand the mechanism more. E. coli has cell-cell communication system called quorum sensing. Utilizing this system, we aimed to generate a periodical colony pattren. This project leads to further understanding of Turing Pattern and development of living things.


Turing Mechanism

Turing Pattern is generated by the interaction of two substances with different diffusion rates. These substances are called activator and inhibitor. Activator promotes its own production and the production of inhibitor, and inhibitor inhibits the production of activator(Fig.1). Inhibitor diffuses faster than activator.

Fig.1 Turing mechanism The interaction of activator and inhibitor generates patterns.

Turing focused on the difference of the diffusion rates, but later Meinhardt and Gierer explained the mechanism from the point of local self-activation and lateral inhibition[5]. This explanation is intuitively easy to understand, so we will explain Turing pattern along this explanation.

Fig.2 The generation of a periodical pattern
The initial condition is (a), and as time passes, a periodical concentration pattern is generated((b)~(f)).

At first, there is no preformed pattern in the field and the concentrations of two substances are homogeneous through the field(Fig2. (a)). However there can be a spot with a little bit high concentration of activator because of perturbation. At such a point, the concentration of activator becomes higher and higher, as well as that of inhibitor, by its function(Fig2. (b)~(c)). This corresponds to local self-activation. Since inihibitor diffuses faster than activator, the concentration of inhibitor becomes relatively high around that point and that of activator becomes low because of the inhibitory effect of inhibitor there(Fig2. (d)~(e)). This corresponds to lateral inhibition. Therefore, a spot of activator is generated. This reaction occurs at different points in the field, and the distance between each spot is regulated by the interaction of the two substances(Fig2. (f)). In this way, a periodical pattern is generated in the field.


Advantage Of Synthetic Biological Approach

Pattens of living things have been investigated for a long time, but it was not easy to prove directly if these patterns are produced by the reaction-diffusion systems or another mechanism. Living systems are so complex that most research was exclusively theoretical. Biologists still face a big problem: identification of proper molecules acting as activator and inhibitor.

We therefore reconstructed a Turing system using two advantages of synthetic biology; controllability and biological directness. We can change the diffusion rate of E. coli and the strength of inhibitory effects of inhibitor by inducing synthetic circuit, which can be a great advantage of the experimental system. Chemical system has a similar advantage, but it is far from living systems. Our system uses cells themselves for pattern formation, so it may be directly applied to developmental studies.

STRATEGY

We designed three different strategies to realize pattern formation.

Strategy 1

The first one is based on Turing’s classic model. (Figure1)

Fig.1 Model 1
Imitation of Turing's classic model performed by E. coli and AHL

Here, we consider the “reaction” between E. coli and AHL. AHL is a type of intercellular signaling molecule that can promote transcription from specific promoter when its concentration in a cell gets higher than certain threshold. This system is called ”Quorum Sensing”. And here, AHL is set to induce protein which inhibits the multiplication of E. coli.

E. coli produces(activates) AHL, and AHL inhibits the increase of E. coli indirectly. E. coli also multiples itself.

Diffusion rate of AHL is fast because its size is quite small. On the other hand, the rate of E. coli should be low considering Turing’s model, and to achieve this, we modified E. coli by knocking out certain gene concerned with its motility.

So, the control loop 1 (Figure1-1) functions as lateral inhibition, and the control loop 2 (Figure1-2) functions as local activation.

Fig.1-1 Control loop 1
This loop is long-range negative feedback, and it drives lateral inhibition when the concentration of E. coli gets higher.

Fig.1-2 Control loop 2
This loop is short-range positive feedback, and it drives local activation when the concentration of E. coli gets higher.

But this model has one potential defect: the order of the growth rate of E. coli is expected to be expressed as linear of the concentration of E. coli, and primary order multiplication may be too small as autocatalysis(local activation) in Turing’s model.

To solve this problem, we introduces another type of E. coli as co-activator, and that is the second strategy. (Figure2)


Strategy 2

Fig.2 Model 2
Improved model from Model 1

E. coli(activator) activates AHL, and AHL inhibits E. coli, as in the model in Figure 1. This loop functions as lateral inhibition. (Figure2-1)

Fig.2-1 Long-range negative feedback loop
Same loop as in Figure 1-1

Two types of E. coli (activator and co-activator) inhibits each other(Figure2-2), and this interaction is achieved by the function of the protein called Colicin. (The detail is shown in the chapter of System.)

Fig.2-2 Short-range positive feedback loop
Mutual inhibition forms positive feedback. And the slow diffusion rate of the agent substance(Colicin) makes the range of this loop short.

In this loop, the increase of $u$ causes the decrease of $v$, and the decrease of $v$ causes the increase of $u$ in turn. And the diffusion rate of Colicin is small for its big size[6]. Therefore, we can regard it as local activation.


Strategy 3

The third strategy is inspired by the pattern formation mechanism of zebrafish[3]. (Figure 3)

Fig.3 Model 3
Two types of E. coli play the role of activator and inhibitor.

In Figure 3, two types of E. coli(activator and inhibitor) reacts each other in two manners according to the distance between them.

When E. coli(activator) and E. coli(inhibitor) are close, they inhibit each other in the same way of strategy 2. (Figure 3-1) This control loop functions as local activation.

Fig.3-1 Short-range positive feedback loop
Same loop as in Figure 2-2

In addition, E. coli(activator) activates E. coli(inhibitor) through AHL in long distance(detailed scheme is explained in the chapter of System). The control loop composed of this activation and the short-range inhibition of E. coli(activator) by activated E. coli(inhibitor) functions as lateral inhibition. (Figure3-2)

Fig.3-2 Long-range positive feedback loop
"Activator" E. coli promotes the multiplication of "Inhibitor" E. coli, and "Inhibitor" inhibits the multiplication of "Activator" near of "Inhibitor"

SYSTEM

In this chapter, we give concrete constructions for the concepts explained at Strategy. First, we explain basic mechanisms which are necessary to understand our construction.

1. Motility control

1.1 CheZ

Fig.1-1 Swimming phase and tumbling phase
The motility of the cell is regulated by CheZ.
(a)Swimming phase: CheZ is expressed and inhibits CheY.
(b)Tumbling phase: CheY binds to flagellar binding proteins.

The movement of E. coli is roughly divided in two phases, swimming(a) and tumbling(b). In order to convert these phases, E. coli use flagellar binding protein called CheY which is controlled by CheZ. Under the expression of CheZ, CheY is dephosphorylated by CheZ and inactivated to bind to the flagellar mortar proteins. Consequently, the rotation of flagellar is changed and cell initiates swimming phase.[7]
In our project, we use cheZ knock out strain (JW18XX), and transfer cheZ under inducible promoter to control cell motility.


2. Cell – Cell interaction

2.1 Quorum Sensing

Fig.2-1 Quorum sensing
Blue hexagon shows AHL molecule. When population density gets higher than certain threshold, the transcription from specific promoter (such as Plux) is enhanced.

AHLs are signal molecules involved in bacterial quorum sensing (Fig 2.1). AHL can easily permeate cell membranes and can regulate the transcription of target cell. For example, in Lux system, AHL binds to LuxR dimer and that complex enhances transcription from PRlux promoter.


2.2 Colicin

Fig.2-2 Mechanism of Colicin release
This figure may a littel bit differ from actual mechanism.

Colicin E3 (ColE3) is a ribonuclease. Colicin E3 Immunity protein (ColI) binds to ColE3 and neutralize it. Colicin Lysis protein (ColL) allow ColE3 to pass thorugh cell membrane. The mechanism of colicin release has not been elucidated.[8][9]
Colicins are a cytotoxins which are released to environment and kill other related strains. Release of colicin involves one protein; Colicin Lysis Protein (ColL). Colicin lysis protein allows colicins to be released. The mechanism how the colicin lysis protein allows colicin release has not been fully elucidated, but it is sure that this protein raise membrane permiability and cause quasilysis.[8] After Colicin released, they diffuse through the medium and bind to the receptor on the target cell membrane. Then, they are imported to the cytoplasm or cytoplasmic membrane of target cell by Tol-system or Ton-system. (If you want to know the mechanism of the Colicin import, see [8])
Colicin producing cells also express Colicin Immunity Protein (ColI) in order to protect themselves from cytotoxity of colicin.
Colicins have variety cytotoxity such as DNase activity, RNase activity or Pore forming across inner membrane.
In our project, we select colicin E3 because it has low risk of safety problems. Colicin E3 specifically digests 16S rRNA which is only in bacteria.


3. Construction

Strategy 1

Fig.3-1 Construction for Strategy 1
Positive feed back loop, which is necessary for local activation, is played by multiplication of E. coli.

Gene modified E. coli is an activator, and AHL is an inhibitor.
Local activation is played by self-reproducing of E. coli.
Lateral inhibition is played by AHL. AHL activates the expression of Barnase (RNase Ba) which causes E. coli to die. Barnase is the RNase from Bacillus amyloliquefaciens.
Since gene modified E. coli is cheZ knock out strain, the difference between the diffusion rate of E. coli and that of AHL is enough to create Turing Pattern.

Strategy 2

Fig.3-2 Construction for Strategy 2
Inhibit means growth competiton. When colicin sensitive cells (in this figure, activator) grow faster than colicin producing cells(co-activator) and initial concentration of activator is higher than a certain threshold, colicin sensitive cells continue to increase and colicin producing cells die.[10]

In order to strengthen positive feedback of local activation, another type of E. coli (co-activator) is added to Strategy 1. Lateral inhibition is played by AHL.
Co-activator produces colicin, which represses the growth of activator. Note that the growth rate of co-activator is modulated to be lower than that of activator. As a result, they compete and repress each other.[10] This feedback loop acts local activation and helps pattern formation.

Strategy 3

Fig.3-3 Construction for Strategy 3

Two types of E. coli play activator and inhibitor.
Local activation is same system as Strategy 2. Activator and inhibitor repress each other in short range.
Lateral inhibition is played by AHL and inhibitor cell. AHL is produced by activator cell and enhances the expression of Barstar, Barnase immunity protein, in inhibitor cell. Thus, the Barnase is inactivated and the growth rate of inhibitor cell is recovered. As a result, inhibitor cell represses activator cell by growth competition.


4. Assay

Fig.4-1 Pattern formation assay

We observed the pattern on semi-agarose gel. Concentration of agarose is low(0.15%), so that E. coli

RESULT

今後一か月の実験の中で出すべき結果(想定)を書く

本来示すべきこと(理想)

  • SSA(semi solid agar)上での大腸菌とAHLの拡散速度測定

  • 大腸菌とAHLの初期濃度によってSSA上で大腸菌濃度が双安定状態をとること(拡散ありで平衡点が不安定)

  • 大腸菌とAHLを試験管に入れて大腸菌濃度が一定値に収束すること(拡散なしで平衡点が安定)

実際に示せそうなこと(想定)

大腸菌とAHLの拡散係数

大腸菌:タイムラプスカメラでコロニーの拡大を撮影。コロニーの境界が動く速さを測定
AHL:pLux-gfp、pconst-luxRを持った大腸菌を一様に培養してあるプレートの中央にpconst-luxIを持った大腸菌を植菌し、GFPが発現している領域が拡大する様子を撮影。 適当に決めた蛍光強度の線が動く速さを測定(どの値を取るのが適切か要検討)
文献値を使ったシミュレーションと比較してfittingにより妥当な拡散係数を出す。
文献値を使ったシミュレーションと一致すれば、その値をそのまま使う。一致しなければ、値を調整して実験結果をよく説明する値を採用する。

測定結果から拡散係数を直接導出する方法もないわけではない。
枯草菌のバクテリアコロニーを研究した論文。
コロニー内での個々のバクテリアの動きをビデオイメージによって時間的に追跡した。
バクテリア細胞が、ある点から出発して時間tの間に移動した直線距離Rを多数回観察し、平均値$&ltR^2>$を測定。
バクテリアの動きをブラウン運動だと仮定してアインシュタインの関係式
$&ltR^2> = 2D*t $ (Dは拡散係数)
に計測値を代入して実効的な拡散係数を導出している。
Experimental Investigation on the Validity of Population Dynamics Approach to Bacterial Colony Formation
wikipediaにアインシュタインの関係式の説明がある
Einstein Relation
ただし、拡散係数が何に依存するかを考察し、適切な近似をすることは必要。
例:濃度依存、温度依存等、すべて考慮すると導出できなくなるので近似をする。


collicinで大腸菌が死ぬこと

試験管でplac-collicinを持った大腸菌を培養しIPTG誘導かけて死ぬかをOD測定で見る。


Luxシステムの機能

pLux-gfp、pconst-luxRを持った大腸菌を試験管で培養し、AHLを濃度振って誘導し、蛍光強度を見る。
luxIを含んだ全体のシステムの確認はAHLの拡散係数測定で同時に行うことにする。


Strategy 2の平衡点安定性アッセイ

strategy 2の二種類の大腸菌をSSA上で培養、培養容器の大きさを変えて、拡散の効果が出てくる大きさを見つける。
容器が小さすぎると一瞬でAHLが拡散して大腸菌の成長が容器全体で抑制された結果、大腸菌密度は容器内で一様になるはず。
容器が大きければ、AHLの拡散によって容器中央はAHL濃度が低く、外縁はAHL濃度が高くなり、容器中央に大腸菌のスポットができるはず。
拡散係数測定のアッセイで出した拡散係数を使ったモデリングと結果を比較する(一致していてほしい)。
一点のスポットのみのチューリングパターン。スポット間の相互作用によるスポット間の位置調整(周期的パターンにつながる)は見られない。


大腸菌-AHLをSSAに播いてみた結果


strategy 3は未定。barstarアッセイ(plac-barstar、pconst-barnaseでIPTG誘導OD測定?)


The experiments results we are planning to show are...

The diffusion rates of E.coli and AHL

E.coli
Take pictures of a colony at regular intervals using a time-lapse camera and measure the expansion speed of the colony.

AHL
Culture E.coli with pLux-gfp and pconst-luxR circuit homogeneously on semi solid agar and inoculate another E.coli with pconst-luxI circuit on the center of a plate. Take pictures of the colony under exciting light with a time-lapse camera and measure the expansion speed of the fluorescent area, where E.coli receives AHL and expresses GFP.

Compare the results of the experiments and the computer simulation, and find the diffusion rates which fits the experiments results.


The cytotoxicity of collicin

Culture E.coli with plac-collicin circuit in a tube, induce IPTG and compare OD before the induction and after the induction.


The function of lux system

Culture E.coli with pLux-gfp and pconst-luxR circuit, induce AHL and measure the fluorescence.
Confirmation of the function of the whole lux system(including luxI) is included in the diffusion rates assay of AHL.


Strategy 2: Equibrium point stability assay

Culture two types of E.coli in strategy 2 in a plate. Change the size of a plate and find the size at which the diffusion become effectual. (Check)

APPLICATION

REFERENCE

[1]Turing, A. M. (1952). The chemical basis of morphogenesis. Philosophical Transactions of the Royal Society of London B: Biological Sciences, 237(641), 37-72.

[2]Sick, S., Reinker, S., Timmer, J., & Schlake, T. (2006). WNT and DKK determine hair follicle spacing through a reaction-diffusion mechanism. Science, 314(5804), 1447-1450.

[3]Nakamasu, A., Takahashi, G., Kanbe, A., & Kondo, S. (2009). Interactions between zebrafish pigment cells responsible for the generation of Turing patterns. Proceedings of the National Academy of Sciences, 106(21), 8429-8434.

[4]Prum, R. O., & Williamson, S. (2002). Reaction–diffusion models of within-feather pigmentation patterning. Proceedings of the Royal Society of London B: Biological Sciences, 269(1493), 781-792.

[5]Meinhardt, H., & Gierer, A. (2000). Pattern formation by local self-activation and lateral inhibition. Bioessays, 22(8), 753-760.

[6]Cascales, E. et al. (2007) Colicin Biology. Microbiology and Molecular Biology Reviews, 71(1), 158-229.

[7]Parkinson, J. S. (2003). Bacterial chemotaxis: a new player in response regulator dephosphorylation. Journal of bacteriology, 185(5), 1492-1494.

[8]Cascales, E., Buchanan, S. K., Duché, D., Kleanthous, C., Lloubes, R., Postle, K., ... & Cavard, D. (2007). Colicin biology. Microbiology and Molecular Biology Reviews, 71(1), 158-229.

[9]Lloubes, R., Bernadac, A., Houot, L., & Pommier, S. (2013). Non classical secretion systems. Research in microbiology, 164(6), 655-663.

[10]Chao, L., & Levin, B. R. (1981). Structured habitats and the evolution of anticompetitor toxins in bacteria. Proceedings of the National Academy of Sciences, 78(10), 6324-6328.