Bacteria expressing genes that code for color will form interesting patterns when spotted onto agar plates in appropriate locations and allowed to spread out and intermingle. But what if we could increase the complexity of the patterns formed by implementing intercellular signaling between different groups of bacteria?

We decided to use parts from the natural quorum-sensing mechanisms of various bacteria to implement our intracellular communication system. If it works, we can for example cause two colonies of bacteria to change color or stop moving as they approach each other, hopefully resulting in interesting patterns.

A brief overview of quorum-sensing: Bacteria, such as V. fisheri, coordinate their gene expression through a system in which each bacterial cell produces a limited amount of signaling molecules, called AHL, which diffuse through the medium and reach other bacteria in its vicinity. AHL molecules bind to receptor proteins which in turn bind to specific promoters that then up-regulate downstream transcription activity. When the bacteria reach a certain density, the amount of AHL in the environment (and thus in the cells) will be sufficient to trigger this increase in promoter activity, and the genes downstream of the promoter will be ‘switched’ on.

[include picture of quorum-sensing system]

Currently we are working on 2 distinct groups of parts: 'Senders' and 'Receivers'. 'Senders' code for enzymes that produce AHL signal molecules, which diffuse out of the cell, through the culture medium and into the receiving cell, where a receptor proteins encoded by the 'Receivers' bind the signals, forming a complex which in turn can bind to and up-regulate transcription from their specific promoter.

We assembled the following devices using parts from the iGEM 2009 Spring Distribution:

Senders:

Lux Sender - produces (Lux-signaling system AHL), includes a double terminator for easy insertion in front or behind any other device/part [include pic of circuit]
Las Sender - same as above, but produces (Las-signaling system AHL)
Rhl Sender
Cin Sender

Receivers:

Lux Receiver - receives signal from Lux Sender transmitted in the form of (Lux-signaling system AHL), which then activates/upregulates transcription downstream of this device [include pic of circuit]
Las Receiver
Rhl Receiver
Cin Receiver

Test Constructs:

Lux Receiver with GFP attached downstream ("X1") - a GFP coding device has been attached downstream of the Lux Receiver described above [include pic of circuit]
Las Receiver with GFP attached downstream ("X2")
Rhl Receiver with GFP attached downstream ("X3")
Cin Receiver with GFP attached downstream ("X4")

Of these parts, the receivers are thought to be the most generally useful as any protein coding region can theoretically be attached downstream of them to be triggered upon induction with the appropriate signal sender (or alternatively by addition of the corresponding AHL). Therefore, we tested those parts extensively, using both PCR sequencing to confirm their nucleotide sequences and fluorimetry measurements to check that the test constructs were functioning as planned.

Unfortunately, we did not have time to properly integrate our sensor modules with the color or motility modules to build the final devices that we initially devised. In future years hopefully we will be able to augment our sensors with further improvements as described under the 'Future Work' section below.

To test the effectiveness of our signaling system, we attached a GFP coding unit (rbs + protein coding region + terminator) behind each Receiver's AHL-activated promoter. (See ‘Test Constructs’ under Parts & Devices above), which theoretically enables us to measure the amount of transcription upregulation that occurs in the presence of AHL.

We then tested our parts in two ways: First we tested the Receivers indirectly by adding chemically-derived AHL into culture solutions of the Test Constructs and looking for GFP fluorescence, which will indicate that transcription has been activated. [include link to results 1]

Following that, we tested the Senders by using the above Receivers. We hoped to determine the amount of AHL produced by the Senders by comparing transcription activity induced by the Senders in relation to AHL-induced transcription activity. However this year we only managed to qualitatively determine that our Senders were functioning. We carried out this test by mixing Senders and Receivers, incubating overnight and comparing the resulting fluorescence of the mixed cultures with that of ‘negative control’ Receivers-only cultures. [include link to results 2]

For a description of the protocols followed during these tests, please refer to the Notes section: [link to protocol 1; link to protocol 2]

To make an effective simulation program we regarded the cell's movement as a diffusion. Therefore we applied the fick's second law of diffusion in modeling the cells movement along with the diffusion of autoinducers. Since the law of diffusion is a differential equation(strictly speaking a partial differential equation) the need for a numerical solution was inevitable. We used the finite difference method. We then applied it and converted it which gave us the equation (1) and (2), each representing the diffusion of the cell(or colony) and the autoinducer. These two equations were the bases of our program.

...(1)
...(2)

Finite difference method-Expicit method

As mentioned above we used the finite difference method (and also the explicit method). Finite difference methods are widely used numerical methods in solving differential equations, by considering the differential equations as a finite difference equations. We also used the explicit method which calculates the future state of a system by using the current state of the system.

Fick's second law of diffusion can be written as

・・・(3)
where C[normalized amount/μm³] is the concentration, t[s] is the time, x[μm] is the position, and D is the diffusion coefficient in dimensions of [μm²/ｓ].(Normalized　amount is an amount of cell divided by maximum.)

If we apply the finite difference method to the above equation, fick's second law of diffusion can be expressed in a finite difference equation written as

・・・(4)
where the concentration[C] of a substance at a time[t], and at a position [i,j] is represented as C^t_i,j.

In this experiment, for convenience we let Δt=h and Δx=1 but for a more accurate result Δt and Δx can be adjusted.

By applying the explicit method to the above equation again we reach the following equation.

・・・(5)
Since the concentration of the substance can not be lower than zero, every term of the equation must be over zero. So we reach the following condition.
・・・(6)

Determination of the values of parameters

For accurate results, the precise determination of the values of parameters used in the simulation was essential. We had to determine the values of the two diffusion coefficients(the cell and the autoinducer) along with the production rate of the autoinducers and the growth rate of the colony itself.

The growth rate of the colony was measured by experiment which took place in our lab[fig_model.1]. As a result the value of the colony's logarithmic　growth rate μ is 0.0024[s^-1].

fig_model.1 growth curve

The diffusion coefficient of the cell is, by definition,
・・・(7)
where v_cell, the average speed of a cell, is 0.02 mm/h(=20 μm/s) and T, the average random walk time, is 1s. So in conclusion the diffusion coefficient of the cell is 300 μm²/s [2].

The diffusion coefficient of serine is known as 1000 μm²/sec. Since the diffusion coefficient is inversely proportional to the square root of the molecular weight by simple computation we were able to figure out the diffusion coefficients of the autoinducers we used. The autoinducer's diffusion coefficients is as follows[1].

D_C4HSL = 784 μm²/s

D_3OC6HSL = 702 μm²/s

D_3OC12HSL = 607 μm²/s

The production of autoinducers and the growth of E.coli

Since the E.coli used in this experiment has itself an ability to produce autoinducers, a term which considers the cell's production of autoinducers must be contained in the equation. The added equation is shown as it follows.

・・・(8)
where δ is the production rate of the individual cell. Because the autoinducer's amount, produced by the colony increases in proportion to the numbers of individual E.colis, a term of multiplication to the cell's density was added.

The growth rate of the colony is slightly more complex. Since the nutrient of the medium is limited, E.coli's density is finite. So we expressed the colony's grwoth rate by an sigmoid curve.

・・・(9)

Triangle model

We try to simulate the triangle model. In this triangle model, three diferrent cells regulates other cells as circulation. Fig_model.2 indicate this model. Red cells regulate green cells movement.Green cells regulate blue cells movement. At the same rule, blue cells regulate red cells.

fig_model.2

Result of triangle model

         

Fig_model.3 and fig_model.4 are simulation results indicating a pattern which cells spread on a petri dish surface. Fig_model.3 indicates our triangle regulation model. Fig_model.4 is the result which cell movements aren't regulated.

Please compare fig_model.3 with fig_model.4. There is slight difference. In triangle regulation model, the cells stops swimming when receive individual signal. Subsequently, cells advance to cells stopped by signal which they send. So, the pattern is like fig_model.3. In no regulation model, The border of different color cells is accurately straight line. This simulation result corresponds with experiments. (see works.)

3-D movie of three colonies

We made the movie captured dynamic change of individual colony.

Adequacy of simulation and experiment

We wanted to get the information we cannot directly know from some experiments. How much length is needed to switch on the transcription and stop cells movement? Oppositely, if we are able to know the length, can we decide the value of some parameter such as AHL production rate? We tried to answer the question using computer simulation. But we haven’t had enough experimental results yet. Though we use some value of parameter on papars, we could not correctly decide other parameters. Simulation is just ideal. We find the fact from experiment. But simulation sometimes becomes powerful tool for helping access the fact. If we were able to make program to fit the experimental results, this approach is significant.

Sensor Function on Agar Plate

Once both Sender and Receiver function has been confirmed and characterized, we will attempt to characterize the function of the whole system in a way that relates to our intended usage. We will spot two colonies, one of Senders and one of Receivers, on a soft agar plate and determine the maximum distance that the Senders can successfully activate the Receivers. Of course, this will only work if the AHL molecules can effectively diffuse through agar medium, which is yet another thing that we have to find out. [include diagram of agar plate distance test]

Amplification

There is a possibility that the upregulation of transcription downstream of Receivers will be too weak, perhaps due to insufficient AHL produced by the Senders or poor diffusion of AHL through the soft agar medium. Therefore, we shall go back to the roots of natural quorum-sensing systems by introducing a positive feedback loop, in which weak detection of AHL will lead to transcription of genes that either enable production of even more AHL, leading to an escalating increase in AHL concentration, or more AHL receptors, leading to higher AHL sensitivity. [include sample circuit diagrams]

Quenching

There is also a possibility that the signals received may be too strong and render the system useless for producing color gradation. For example if the signal diffuses too efficiently the receiving cells may all uniformly change into the secondary color instead of showing a gradual change in color at the edges near the senders.

Therefore, we may consider implementing a signal-quenching system, perhaps using the aiiA system which involves the AHL-degrading aiiA enzyme. The aiiA gene can be made to express in limited quantity, such that the amount of intracellular AHL will be reduced but not completely erased. This causes only the receivers closer to senders, which receive a higher amount of AHL by diffusion, to be activated while the receivers further away receive too little AHL to pass the quenching threshold and are thus kept inactivated. [include sample circuit diagram]

Extensive quantitative characterization of parts will be required to determine which of these future work ideas should be pursued and to what degree. As we have unfortunately ran out of time this year, these ideas will have to wait till the next iGEM Competition.

[1]H.C.Berg, D.A.Brown, Chemtaxis in Esherichia coli analysed by three-dimensional traking, Nature 239 (5374) (1972) 500–504.

[2]H.C.Berg, Random Walks in Biology, Princeton University Press, 1983.